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ABSTRACT 

The  present  thesis  will  allow  Users  to  analyze  rectangular 
plates  with  or  without  a  rectangular  hole  with  its  edges  parallel 
to  the  edges  of  the  plate,  by  using  an  approximate  grid  analogy. 

The  method  is  completely  computerized  and  the  user  need 
only  specify  properties,  geometry,  loading  and  kind  of  output  de- 
sirable. The  computer  programs  then  will  generate  a  deck  of  cards, 
that  form  a  complete  STRUDL  input.  STRUDL  itself  is  an  equivalent 
of  STRESS  but  with  much  more  useful  capabilities. 

Initially  it  was  intended  to  make  everything  an  one  step 
process,  i.e.  to  read  injxit  in  a  free-field  form,  generate  all 
necessary  statements  of  STRUDL,  store  them  on  a  file  and  let 
STRUDL  read  the  file  and  give  the  desired  output.  However  due  to 
several  delays  in  the  Civil  Engineering  system  Laboratory,  STRUDL 
is  not  yet  available  and  CDL  (Command  definition  language)  is  now 
in  the  process  of  experimentation.  In  consequence  the  user  has  to 
supply  his  input  according  to  specific  format  and  the  linking  be- 
tween this  subsystem  and  STRUDL  is  not  automatic.  In  other  words 
the  process  is  a  two-step  one  and  the  user  has  to  take  the  output 
of  this  program  and  use  it  as  a  STRUDL  job. 

The  solution  to  the  flexural  plate  problem  contained  herein, 
is  extremely  general  regarding  boundary  conditions  and  loadings 
so  that  almost  every  engineering  problem  of  this  nature  can  be 
taken  care  of. 


Ihe  particular  connection  of  the  present  job  with  Naval 
Architecture  is  that  solutions  to  bottom  and  side  plating  with 
hydrostatic  loads »  bulkheads  with  openings  and  decks  with  car,  o 
holes  will  become  readily  available  with  a  ainiirsira  of  human  effort 
and  intervention. 
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CHAPTER  I 

-i:o^cy*,:.  ~..v  ..vX-Ffi 
The  development  of  tee  large  digital  computers  was  a  deci- 
eive  step  towards  the  solution  of  many  engineering  problems.  Indeed 
during  the  past  few  years  people  looked  back  and  reexa* 

but  it  was  impossible  to  carry  it  out  by  hand.  it.  was  *  »**«».u 
simple  matter  to  solve  all  these  problems  then  by  a  digital  com- 
puter. 

Among  all  these  other  problems  one  is  the  structural  analy- 
sis of  plates  a  special  case  of  which  is  the  rectangular  plate  in 
flexure.  Several  scientists  and  engineers  have  done  a  marvellous 
job  in  solving  special  cases  of  it. 

It  seems  that  there  would  be  two  distinct  approaches  in 
coding  the  solution  of  the  problem  into  a  computers 

(i)  Take  each  special  case  and  code  its  exact  solution. 

(ii)  Invent  a  unified  approach  applicable  in  all  or  most 
special  cases* 

The  first  approach  requires  a  tremendous  ammount  of  work. 
Besides  it  could  never  be  called  general*  For,  no  matter  how  many 
special  cases  one  solves , there  will  always  be  some  other  special 
case  left  out. 

The  second  approach  seems  quite  attractive,  with  the  only 
disadvantage  that  it  is  an  approximate  one.  However,  in  cases  where 


2. 


an  exact  solution  is  not  known,  an  approximate  one  is  best  than 
nothing.  Besides,  very  often  engineers  are  interested  in  getting 
a  first  estimate  of  what  is  happening  in  a  system* 

In  the  structural  analysis  of  rectangular  plates,  the  grid 
analogy  which  will  be  described  in  detail  in  tho  next  chapter,  is 
among  the  best  approximate  methods  to  attack  tho  general  problem. 
In  addition  since  several  methods  to  solve  skeletal  structures 
have  now  been  developped,  one  has  only  to  transform  the  plate 
into  such  a  structure,  that  existing  computer  formulations  can 
solve* 

The  goal  set  forward  by  the  author  has  been  to  develop 
a  set  of  computer  programs,  able  to  reduce  any  rectangular  plate 
to  its  skeletal  equivalent  and  then  solve  it  as  such*  Although 
the  sample  problem  included  in  the  present  work  is  a  very  simple 
one,  the  capabilities  of  the  programs  developed  are  quite  general* 
In  fact  the  choice  of  the  sa.r.ple  problem,  made  strictly  because 
its  solution  already  exists,  is  a  very  unfortunate  one  in  illustra- 
ting the  application  of  the  method  into  the  more  general  problem* 

The  case  of  a  rectangular  hole  on  the  plate  was  included 
here,  due  to  the  interest  of  Naval  Architects  on  the  effect  of 
cargo  openings  on  the  strength  of  ship*s  decks*  Also  man-holes 
on  vertical  bulkheads  forming  boundaries  of  fuell-oil  tanks  con- 
stitute another  field  of  application* 

What  really  "lakes  the  method  powerful  is  its  generality 
regarding  the  boundary  conditions  and  the  loading. 


Indeed,  point  supports  at  any  point  over  the  are*  of  the 
plate,  line  boundary  conditions  over  all  or  any  part  of  the  edgos, 
concentrated  point  loads,  uniform  or  linearly  distributed  line 
and  area  loads  or  any  combination  of  them,  are  all  cases  taat  can 
be  taken  care  by  the  programs  with  a  minimum  of  effort.  In  addition 
the  ability  to  get  any  kind  of  output  desired  might  be  vexy  useful 
in  some  cases* 

In  the  following  chapters  the  analogy  will  be  described 
and  the  method  of  its  computer  application  will  be  explained.  In 
addition  directions  will  be  given  to  the  users  on  how  to  use  the 
programs. 
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CHAFER  II 

qftp  AJ^/m 

2*1*-    General 

Perhaps  one  could  say  with  complete  justification  that  the  gene* 
rel  problem  of  a  plats  In  flexure*  has  bean  solved,  at  least  in  principle. 
Indeed  the  differential  equation  describing  the  phenomenon  is  known  for 
asny  years  and  it  is  only  necessary  to  consider  the  loading  and  boundary 
conditions  applicable  in  each  particular  problem.  However  this  last  step* 
namely,  the  solution  of  the  differential  equation  for  a  given  set  of 
boundary  conditions  and  a  given  loading*  ecmetines  presents  innumerable 
difficulties  from  the  analytical  point  of  view*  Sven  in  the  case  where 
a  solution  can  be  found  in  a  reasonable  number  of  steps*  it  is  vory 
often  expressed  in  the  form  of  an  infinite  series*  which  especially 
in  the  case  of  moments  and  afeoars  converge  rather  slowly*  a  fact  that 
reduces  their  value. 

It  is  the  above  difficulty  that  «ade  engineers  search  for  other* 
perhaps  less  accurate  but  wore  easily  developed  solutions  of  the  plate 
problem.  There  have  been  three  approaches  that  yield  approximate  solu- 
tions to  the  plate  stresses  problem. 

(i)       Finite  differences  methods. 

(11)     Orid  analogies. 

(ill)  Finite  element  methods. 

Of  these*  the  finite  differences  approach,  depending  upon  the 


boundary  conditions,  very  often  requires  fictitious  extensions  of 
asabers  or  consideration  of  members  with  fictitious  properties.  This  ma- 
kes a  rather  complicated  task  the  coding  of  the  general  rectangular  plate 
solution  to  a  digital  computer.  There  is  also  some  suspicion  that  Inif 
method  gives  results  deviating  from  the  exact  ones  (those  that  they 
exist)  more  than  those  computed  by  the  other  tvo  methods. 

The  finite  element  method  seems  to  be  quite  promising  not  only  for 
the  rectangular  plate  case  but  for  shells  of  arbitrary  shape.  Conside- 
rable number  of  people  is  working  on  that  method,  in  the  civil  Engineering 
department  of  M,I,T. ,  but  there  will  be  some  time  before  this  method  be- 
comes available  to  users. 

The  grid  analogy  was  selected  by  the  author  for  the  following 
reasons. 

(i)  It  is  rather  simple  and  gives  very  satisfactory  results 
regarding  moment  intensities  and  deflections  as  it  was  indicated  in  (i), 
(2)  and  (5). 

(ii)  It  can  be  rather  easily  programmed  to  a  digital  computer. 

(iii)  The  solution  to  the  equivalent  grid  can  be  very  easily 
found  by  the  existing  STRKSS  program  or  by  its  new  equivalent  STKOGL, 

(7), (8) ,(9). 

(iv)  The  present  work  int  nded  riaarily  for  the  field  of  Naval 
Architecture  as  well  as  other  professions  interested  in  structural 
engineering  deals  almost  exclusively  with  rectangular  flat  plates  or 
plates  with  small  curvature  w  :lch  can  be  neglected,  without  considerable 
reduction  in  accuracy. 


2.2.  Description  of  grid  analogy 

The  grid  analog  of  a  flat  plate  has  been  represented  differently 

by  different  people*  lightfoot  in  (3)  uses  a  square  grldwork  consisting 
of  grid  beams  with  flexural  rigidity  and  a  set  of  torsion  bars  with  tor- 
sional rigidity  with  no  diagonals,  Efiaba  in  (6)  in  a  program  written 
to  analyse  bridge  decks  only,  uses  TambergSs  metnod  of  analysing  beam 
and  slab  type  bridges  as  torsicnally  weak  grids  (1C). 

The  approach  taken  by  fcrennikoff  in  (2)  and  lettram  and  thisain  in 
(1)  is,  that  the  plate  can  be  represented  by  a  plane  structure  composed 
of  articulated  bars,  joined  by  pins  normal  to  the  plane  of  the  structure. 
The  pattern  has  to  repeat  itself  throughout  the  area  of  the  plate  and 
the  outline  and  external  dimensions  of  the  plate  should  be  the  same 
with  the  corresponding  ones  of  the  equivalent  structure.  Repetition  of 
the  pattern  does  not  however  mean  that  the  size  of  the  unit  has  to  re* 
main  the  same*  In  fact  by  changing  the  size  of  the  unit  we  can  satisfy 
any  exist ingdeviat ions  from  the  simplest  case  of  a  simple  rectangular 
plate,  like  differsnees  in  thickness,  holes  or  differences  in  material 
properties. 

Hrennikoff  set  up  the  problem  in  a  very  detailed  form  and 
developed  the  equations  that  will  define  the  properties  of  the  elements 
of  the  unit.  The  restrictions  that  he  imposes  in  order  to  solve  for 
the  properties  of  the  members  are,  that  the  frame  work  has  to  deflect 
properly: 

(i)  When  subjected  to  pure  bending  in  the  X  direction  with 
no  bending  deformations  in  the  T  direction. 

(ii)  When  in  pure  bending  in  the  I  direction  with  no  deforms- 
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tion  in  the  X  direction. 

(iii)  When  subjected  to  pure  torsion  on  the  X  and  Y  planes. 
These  are  general  restrictions  applicable  in  any  case  and  are 

useful  in  examining  the  suitability  of  any  particular  unit  to  be 

investigated.  The  third  condition  falls  out  whenever  the  pattern  has 

the  same  form  in  the  direction  of  both  axes. 

Subsequently  Hrenikoff  investigated  two  patterns  as  shown  in 

fig's  (2.1)  and  (2.2)  and  proceeded  to  develop  their  properties. 
Of  these  two  patterns  the  first  one  restricts  the  solution 

to  values  of  Poisson*s  ratio  of  1/3  only. 

The  second  one  although  valid  for 
any  Poisson's  ratio  has  the  disadvan- 
tage that  is  a  square  one.  This  causes 
a  serious  problem  in  the  case  that 
we  cannot  satisfy  the  boundary  don- 
ditions  with  a  square  pattern.  We 
also  notice  that  for  each  square  we 
have  10  slender  members,  which  would 
make  a  hand-solution  more  difficult 
or  a  computer  solution  longer. 
A  rectangular  pattern  with  diagonals 
only  is  much  more  general  and  simple. 
This  type  of  pattern  has  been  very 
recently  investigated  by  Yettram 
and  Husain  in  (1).  Their  pattern 

fo.S.g  Squty^e   Putters 
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looks  to  bo  the  best  of  the  existing  ones  because  it  is  simple,  easily 
developed  and  co:isiderably  general*  Due  to  the  fact  that  the  method 

is  relatively  new  and  espesially  useful  in  other  kinds'  of  .roblems 

ose  satiating  the  Langrangian  equation  like  thermal  stresses  in 
a  plate  ot  special  cases  of  fluid  flow),  it  is  believed  that  a  rela- 
tively detailed  discussion  ot  it  must  be  given.  This  will  be  dome  in 
the  next  chapter. 
2*3*     SaflU:iguiar  jjattem  grid  analogy 

The  idea  of  the  ^odel  is  that  we  can  substitute  a  rectangular 
network  consisting  of  slender  members  only,  in  place  of  a  particular 
plate  element*  Fig*  (2,3a)  shows  a  roctangula?  element  of  plating  and 
Fig* (2.3b)  its  equivalent  grid  element*  The  unknown  properties  of  the 
grid  element  members  are  the  flexural  rigidity  of  the  side  beans  Ie 
and  Is  ,  the  torsional  rigidity  of  the  side  beams     &   *       and 


and  the  flexural  rigidity  of  the  diagonals  Z^  •  That  gives  5  unknown 
properties  i.e*  Ia»  I*  •  Id  ,    &?*     and        -^ —  - 

E  E. 

The  equations  giving  the  above  5  properties  can  be  developed 
by  equating  the  rotations  at  the  nodes  of  the  grid  element  with  those 
of  the  plating  element  of  equal  else*  Although  the  method  doesn't  give 
the  cross-sectional  area  of  the  beams,  it  gives  all  tae  necessary  in- 
formation to  make  the  solution  of  the  tfrid  framework  possible*  Of 
course  the  plate  a  xi  the  grid  elements  nust  be  subject  to  statically 
equivalent  torques  and  moments* 

Considering  first  the  plating  element  we  examine  three  cases* 
(1)  i^qual  and  opposite  bending  moment  intensities  on  the  x  sides 
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(id)  idqual  and  opposite  inlensikes  on  y  sides  . 

(iii)  Self-equilibrating  torque  intensities  along  4  edges. 

Fig's  (3d),  (3©)  «»d  (3g)  show  clearly  the  WHMMfc  and  torque 
distributions  and  the  corresponding  angles.  The  values  of  the  angles 
in  terms  of  the  moment  intensities  can  be  found  by  differentiating  the 
expressions  for  the  deflections,  (11).  This  was  done  in  (1)  where 
they  came  up  with  the  following  expressions  for  the  angles. 


•■"?       *3",'    ;  U-i) 


£(^/.<t) 


®»-      777T7- T  *2*3' 


-.^  «"> 


»»-  ^"V/?  (2.5) 


For  the  grid  element  loaded  as  shown  in  Fig's  (3d),   (3**)  and 
(3fa).  Ref  (1)  gives   the  flexibility  matrix  (shown  here  in  Fig  2.  4) 
and  concludes  that  only  2  out  of  the  12  equations  are  independent  when 
we  have  only  side  bending  moments  or  self •equilibrating  torques, 
§y  selecting  properly  these  2  equations  in  each  case  we  come  up  with 
the  following  expressions  for  the  grid  angles. 
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By  equating  the  rotations  for  the  plate  and  tho  grid  elements 
wo  can  got  tho  following  equations  giving  the  required  properties  of 
tho  grid  suaaber.   [The  motaod  i*  explained  in  >nore  dotail  in  (1)J . 


wJN3 £-  (e.i2) 


9.0-rS    '      ^ 


1 


(2.13) 


E      "       9.ci-^)'    »«.  (2»i5) 


1*. 


Where    k    is  the  ratio  of  the  side  beams 
|X   is  tii*  Foiseon'e  ratio 
h     is  the  thickness  of  plating 

is  the  ioung*s  atiotkius  of  elasticity 
G     is  the  taoduius  oi  rl&kiit$ 
r    is  tae  diagonal  to  side  beam  ratio 
For  future  use  we  will  transform  the  above  relations  by  usiagt 


*s   — * (2.17) 


into  the  following  forat 


ajK^HJj    .  (2.18) 
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l^f.J^L  (2"20) 


lc=   iHi^)0"^  (2#21) 


^--^H^^C-^)  (2<22) 


1, 


viWexe  $=     W 


is  *  cowrtant  for  a  given  isotropic  plate  of  constant  thickness. 

The  potentialities  of  the  above  equations  are  immediately 
apparent,  they  can  be  used  for  the  solution  of  any  plating  problem, 
with  any  boundary  conditions  and  load  and  any  number  of  holes  provided 
that  we  can  form  a  grid  fine  enough  so  that  it  will  represent  a  parti- 
cular plating  element,  i.e.  provided  that  we  will  satisfy  equations 
(2.18)  through  (2.23)  at  each  particular  element  in  the  division  of 
the  plate.  Also  by  forming  a  proper  grid  we  can  take  care  of  non- 
isotropic  plates. 
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In  the  following  chapter  some  limitations  and  principles  of 
application  of  the  above  ssetnod  will  bo  discussed  bafore  the  computer 

approach  is  explained* 

Mi  ifc4t*t*ffn»  i  mkaMm  ?t  mtMmtim  9?  yr\4  MUM 

axpresslcns  (2*16)  through  (2,23)  ssight  give,  under  special 
circumstances  negative  values  for  the  l*s  or  J*  a.  Altuou^b  Ref.(l) 
claias  that  this  is  :»t  a  problea,  existing  computer  programs  for  the 
solution  of  a  grid  expect  positive  raojaents  of  inertia.  Xhis  limits 
the  range  of  applicability  of  the  a&ethod  in  socae  mannar. 

Ihe  range  of  the  Polsson9s  ratio  for  engineering  purposes 
varries  between    .15  and     .3    •  looking  carefully  into  equations 
(2.18)  through  (2.13)  we  conclude  that  in  order  to  assure  positive 
values  of  all  J*s  and  I*s  tne  value  of  X    *ust  be: 


£  o-c      ^   -  .-5 


U.24) 


5"^  ^   Vn£    1-81  (v#25) 


For  higher  values  of  u  than  1/3.  negative  properties  are 
unavoidable,  however  the  authors  of  Ref.(l)  state  that  for  good  re- 
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suits  the  value  of  k  must  be  between  2  and  0,5* 

The  computer  programs  will  be  such  that  will  automatically 
form  a  grid  with  k*a  (if  more  than  one)  as  much  close  to  unity  as 
possible*  It  is  obvious  here  that  the  best  results  are  obtained  for 
k»l  but  sometimes  the  boundaries  are  such  that  it  is  quite  difficult 
to  have  a  square  pattern. 

In  the  actual  application  of  the  method  the  following  points 
must  be  taken  into  account  . 

(1)  The  properties  found  above  give  the  contributions  of  only 
one  element  in  the  division  of  the  plate*  Thus  for  the  side  beam  we 
have  also  to  add  the  properties  of  the  neighbouring  elements*  Of 
course  this  is  not  true  for  the  diagonals  and  the  members  lying  along 
the  edges  of  the  plating  or  along  the  edges  of  any  existing  ooles. 

(ii)  The  boundary  conditions  along  the  edges  of  the  plating 
will  be  considered  as  applicable  on  the  nodes  of  the  grid  at  that 
neighbourhood  * 

(ill)  Any  isolated  point  supports  in  the  interior  of  the  plate 
are  considered  applicable  on  the  nearest  node* 

(d)  The  loading  is  considered  as  acting  on  the  nodes  or  on 
the  members  of  the  grid* 

The  exact  manner  in  which  the  loadings  and  the  boundary  condi- 
tions are  treated  will  be  explained  in  the  description  of  the  parti* 
eular  SUBROJTIftES  that  handle  them*  Of  course  most  of  the  assumptions 
made  here  require  that  we  have  a  fine  enough  grid  so  that  the  effect 
on  changing  the  problem  from  the  original  one  to  another  problem  will 
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be  minimized. 

Having  now  developed  the  principles  on  which  the  program  will 

operate  upon*  we  are  ready  to  embark  in  the  actual  description  of  the 

program. 
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CHAPTER  III 

A  PRIEF  DISCUSSION  OF  ICES 

3a.    overall 

Although  this  chapter  seems  fairly  obvious  to  people  in  the 
Civil  Engineering  Department  of  M.I.T.,  related  to  computer  work, 
it  is  strongly  believed  that  at  least  a  brief  description  of  ICiSS  must 
be  given  here  for  people  not  fa.riliar  with  it. 

1CSS  (Integrated  Civil  Engineering  sy stent)  is  a  computer  system 
designed  for  the  efficient  solution  of  a  variety  of  Civil  Engineering 
problems  among  which  one  is  the  structural  Analysis  Problem. 

Although  it  looks  as  the  system  is  applicable  only  to  Civil 
Engineering  problems  the  generality  of  some  SUBSYSTEMS  of  it  makes 
them  quite  useful  to  any  other  field  of  Engineering, 

The  system  has  been  espesially  designed  in  such  a  way  that  the 
user  (A  user  is  a  person  that  uses  programs  written  by  other  people 
called  the  programmers,  for  engineering  computational  work)  communi- 
cates with  the  computer  by  informing  the  computer  on  t:ie  nature  of 
the  problem  and  recommending  the  steps  to  be  taken  for  the  solution 
of  the  problem.  In  addition  he  doesn't  have  to  provide  the  input  at 
any  particular  format.  In  that  way  the  engineer  instructs  the  computer 
instead  of  following  instructions  imposed  by  the  machine. 

ICSS  has  been  made  to  fit  into  the  j€0  &&  machine  and  uses 
as  a  programming  language,  ICgTRAK.  ICETRAM  (ICES,  FCRTRAJi)  is  a  com- 
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bination  of  FORTRAH  IV  plus  powerful  new  capabilities  useful  in 
performing  engineering  work.  All  ICES  programs  are  supposed  to  be 
written  in  IC8TRAS.  For  the  purposes  of  this  work  it  suffices  to  say 
that  a  program  in  I  :£TRA2i  is  precompiled  by  a  special  precompiler 
(the  ICSTRAN  precompiler)  into  a  FORTRAN  program  which  is  finally 
processed  by  the  FORTRAN  compiler, 

A  SUBSYSTEM  in  ICES  is  a  set  of  programs  that  can  be  used 
independently  or  collectively  to  solve  engineering  problems.  One  of 
the  subsystems  is  STRUDL  (STRUCTURAL  ANALYSIS  *  DESIGN)  which  includes 
all  the  capabilities  of  STiiSSS  plus  new  ones.  In  particular  the  work 
done  herein  will  hopefully  initiate  another  subsystem  called  PLATSS, 
which  eventually  will  be  able  to  solve  all  kinds  of  problems  related 
to  plates, 

3.2.     Seme  capably, j^ea  of  WmW 

There  have  been  several  capabilities  of  ICES  and  1CETRAN  that 
simplify  very  much  this  work  and  in  addition  make  It  quite  general. 
These  capabilites  are: 

(a)  Interaction  between  subsystems.  Since  the  approach  to  this 
problem  has  been  to  reduce  it  into  an  equivalent  grid  and  then  let 
SSSHDL  solve  it,  the  interaction  between  subsystems  is  very  important. 
At  the  time  that  this  work  was  untertaken  this  feature  was  not  avai- 
lable and  that  is  wny  it  was  decided  to  make  the  output  of  tela 
program  a  deck  of  cards  that  form  a  complete  STftUCL  input.  However 
it  is  a  simple  matter  to  modify  this  program  to  assure  one-step 
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processing  when  interaction  between  subsystems  bocoraes  available, 

(b)  The  XCKS  executive  ;-rorra?q.  it  is  the  ICSS  executive  system 

that  gives  the  programmer  the  capability  of  allowing  the  user  to  write 
his  input  in  an  iinglish-like  form.  The  executive  process  consists  of 
reading  the  input,  finding  the  command  in  a  dictionary,  storing  the 
data  in  the  proper  place  and  executing  associated  subroutines. 

The  rules  for  defining  commands  and  associated  procedures  are 
handled  by  the  command  Definition  Language  (CDL)  compiler.  More  in- 
formation about  how  to  define  commands,  what  the  IC£X(ICES  EXECUTIVE) 
input  structure  looks  like  etc.,  may  be  found  in  (9)»  chapter  3, 

(c)  ifrnamle  memory  allocation 

The  characteristics  of  the  equivalent  Grid  wight  change  from 
one  problem  to  another.  Wo  cannot  then  say  that  we  will  arbitrarily 
set  the  number  cf  nodes  to  500  for  examble,  since  t:  is  will  depend  on 
the  kind  of  the  problem  in  general  and  tr.e  accuracy  that  the  user 
desires.  That  indicates  that  the  arrays  to  be  used  throughout  the 
problem  are  not  of  a  given  size  but  their  sise  varries  from  problem 
to  problem.  As  it  is  well-known  in  FORTRAN  one  would  have  to  indicate 
the  sise  of  the  different  arrays  at  the  beginning  of  the  program. 
Then  if  the  problem  is  small  (relative  to  the  defined  sise  of  arrays) 
memory  storage  is  waisted.  If  the  problem  is  large  it  cannot  be  exe- 
cuted. 

3f  using  the  dynamic  memory  allocation  in  ICSTRAN  the  : ro;  rammer 
can  define  the  sixe  of  the  arrays  during  the  execution  time.  This  in 
turn  has  the  effect  that  at  no  time  the  memory  is  waisted. 

In  addition  to  bercg  able  to  define  the  sise  of  the  arrays 
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during  the  execution  period  smother  advantage  of  the  dynamic  memory 
allocation  is  tho  fact  that  secondary  memory  it  considered  a  logical 
extension  of  the  primary  memory.  Thus  at  any  tine  that  the  primary 
WBiory  is  filled  up,  the  secondary  memory  automatically  is  uead  to 
store  additional  information. 

xhe  important  thing  then  for  the  user  is  that  the  sise  of 
his  program  is  not  restricted  by  the  primary  memory  siae.  tfhat  makes 
the  dynamic  memory  allocation  powerful  and  nice  is  the  fact  that 
neither  the  programmer  nor  the  user  has  to  worry  about  siae  of  arrays, 
provided  that  they  are  dynamic  and  their  sise  has  been  defined  by 
an  arithmetic  expression. 

(d)     OrfifrnlfaUfrft  o£  W^^W  MOT 

The  last  extremely  useful  feature  of  ICETRA.K,  regarding  this 

work  is  the  organisation  of  secondary  memory.  Since  the  data  generated 

by  the  PLATE  program  will  be  finally  used  by  STRODL  it  looks  quite 

right  to  write  all  of  it  on  a  file  and  then  let  STRUDL  read  that 

file  and  proceed  to  the  solution. 

this  has  not  been  done  at  this  stage  because  of  the  following 
reasons. 

(i)  When  this  work  was  done  STRUM,  didn't  have  the  capability 
to  accept  information  from  a  file  in  an  internal  form, 

(ii)  The  system  Ri&TUJSS  that  can  be  used  to  write  information 
on  a  file  cannot  do  it  under  a  specified  FORTRAN  type  FGRHaT,  a  fact 
that  would  allow  the  writing  of  inf  orraation  on  a  FILE  in  an  external 
form. 
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(ill)  Thare  are  statements  in  FOHTRAil  that  can  bo  used  to  write 
information  on  a  FILE  under  a  3,  ecif ied  it  statment  but  the  siee 

of  the  file  -nust  be  specified  in  advance  (Like  the  DLMSNSIGS  statement) 
which  makes  them  inadequate  for  our  case* 

However  when  the  capabilities  (a)  and  (b)  will  become  aifiailable 
it  will  be  a  simple  matter  to  modify  the  programs  in  such  a  way  that 
information  will  be  directly  written  on  a  file  wnich  S7P.U  0L  will  be 
able  to  read. 

This  concludes  the  brief  discussion  of  ICES  given  in  tiis 
chapter.  The  wir^-ose  of  it,  was  to  c-ive  to  the  non-familiar  reader 
a  slight  idea,  why  this  work  has  been  done  in  ICETftLM.  A  complete 
description  of  ICSS  is  given  in  (9)  . 
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CHAPTER  IV 


CuMPUTER  PROGRAMS 


4.1.  Introduction  and  capabilities 

Before  a  general  description  of  the  computer  programs  is  given, 
it  is  necessary  to  give  a  rather  detailed  description  of  the  capabi- 
lities of  the  new  subsystem. 
>  4 


FIG. 4.  I      -    LOClVTtOM         UNO       ORHNTMIOH        Of    ftX  E-S 

Fig.    (4.1)   indicates  the  orientation  and  location  of  the  axes 
of  coordinates  relative  to  the  plate.   The    I    -  axis  is  positive  upwards. 
This  subsystem  will  be  able  to  handle  constant  thickness  rectangular 
plates  of  constant  Young's  modulus  and  Poisson's  ratio.   The  capabilities 
of  the  subsystem  are: 

(i)   Number  of  holes:   There  can  be  either  no  hole  at  all  or  a 
rectangular  hole  with  its  edges  parallel  to  the  edges  of  the  plate,  lo- 
cated anywhere  on  the  surface  of  the  plate. 

(ii)   Boundary  conditions:   The  boundary  conditions  that  this  work 
can  handle  can  be  divided  into  two  categories.   The  first  category 
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includes  line  boundary  conditions  (nin^Qd.free  or  fixed)  arround  all 
or  any     :on  of  the  edges  cf  the  gating,  iieodloss  to  say  that  any 
one  of  the  four  edges  cf  the  pitting  night  be  hinged  along  one  portion 
cf  it,  free  along  another  and  fixed  over  the  rest  of  it.  The  second  cate- 
gory includes  simple  point  supports  over  the  entire  area  of  the 
plate  except  t  e  aole  if  there  is  one.  So  line  supports  are  allowed  at 
any  other  locations  except  Um  four  edges,  however  any  line  support, 
in  a  position  other  than  the  edges,  can  be  perfectly  well  treated  as 
a  set  of  equally  spaced  roint  supports, 

(iil)  Loading:  The  loading  can  be  of  one  of  three  types.  A 
vertical  force  along  the  *  -  axis  or  a  moment  the  vector  of  which 
is  parallel  to  trie  x  or  y  axis.  If  the  vector  of  the  moment  has  an 
orientation  other  than  parallel  to  the  x  or  y  axis  it  can  be  treated 
by  considering  its  components.  In  any  of  the  above  throe  cases  the  loa- 
ding can  be  a  concentrated  one  acting  on  a  point,  a  uniform  one  along 
a  straight  line  or  a  linear  one  along  a  straight  lino.  Forces  not 
parallel  to  the  i-  axis  or  bending  moments  having  a  component  parallel 
to  the  i   •  axis  connot  be  treated  since  this  requires  a  different 
grid  analogy, 

(*v)  U>t.xiti  Moments  forces  or  reactions  at  specified  toints, 
along  straight  lines,  or  over  the  entire  area  of  the  plating  can  be 
requested, 

A  more  detailed  description  explaining  how  each  case  has  been 
taken  care  of,  will  be  given  in  the  documentation  chapter  of  the  programs 
that  perform  specific  tasks. 
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M«  A  macroscopic  view  of  the  program 

'ihe  entire  work  can  be  logically  divided  into  4  steps,  namely: 

(i)  Formation  of  the  equivalent  grid  . 

(ii)  Reduction  of  the  boundary  conditions  into  equivalent 
grid  supports. 

(iii)  Reduction  of  the  loads  into  equivalent  grid  loads  acting 
on  the  joints  or  the  members  of  the  grid. 

(iv)  Reduction  of  the  specified  output  into  an  equivalent 
legitimate  grid  output. 

Each  of  the  U  steps  outlined  above  is  performed  by  a  set  of 
subroutines*  Some  of  the  subroutines  are  equally  applicable  in  both 
cases  where  a  hole  does  or  doesnot  exist.  Whenever  tills  is  not 
possible  a  different  subroutine  is  used  for  each  case.  The  USER 
however  doesnot  have  to  worry  about  these  details.  If  he  specifies 
that  there  or  there  is  not  a  hole  the  proper  subroutines  are  exe- 
cuted automatically. 

In  the  next  chapters  a  macroscopic  description  of  each  of 
the  above  ^  steps  will  be  given. 

4.3-  Formation  of  the  «rrid 

The  first  step  in  the  solution  of  the  problem  is  the  gene- 
ration of  the  equivalent  grid.  This  step  is  a  combination  of  geo- 
metry and  particular  boundary  conditions  in  order  to  produce  the 
proper  JOINT  COORDINATES,  JOINT  RELEASES  and  MfsMBER  INCIDENCES 
commands  in  SFRUDL.  This  part  of  the  job  is  handled  by  the  following 
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subroutine*  i 

Dli  GHC     (J^vlsion  joints  -  Hole) 

DI?Om     (Division  j&ints  -  go  iyiole) 

GRID     (Geometric  Formation  of  Grid  -  no  hole) 

MUGRID     (Geometric  Foraation  of  l££lti  Grid  -  there  is  a  hole) 

HEXKNH     (Meaber  Incidences  -  Ko  hole) 

MalKHO     (£gnber  J&cidences  -  M*-«) 

KKMPRO     (Member  properties) 

JOINRE     (Joint  pleases) 

A  generalized  flow-dia  :ra»  is  given  in  Fig^.^l.  explaining 
the  sequence  and  interaction  of  t  is     phase  as  well  as  the  other 
three    phases. 

*•*•    Bouqfory  OTtoUopt 

The  boundary  conditions  of  the  plate  are  transformed  into 
an  equivalent  set  of  supports  for  the  grid  by  the  sase  subroutines 
no  natter  if  there  is  or  there  is  not  a  hole. 

the  specific  routines , that  do  the  job  here, are* 

FOBOC  (jPo.int  b&undary  Conditions) 

LIBGC  (L^ne  boundary  Conditions) 

4«5.  Loading 

There  are  four  subroutines  that  perform  the  transformation 
of  the  plate  loading  into  an  equivalent  grid  loading.  These  are: 

UN.oLo  (Line  Loadings) 
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..sLG     ( Point  j^oadings) 
LIKE     (A  logical  extension  of  L1N2LQ  to  get  arround  the  difficulty 
of  tin  I1£.E  to  treat  large  number  of  symbols) 

■^     (Area  fading) 

There  is  only  one  subroutine  (OUT)  which  in  conduction  with 
I  -..18L0  and  LI^SLO  is  taking  care  of  the  output. 


*9. 


CHAPTER     V 

HOOKAH     DOCUK&vIATIQN 
5.1.     General 

This  chapter  contains  adequate  information  on  how  the  c\ li- 
ferent programs  function  and  how  they  are  linked  to  each  other. 
This  information  is  necessary  especially  for  other  programmers 
intending  to  modify  a  particular  subroutine.  In  fact  this  chapter 
is  mostly     useful  to  other  subprogrammers  but  it  can  also  be  of 
some  help  to  the  users  of  the  program. 

Especially  for  users  the  chapter  on  a  sample  problem  and 
form  of  data  is  extremely  useful ♦  in  guiding  theft  how  to  write 
their  own  input. 

In  writing  the  documentation  of  the  KOBOOTZm  the  author 
conformed  with  the  recommended  methods  of  (9)«  Particular  attention 
has  been  paid  on  including  all  the  fine  points  of  the  programs 
into  either  their  long  description  or  their  flow  diagrams,  however 
if  any  point  is  not  made  clear  in  the  documentation  of  a  program 
the  reader  is  referred  to  Appendix  8  which  contains  all  the  ICETRAN 
listings. 

Some  general  points  applicable  to  the  entire  documentation 
program  are: 

(i)  All  COMMONS  are  the  same  throughout  all  flmOBTZla*, 

(ii)  All  global  dynamic  ARflAUS  are  defined  in  either  DIPQMH 
or  DZraO  and  are  destroyed  by  the  MAIN  PRXKAM. 
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(ill)  All  LOCAL  DXHAftXC  ARRAYS  are  defined  in  the  particu- 
lar SUBROUTINE  that  uses  them  and  destroyed  before  a  return  is 
executed  from  the  subroutine, 

(iv)  All  transmitted  dyn«  arrays  are  defined  in  the  CALLING 
program  and  destroyed  before  a  return  is  executed  from  the  calling 
program. 

(v)  Since  it  is  impossible  at  this  time  to  have  I/O  state- 
ments with  dynamic  array  elements,  all  output  is  written  by  first 
storing  the  dynamic  element  into  a  non-dynamic  one  and  then  using 
I/O  statements  with  this  last  clement. 

In  tho  following  chapters  a  description  for  each  subrou- 
tine is  given.  Appendix  A  contains  the  flow  diagrams,  Appendix  B 
the  Icetran  listings  and  Appendix  S  a  dictionary  of  all  names  in 
COMMON. 
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1  -  NAME  -  Subroutine  BXFOJH 

2  •  AOtSOi  -  John  KollinUtis 

3  -  BATE  -Hay  7.  1966 
Jf  -  SCn  -  Specific 

5  -  S;  CRT  DESCRIPTION  -  This  subroutine  is  used  in  the  case  that 

there  is  no  hole,  in  order  to  compute  the  nuraber  of  division 
points  along  each  side  for  a  specified  naxttnun  of  joints. 

6  -  LMmhGK  m  Icetran 

7  -  FLOW  CHART  -  As  shown  in  Fig  (A5.2)  Appendix  A* 

8  m  USAGE  -  CALL  DIPCMh  with  the  necessary  input  values  in  CQ&KUI, 

Input  variables  are: 

VSLEN  *  mm  •  Span  of  plate 

HOLES  ■  HGLEtt  «  Width  of  plate 

HAf&JJO  «  Desired  ssaxisaxm  number  of  joints. 

9  •  OUTHJT  -  Output  in  COHKON  includes: 

MBhOJl  «  timber  of  division  points  along  the  x«&xis 
Mm  »  Number  of  division  points  along  the  y-axis 

aj         "\ 
MMNI 

With  the  same  meaning 
BUEIDI  as  in  subroutine  DIPCSiO 

7 
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10-  U 

Source  Program  -     60  Ieetran  instructions 
.ect  Progress  -     6?  Binary  cards 

11-  CALLS  -  Nona 

ia.  ions  -  aoasiiBH  futss 

13  ~   .  ,  DSSCRIi-TIOS  -  This  subroutine  computes  the  division  pts 
along  either  side  of  tho  plate  in  such  a  way  that  the  spacing 
of  thee?  will  be  equal.  Since  tho  values  found  will  be  fractio- 
nal it  truncates  the  values  computed  to  the  lower  integral. 
Subsequently  it  axar&nes  whether  it  is  possible  to  Increase 
the  division  points  along  only  on©  side  so  that  we  approach 
more  to  the  specified  number  of  division  pts. 
In  addition,  since  this  subroutine  is  the  first  one  to  be 
executed  if  there  is  no  hole  it  defines  all  the  global  dyna- 
mic arrays. 
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5.3.  SUBROUTINE  DIPOHO 

1  -  NAME  -  Subroutine  DIPOHO 

2  -  AUTHOR  -  John  Kolliniatis 

3  -  DATE  -  April  20,  1966 

4  -  TYPE  -  Specific 

5  -  SrfORT  DESCRIPTION  -  This  subroutine  is  used  in  the  case  of  a 

rectangular  plate  with  a  hole  and  a  specified  maximum  number 
of  joints  for  the  equivalent  grid,  in  order  to  find  the  sub- 
division points  for  each  of  the  eight  sub-panels  that  the 
plate  is  considered  to  be  consisting  of, 

6  -  LANGUAGE  -  Icetran 

7  -  FLOW  CriART  -  As  shown  in  Fig.   (A5.3)  Appendix  A. 

8  -  USAGE  -  Call  DIPOHO,  with  the  necessary  input  values  in  COMMON, 

Input  variables  are: 

HOLEN  m  Length  of  plate  along  x-axis 

VELEN  ■  Width  of  plate  along  y-axis 

HOLEN1  •  Distance  along  the  x-axis  between  the  center  of  co- 
ordinates and  the  lower  left  corner  of  the  hole. 

VELEN1  ■  Distance  along  the  y-axis  between  the  center  of 

coordinates  and  the  lower  left  corner  of  the  hole 

H0LEN2  *  Length  of  hole  along  the  x-axis. 

VELEN2  ■  Width  of  hole  along  the  y-axis 

MANUJO  ■  Desired  maximum  number  of  joints  for  the  equivalent 
grid 
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9  -  OUTPUT  -  Output  in  COMMON  includes: 

NUHOJ1  :  Number  of  division  points  along  the  x-axis  of  sub- 
panels  A,  D  and  F  (see  Fig.  5«1  and  5»2). 

NUH0J2  :  Number  of  division  points  along  the  x-axis  of  sub- 
panels  B  and  G. 

NUH0J3  :  Number  of  division  points  along  the  x-asis  of  sub- 
panels  C,  E  and  H. 

NUVEJ1  :  Number  of  division  points  along  the  y-axis  of  sub- 
panels  A,B  and  C. 

NUVEJ2  :  Number  of  division  points  along  the  y-axis  of  sub- 
panels  D  and  E. 
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NUVEJ3  :  Number  of  division  points  along  the  y-axis  of  sub- 
panels  F,  G  and  H. 

KJ  :  Number  of  joints. 

NUHOMB  :  Number  of  horizontal  members 

NUVEMB  :  Number  of  vertical  members 

NURIDI  :  Number  of  right  -  diagonals 

NULEDI  :  Number  of  left  -  diagonals 

NUME  :  Total  number  of  members 


'   ^*^ 

rt 

\ 

NUH0J1  =  h 

NUH0J2  =  4 
KUHQJ3  =  2 
NUVSJ1  =  4 
NUVEJ2  =  3 
NUVEJ3  =  3 
NUHOMB  =  53 
NUVEMB  =  52 
NURIDI  =  43 
NULEDI  =  43 
NUKE   =  191 
NJ    =  62 
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10  -Li 

Source  Program  -  70  ICEFSUL&  instructions 
Object  Progress  -  SO  Binary  cards 

11-  CALLS  -  Son© 

12-  APPLICATIONS  -  SOBSTSTES  PLATES 

13-  LONG  BMQMH  -  Initially  this  subroutine  exas&nes  whether 
it  is  possible  to  select  values  for  the  division  points  of  the 
sub-panels  in  sucto  a  way  that  any  grid  element  is  a  square  one. 
Since  these  results  are  in  fractional  fore  in  general  a  trunca- 
tion is  made  to  the  lowest  integral  value*  In  ease  the  hole  is 
very  close  to  one  of  the  sides  of  the  plate,  the  pregraa?  over- 
rides the  restriction  imposed  by  the  V$M  by  his  specification 
of  KAJ2IJ0,  and  arranges  the  division  points  la  such  a  way  that 
ensures  consideration  of  the  hole  effect*  Of  course  in  this 
case  II  is  larger  than  HAi&VO, 

In  addition  since  this  subroutine  is  the  first  one  to  be  ex- 
ecuted in  the  case  of  a  plate  with  no  holes  it  is  used  to  de- 
fine all  global  dynamic  arrays. 
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5.a.      :■  •■?!.•       ■  .:, 

1  -  -  Subroutine  GRID 

2  -  ft  -  John  Kolliniatis 

3  -  OATS  -  :iay  ?•  19^6 
k  -  Ti?S  -  Specific 

5  -  MOM  .    rSSM  *  This  subroutine  is  executed  after  BSPM 

in  order  to  form  a  list  of  the  joints  and  thoir  coordinates 
for  later  use  (Tnsre  is  no  hole), 

6  -  LASGCAGE  -  Icotran 

7  -  FLOW  -  As  shown  in  Fig,   (A5.5)  Appendix  A, 

8  -  USAGE  -  CALL  GRID  with  the  following  input  variables  in  WMBWi 

■I 

KBiCJl 

S6H0J2 

IB 

VSLEK 

9  -  OUTPUT  -  Output  in  CGJttON  includes* 

OOH(I)  ■  x-coordinate  of  Joint  I 

VSCOR(I)  •  y«coordinate  of  joint  I 

NUJG(I)     *  fAusber  of  joint  1*1 
XQ-  LSNGfc    - 

Source  Program  -         ^3       Ieetran  instructions 

Object  Progress  -         30       Binary  cards 
11-  CALLS  ~  Hone 


with  the  same  meaning 
as  in  DIPONH. 


33. 


12-  -  n       PLATES 

13-  JO  SCftXFfSQN  •  All  that  It  has  to  be  said  her©  Is  that 

the  nuieberijig  of  the  joints  proceeds  from  the  bottom  to  the 
top  and  from  the  left  to  the  right.  Although  the  flow  diagram 
Is  the  saree  as  In  MKA2S  the  procedure  followed  differs  appre* 
ciably. 
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5*5.     SilLH-UTINS  MUGR1D 

1  -  NAME  -  Subroutine  mQRlD 

2  -  AUTHOR  •  John  Kolliniatis 

3  -  DATE  -  April  20,  1966 

4  -  B91  -  Specific 

5  -  SHORT  DESCRIPTION  -    This  subroutine  is  used  in  the  case  of  a 

rectangular  plate  with  a  hole  in  order  to  number  the  joints  of 
the  equivalent  grid  and  find  their  horizontal  and  vertical 
coordinates. 

6  -  LANGUAGE  ~  Icetran 

7  m  FLOW  CHART  -  As  shown  in  Fig.   (A^5) ,  Appendix  A. 

6  -  USAGE  -  CALX.  HOGRID.  with  the  necessary  input  values  in  COMMON 
Input  variables  are: 
HUB 
VELENI 
hOLE82 

mm 
mm 
mon 
mm 
mai) 

NUVEJ3 


having  the  sane 
meaning  as  in 
Subroutine  DIPOUO 


H0LKN3  ■  iiOLEN  -  iiOLEMI  -  hOLOB 


*o. 


VSLi^3  ■  nUM  -  VHLENI  -  VKLEN2 
9  -  GOTKiT  -  Output  in  COMMON  includes  four  dynamic  arrays  of  siae  NJ 
where  U  being  in  OQMKM  itself,  has  the  same  meaning  as  in 
DIFOHC*  These  are: 
mJQ(l)=  Number  of  Joint  I. 
hGCOR(I)  «  Horizontal  coordinate  of  joint  I. 
VSCOR(I)  ■  Vertical  coordinate  of  joint  I. 
JtRSMJ(I)  -  Release  number  of  joint  I. 

10-  LENGTH 

Source  Program  -107  ICETRAN  instructions 
Object  rrogram  -  ^5    Machine  instructions 

11-  CALLS  -  hone 

12-  AH-LlCATIOtiS  -  SUBSYSTEM  PLATES 

13-  LQI&  ^CRXfTIOS 

This  subroutine  which  is  called  by  the  Main  program  after 
subroutine  BXPCttQ  has  been  executed  is  used  in  order  to  find  all  the 
nodes  of  the  equivalent  grid  and  their  coordinates  relative  to  a 
system  of  coordinates  located  at  the  left  lower  corner  of  the  plate, 
the  numbering  of  the  joints  is  as  shown  in  Fig*   (  C»i      )  i.e.  it 
proceeds  from  the  bottom  to  the  top  and  from  the  left  to  the  right. 

In  addition  if  fills  with  aero»s  the  elements  of  the  dyna- 
mic array  JGRENU  walch  is  supposed  to  be  filled  later  by  aether 
subroutine  by  characteristic  numbers  indicating  the  kind  of  the  joint. 
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5.6.  Z    Mffiflgff 

1  -  mm  -  Subroutine 


2  -  AUTHOR  -  John  Kolliniatis 

3  -  DATE  •May  7,  1966 
k  -  TIPS  •  Specific 

5  -  ShORT  DESCRIPTION  -  This  subroutine  is  called  after  DIKafc 

and  GRID  have  been  called  in  order  to  supply  to  other  subrou- 
tines called  later  information  about  the  number  of  ary  member 
and  its  minus  and  plus  Joint 

6  •  LANGUAGE  -  Icetran 

7-FLOW  CHART  -  As  shown  in  Fig  (A5.7)  Appendi*  A, 

8  •  USAGE  •  CALL  8ELNNH  with  the  following  variables  and  arrays 

(defined  in  previous  subroutines). 
HGCOR  W 

VECCR  MJiiGKB 

i«jo  mrvEHB 

Hi  NUREDI 

HULEDI 

9  •  OUTPUT  -  In  COKMON  includes  the  following! 

flUMSKB(I)  a  Number  of  member  1  =  1 
WMIJO(I)  »  fJumber  of  alms  joint  of  member  I 
i»PUO(I)  a>      •        •     plus        •      •        •        ■ 
KBiOJO(l)  ■  kind  of  member  I 

KimJO(I)  «  o    side  beam 

KIKDJO(I)  *  i    diagonal 
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10*  - 

urco  Program  -    71   Icetran  instruct  ions 
Object  Program  -    53   Binary  cards 

11-  CALXS  -  Son* 

12-  imJCaTIOMS  -  SUBSIST**  fUSM 

13-  3&FTC08  -  The  subroutine  examines  whether  or  not  a 
joint  is  eligible  for  being  a  "»ims  joint  of  a  member.  If  it 
is  the  number  of  the  corresponding  plus-joint  is  found  and 
stored  together  with  information  indicating  the  kind  of  the 
member. 

Although  the  flow  diagram  for  this  subroutine  is  the  sacie 
(ma^rosco.ically)  with  the  one  for  ME1HH0,  the  detailed  pro- 
cedure is  different.  That  was  the  r<sascn  for  writing  two  dif- 
ferent subroutines. 
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1  -  hAME  -  iiubrcutine  NSXMBQ 

2  -  AUTttOh  -  John  .^olllnlatis 

3  -  DATE  -  tfay  1,  1966 

4  •  im  -  Specific 

5  -  flion  MGtXFTXV  •  This  subroutine  is  used  in  the  case  of  a 

rectangular  plate  with  a  hole.  Its  task  is  to  assign  en  identi- 
fication number  on  each  somber*  Find  its  plus  and  minus  joint 
and  store  a  characteristic  number  (side  beam  or  diagonal)* 
This  information  is  essential  for  generating  the  necessary 
information  for  STHUDL  (M&fBER  ItfC^  ^)  as  well  as  for 

supplying  data  to  subroutine  KEtfPRG  for  the  purpose  of  compu- 
ting the  properties  of  the  members* 

6  •  UMM8I  -  Icetran 

?  •  FLOU  ChART  -  As  shown  in  Fig*  <A5,7)  •  Appendix  A. 
8  -  ~c;AG£;     -        CALL  ME124HO,  with  the  necessary  input  values  in 
COMMOU,  Input  variables  ares 
&W0J1 

■ten 
matt 

MM 

MM 

MQ 

MM 
,VEHB 


having  the  same 
meaning  as  in 

SDBftGuTlNS  DIPOHO. 
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NUKIDI 

mm 
m 


HOCOR  having  the  same 

VECOH     r  meaning  as  in 

mjQ  J  SUBROUTINE  MUGRID. 

9  *  QOTHJt  •  Output  in  COMHON  includes  four  dynasdc  arrays  of  siso, 
NUKE,  namely* 

NUHEMB(I)  *  Number  of  member  I  •  I . 
HUMWO(I)  «  Number  of  minus  joint  of  member  I . 
NOPUO(I)  ■  Number  of  plus  joint  of  member  I . 
KJJU}MS(I)  ■  Kind  of  member  I  . 

ONDHE(I)  *  0  if  member  I  is  a  sidebeam 
KINDHg(I)  *  1  if  member  I  is  a  diagonal 
10.  LENGTH 

Source  Program  -      139     ICSfRAN  instructions 
Object  Program        85      Machine  instructions 

11-  CALLS  .  None 

12-  APPLICATIONS  -  SUBSISTS*  PLATSS 

13-  LONG  DESCRIPTION 

The  task  of  this  subroutine  has  been  explained  in  the  short 
description  above*  For  the  purposes  of  explaining  how  this  program 
functions,  the  members  have  been  characterised  as  belonging  to 
one  of  four  reams  depending  upon  the  direction  of  their  minus- joint 
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to  plus- joint  vector*  These  four  teams  are: 
(i)   vertical  members. 

(11)  horizontal  members* 

(ill)  right  diagonal  members. 

(iv)  left  diagonal  members. 

The  numbering  of  the  members  proceeds  from  the  bottom  to 
the  top  and  from  left  to  right  in  each  one  of  the  above  teams. 
The  sequence  of  the  different  teams  is  as  listed  above*  In  other 
words  if  there  are  20  horizontal  members  and  30  vertical  ones*  the 
bottom  right-most*  right  oriented  diagonal  member  will  have  the 
number  51  and  so  forth* 

The  method  that  this  subroutine  follows  is  to  consider 
all  joints  that  are  legitimate  minus  joints  for  a  member  in  se- 
quence (joints  on  the  left  edge  of  a  hole  are  not  legitimate  minus 
joints  for  tight-diagonal  members  except  the  upper  most  one  for 
enable)  and  compute  the  corresponding  plus- joint*  The  number  of 
the  joint  as  well  as  the  minus  and  plus  joint  numbers  and  the  kind 
of  the  member  information  are  stored  in  the  appropriate  array 
elements  and  the  program  proceeds  with  the  next  member. 
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1  -  IMS  -  Subroutine  iObOC 

2  -  AUTHOR  -  John  Kolliniatis 

3  -  DATS  -  May  2t  1966 

4  •  TYPE  -  Specific 

5  -  SHORT  DESCRIPTION  -  This  subroutine  is  used  for  the  purpose 

of  matching  the  actual  supports,  with  the  nearest  node  of 
the  grid. 

6  -  LANGUAGE  -  Icetran 

7  -  FLOW  GttH  -  As  shown  in  Fig.  (A5.8)  Appendix  A 

8  -  USAGE  -  CALL  iQBGC,  with  the  necessary  input  values  in  COKMOK. 

Input  variables  are: 

10008   * 

VECGR  same  meaning  as  in  previous 

HUJO  SUBROUTINES . 

EX  ■  Distance  along  the  x-axis  between  the  center  of  coordi- 
nates and  the  particular  point  support. 

WHY  *  Distance  along  the  Y-axis  between  the  center  of  coordi- 
nates and  the  particular  point  support. 

KIND1*  Kind  of  point  support! 

1  -  if  support  is  hinged 

2  -  if  supi-ort  is  fixed 

9  -  GUTOT  -  Output  in  COMMON  includes 

KINUO:  An  array  of  sise  &J  containing  numbers  with  the  same 
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meaning  as  above  for  each  node,   except  that  zero  indicates  a 
regular  node  (no  releases).  This  information  is  necessary  for 
subroutine  JOBiRE  in  order  to  produce  the  proper  m4(XOH  COORDI- 
NATES" And  "JC.IHT  — ffftfflf  statements. 

10-  Hi 

Source  program  -  27  E39UH  instructions 

Object  program  -  30  Machine  instructions 

11-  CALLS  -  rione 

12  -APPLICATION  -  SUBSISTS*  PLATES 
13-  LONG  DESCRIPTION. 

For  a  given  point  support  the  distance  from  oach  joint  is 
computed  and  the  nearest  joint  is  located.  If  this  joint  is  not 
already  used  as  a  support,  it  is  considered  as  such. 

If  the  joint  already  i3  a  simple  support  and  the  point  sup- 
port is  fixed,  the  joint  is  reclassified  as  a  fixed  one* 

After  the  appropriate  joint  has  been  located,  information 
is  stored  in  .-LlrtDJO  regarding  its  kind. 
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1  -  mm  -  Subroutine  LIBOC 

2  *  A8ThQ8  -  John  Kolliniatis 

3  -  Ditt  -  Hay  2,  1966 
k  -  TH  I  -  Specific 

5  -  8SD83  b^OBU'TlOU  -  This  subroutine  takes  a  line  boundary 

condition  and  identifies  it  with  a  set  of  grid  nodes.  Onlj 
lines  along  the  4  edges  of  the  plate,  are  considered*  Any 
other  line  boundary  condition  can  be  approximated  by  a  sot  of 
point  supports,  but  the  User  has  to  do  it* 

6  -  LA&UtfAG£  -  Icetran 

7  -  FLOW  ChAivI  -  As  shown  in  Fig,   (A3* 9)  Appendix  A* 

3  -  USAGE  -  CALL  LI30C    with  the  necessary  input  in  CSmOH  (or 
S#IVAL£;CI2)  to  other  variables  la  €  ).  Input  variables 

are) 
KXhOJO 
hOCOR 
VECGR 
3BWQ 
0 
KX1    *    Starting  (left  side)  x-eoordinate  of  line  boundary 

condition . 
WH3EL  «  Starting  y-coordinate  of  line  boundary  condition. 
EX2    »  ijiding  (right  side)  x~coordinate  of  line  boundary 

condition  (l&UIVALSKCED  to  EX). 
WHX2  ■  Ending  y-coordirate  of  line  boundary  condition . 


same  meaning  as  in  iGBQC  . 
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(aSOIVALEMCED  to  WaT)  . 

KINDS  ■  Kind  of  lino  boundary  condition  (sWRfttOKM  to  HUM) 
KINB2  has  the  same  meaning  as  siimi  in  POBOC. 
9  •  OUTPUT  -  Output  in  COMMON  includes  KIUOJO  with  sane  meaning 

as  in  POBOC  . 
10.  LENGTH 

Source  Program  -  61  Icetran  instructions 

Object  Program  -  32  Machine  instructions 

11-  CALLS  -  None 

12-  APPLICATIONS  -  SUfiSISTEK  FLATUS 

13-  LUNG    DESCRIPTION  -  LIBOC  is  to  be  executed  after  POBOC  if 
there  are  any  line  boundary  conditions.  It  examines  what  joints 
satisfy  the  equation  of  the  line  boundary  condition  and  stores 
the  value  of  KINDS  into  the  corresponding  elements  of  KIf«BJ0(X). 
Of  course  before  it  does  that  it  examines  whether  the  line 
boundary  condition  coincides  in  all  or  in  part  with  one  of 

the  fc  edges  of  the  plate*  If  at  any  time  information  is  to  be 
stored  into  a  particular  element  of  ISssVO  which  already  con- 
tains a  number  other  than  sero  the  larger  val»e  if  finally 
kept  into  the  element*  In  other  wortis  If  we  have  two  line 
boundary  conditions  on  the  same  or  different  edges  tou ehing 
each  other  at  only  one  point  and  one  of  them  is  hinged  and 
the  other  is  fixed,  the  common  joint  will  be  considered  as 
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5.10,  ij'jjz^  * ..     -j:-a,o 

1  -  MM  -  Subroutine  rClNLC 

2  -  mOm  -  John  Kolliniatis 

3  -  Bftfl  -  May  3,  1966 
k  -  TOE  -  Specific 

5  -  ftMtt  UUNUflXQI  -  Initially  this  subroutine  was  written  in 

order  to  take  care  of  the  point  loads,  Lator  the  subroutine 
has  been  augmented  in  order  to  take  also  oare  of  the  displace- 
ments at  specified  points. 

6  -  LANGUAGE  -  Icetran 

7  -  FLOW  GUIS  -  As  shown  in  Fig.   (A5.10)  Appendix  A. 

8  -  USAGE  -  CALL  POBGC,  with  the  necessary  input  variables  in 

C0KK0W»  input  variables  are* 
80  JO 

ith  their  usual 
meaning. 


> 


m 

vscoa 

Ml 

P0L0H1  »  Point  load  magnitude. 

^QAD  *  x-coord    of  point  of  application  (i&l-IVALEJJCED  TO  EX) 
WILQAD  «  X-coord  of  point  of  application  (KJUIVALKSCED  TO  WbT) 

LO  *  iJustber  of  point  loads 
NLOTX  •  Type  of  loading 
1  ■  Z  FOHCB 
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2  *  x  ■■xm^m 

3  ■  Y  NOKBf 

.   =  A  variable  that  is  set  by  OR  and  explained  there. 
9  -  vUTHJT  -  Output  in  a  punched  form  consists  of  a  set  of  "JwINT 
ADS"  "PHIKT  DISPLACEMENTS"  and  "PRIOT  LOADS"  statements. 

10-  MM 

Source  Program  -  60  lOTWil  instructions 

Object  iTogram  -  46  Machine  instructions 

11-  CALLS  -  KQfOH 

12-  JtlfcTQftf  -  oU3SYST3H  PLA?v 

13-  14  iPTIOH  -  The  nearest  non-support  joint  is  found 
in  the  sane  way  as  in  POBOC.  Ilion  according  to  the  value  of 

-OUT  the  appropriate  statements  are  punched  out,  Tho  subrou- 
tine itself  can  handle  more  than  one  point  loads  but  in 
SUBSISTS*  plates  it  is  used  for  one  load      each  time. 


5.11.  .  .  -JJIWE  :jOS0JO 

1  -  NAME  -  Subroutine  KOSUJO 

2  -  AUTaOR  -  John  Kolliniatis 

3  -  DATE  -  May  2,  1966 

4  -  TXPE  -  Specific 

5  -  SHORT  DESCRIPTION  .  This  subroutine  is  called  after  LIBOC  and 

POBOC  have  been  called  in  order  to  generate  a  set  of  non-support 
joints 

6  -  LAIOAGE  -  Icetrar, 

7  -  FLOW  CHART  -  As  shown  in  Fig.   (A5.ll)  Appendix  A. 

8  -  USAGE  -  CALL  KOSUJO  with  the  necessary  input  variables  in 

COMMON  Input  variables  are* 

MOJO 

HOCOR  with  the  sane  meaning  as  in 

VECOR  every  other  subroutine  . 

KINDJO  , 

m 

9  -  OUTPUT  -  Output  in  COMMON  includes  : 


> 


NUJ01 

H0C0R1 

VEC0R1 
IZE 
10»  LENGTH 

Source  Program  • 

Object  Program  - 
11.  CALLS  -  Hone 


with  the  same  meaning  as 
MDJO,  HOCOR,  VECCR,   NJ 
respectively  but  for 
non-support  joints  only. 


21 

25 


Icetran  instructions 
Machine  instructions 
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12-  APPLICATIONS  -  Subroutine  POXAO  OF  SUBSISTS  PLATES. 

13-  The  value  of  KIMDJO  is  examined  for  each  joint.  If  it  is  aero 
the  joint  is  stored  as  a  non-support  joint. 


54. 


5-12.  -  -  -a   y  — 

1  -  NAME  -  Subroutine  JOINRE 

2  -  AUTHOR  -  John  Kolliniatis 

3  -  BATE  -  May  2,  1966 

4  -  TCPE  -  Specific 

5  -  SiiGRT  DESCRIPTION  -  Subroutine  JOINRS  is  used  in  order  to 

write  all  the  introductory  STKUDL  statements,  the  JOINT 
COORDINATES  table,  the  JOINT  RELEASES  statement  and  the 
MEMBER  INCIDENCES  table. 

6  -  LANGUAGE  .  Icetran 

7  •  FLOW  CnART  -  As  shown  in  Fig  (A5.12)  Appendix  A. 

8  -  USAGE  -  CALL  JCItfRE,  with  the  necessary  input  variables  in 

COMMON,  Input  variables  arei 

NUJO 

tiOCOR  with  the  sane  meaning 

VECOR        r  as  in  all  previous 

NJ  SUBROUTINES . 

KIKDJO      J 

JURENC  «  Joint  Release  number 

9  ■  CUTHJT  -  Output  in  COMMON  includes  all  the  statements  mentio- 

ned in  5  above,  in  a  punched  fona . 
10*  LENGTH 

Source  Program  -       89      ICETRAN  instructions 

Object  Program  -      55     Binary  cards 
11.  CALLS  -  None 


55. 


12-  APPLICATIONS  -  .STEM  FLAT 

13-  '•-  The  introductory  statement*  Are  punched 
first.  Then  according  to  the  value  of  EZWQ  a  table  of  JOIST 
COOHDIiiATi.S  is  punched  with  the  proper  EXIST.  ffifl)    sup- 
port#if  joint  is  a  support.  Subsequently  according  to  the 
value  of  JOREIAi(J)  set  by  LIBOC  the  proper  JOINT  R&fiASES 
statements  are  punched. 

In  particular  if 

JORSNU(J)  =  0  {  Joint  J,  releases  Moments  XVI 

JQRENU(J)  ■  1  ;  Joint  J  releases  Moment  I 

JORENO(J)  *  2  ;  Joint  J  releases  Moment  X 

JOREMJ(J)  ■  3  ;  Joint  J  releases  Moment  X,X 
Finally  the  "MEMBER  INCIDENCES"  table  is  punched  out. 
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5.13.  -   •■  ■-,:•  . .,, 

1  -  HAli£  m  Subroutine  H£KPRO 

2  -  MWB  -  John  uolliniatis 

3  -  DAT£  -  ^5ay  2,  1966 
^  -  TIPS  -  Specific 

5  -  WW  DESCRIPTION  •  Ibis  subroutine  computes  and  punches  STRUI& 

statements  for  the  properties  of  the  members  of  the  grid. 

6  •  la*—aT  -  Icetran 

?  -  FLOW  CHART  -  As  shown  in  Fig.  (A5.13)  Appendix  A. 

8  -  USAGE  -  CALL  MEMPRO,  with  the  necessary  input  variables  in 

COMMON.  Input  variables  are; 

iBHCKB 

WWi 

"WB 

wa 
wwt 

laipwo 


Have  same  meaning  as  in 
any  other  subroutine. 


> 


9  - 


;0R 
VECOR 

PRATI0  «  Poisson*s  ratio 
T..ICK     *  Plating  thickness 

OQTPUT  -  The  output  conists  of  a  set  of  acceptable  STRBDL     sta. 
tements  that  describe  the  properties  of  the  various  members. 
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10-  LEHGTH 

Source  nrogram  -  S3  ICKTRAN  instructions 

Object  rrot;ran  -  92  Binary  cards 

11-  CALLS  -  hone 

12-  Arf LIGATIONS  -  SUB3XSTEX  PLATSS 

1>  LCflG  DESCHimOH  -  Ike  kind  of  each  masher  is  exardnod  (3ide 
beam  or  diagonal)  and  then  after  computing  the  length  of  it 
and  all  the  necessary  dimensionless  ratio's,  the  properties 
of  the  member  are  computed  by  using  equations  2.18  through 
2.23  of  chapter  II,  Finally  all  the  above  information  is 
punched  out,  forming  part  of  the  STRUBL  input. 
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5.14.      --'■    ..-c-r.    -g--  ~- 

1  -  RAMI  -  Subroutine  LIHSLO 

2  -  MM  m  John  iolliniati* 

3  -  DATE  -  ?lay  3 1   i960 

4  -  fm  -  Specific 

5  -  Ml  HMBflM  -  Subroutine  LINE  is  used  in  order  to  find 

tne  appropriate  point  (nodes  or  points  along  the  side  meobere), 
where  concentrated  loads  have  to  be  applied  in  order  to  appro- 
ximate a  pacified  line  load.   It  also  handles  cases  of  output 
along  a  straight  line. 

6  •  lAMMtt  -  Icetran 

7  -  FLOW  ChART  -  As  shown  in  Fig,   (A5.14)  Appendix  A, 

6  -  USAGE  CALL  -  LZMLO  with  the  necessary  input  variables  in 
CO^iON.  Input  variables  are: 

HOCCR,  VKCOR,  &INDJO,  SDMEMB,  itfJMUO,  MUM*  NOVEHB,  ftUHOMB, 
KJ  with  the  sasre  meaning  as  in  previous  subroutines.  In  addition: 
U  -  SUrt  (left)  x-coordinate  (EQuIVALffiiCED  to  EX) 
I A  ■      **  "        y-coordinate  (RSOIVALENCeD  to  WliX) 

IB  -  end  (right)     x-coordinate  (&&IVALEMCSD  to  EX2) 
XL  *     ■         *  y-coordinate  (iS^IVALENCED  to  WStt) 

STALD  ■  Start  loading  Magnitude  (&1UIVALKNCED  to  F0L0K1) 
SNDL0  1  End  loading  Magnitude 
RLCTX  »  Same  meaning  as  in  rulNLO 
LTYPfi  •  0  if  uniform  load 
1  if  linear  load 
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9  -  OTPOT  *  Output  in  a  punched  fona  Includes  STRUDL  statements 
that  will  be  explained  in  the  chapter  for  LINE  since  this 
subroutine  performs  all  I/c  operations  for  LlftELQ 

10-  LLiGTH 

jrce  Program  -  230  Icetran  instructions 

Object  Pro£:rasi  -  2&I  Binary  cards 

11-  CALLS  -  LIRS 

12-  Ai'rT  ICATIONS  -  SUL.  •    PLATES 

13-  LONG  DESCRIPTION 

If  the  line  load  does  not  coincide  with  a  set  of  members,  then 
the  intersections  of  the  line  load  projection  with  all  members 
intersected  by  it  are  found.  The  number  of  the  intersected 
member  and  the  local  and  global  coordinates  of  it  are  stored 
for  later  use  by  LINE,  Special  consideration  is  also  given  to 
the  case  where  the  line  passes  thru  a  joint  or  a  small  area 
arround  it  in  order  to    be  considered  as  a  JOIiiT  LOAD  instead 
of  a  LINE    LOAD. 

In  case  that  tie  line  coincides  with  a  set  of  members  then 
the  nu-ber  of  the  coinciding  (with  the  line)  members  as  well 
as  local  and  global  coordinates  of  the  start  and  end  loading 
point  for  each  member,  are  stored  for  later  use  by  LINE, 
In  case  tnat  this  subroutine  is  used  for  a  line  output  certain 
parts  of  it  are  by-passed. 

Finally  subroutine  LINE  is  called  which  in  addition  to  other 
tasks  performs  also  all  punching  necessary  for  LINELO, 
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5.15.  ^  -  -^w  vm 

1  -  -  £ubx*o-itine  LI 

Z  ~  3ft  •  John  lolliniatis 

3  -  SATE  -  toy  3.  1966 
4-  -  tm  -  Specific 

5  -  SilORT  NbMftlPfltt  -  Subroutine  LI3E  is  a  logical  extension  of 

IXB3BA  sepai\*tedfrom  it  in  order  to  ;;et  arround  the  difficulty 
of  t^e  precompiler  to  accept  a  larje  number  of  symbols.  It  is 
called  by  LINELC  in  order  to  compute  the  different  grid  loads 
that  are  equivalent  to  the  plate  line  loads  specified.  In  addi- 
tion it  performs  all  the  I/O  functions  for  LIMELO  for  both 
the  case  of  line  loads  a  4  line  outputs. 

6  m  LANGUAGES  -  Icetran 

7  -  WM  (MM  -  As  shown,    in  Fig.   (A5.15)  Appendix  A. 

8  -  USAGE  -  GALL  LINK  (KiriDLO,  XGLLD,  YGLLD,  13UL0KE,  XLOCAL, 

In  addition  all  tne  input  in  COMMON  in  LBJELO  must  also  be  in 

lO.'J  in  LINE  plus  the  variable  KIOUT. 
The  values  transmitted  between  LINELO  and  LINE  as  arguments 
have  different  meanings  according  to  the  value  of  the  variable 
ffTWMff.  KWO  itself  is  an  array  of  aiae  :ftiLDS  which  contains 
numbers  0,1,2,  selected  as  it  was  indicated  in  the  flow  dia- 
gram for  fiTHteffii   (See  Appendix  E  for  meaning  of  NULOS) 
In  particular: 
(i)     if  KINDLO(J)  ■  0  (Member  concentrated  load) 
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XGLLD  -  ]  x  coordinate  of  concentrated  mepr-er  load 

YGLLD  *    .lobal  y  coordinate  of  ■  •  .« 

rQME  -  ftmber  of  loaded  member 
XIOC&L  -  Local  x  coordinate  of  concentrated  member  load 
tXOOtt  m  Local  y        *  •  •  •  « 

(ii)   if  amoto  «  i  (joirt  load) 

XGLLD  | 

r  *  Global  coordinates  of  loaded  joint 

HUB  J 

LOT  ■  dumber  of  loaded  joint 

(iii)  If  JUNOLO(J)  »  2  (Member  uniform  or  linear  load) 

XGLID  ■  Global  start  of  raerber  line  load 
YG1XD  «      *         end      "        ■  •  .     « 

XLOCAL  *  Local  start  of  member  line  load 
YLCCAL  ■  Local  end      *        "  ■.       • 

SULOME  •  Kumber  of  loaded  member 

In  this  case  it  is  necessary  to  say  that  XGLLB,  XGLLD, 
XiOCAL,  Y LOCAL,  all  refer  to  the  changing  coordinate, 
since  in  this  case  either  the  x  or  y  coordinate  remains 
unchanged. 
9  -  OUTPUT  -  Output  in  a  punched  form  depends  upon  the  values  of 
different  key  parameters.  Table  (515«1)  gives  the  different 
combinations  of  these  parameters 
10- 

Source  Program  -  269  Icetran  instructions 

Object  Program  •  25O  Dinary  cards. 
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11—  CALLo  «■    -cne 

13-  -    r"     T.':CL0(J)  ■  2  the  list  of  the  point*  that 

will  be  used  to  approximate  the  loading  is  shorted  according 
to  the  distant*  from  the  left  end  (or  lower  end  if  load  line 

is  parti  11  el  to  y  a:g±s)# 

If  (J)  *  Z  this  step  is  unnecessary  and  it  is  skippod. 

Each  one  cf  the  above  points  takes  so   tuch  of  the  :oad  as  it 
is  applied  between  the  middle  points  of  the  segments  joining 
it  with  it*  neighbouring  points.  Of  course  this  is  not  the 
case  when  OKDLC(J)  ■  2,  for  ail  J's,  Ir.  this  case  each  member 
is  loaded  with  the  rocact  linear  or  uniform  load  acting  over 
it.  Subsequently  the  necessary  output  is  punched  according 
to  the  values  of  the  control  variables  as  in  table  (5»15»1). 
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1  -  «AMii  -  SotootttllM  A&£4}#0 

2  -  iwi   u„  -  John  &oll litis  wis 

3  -  GATS  -  way  3,  i>66 
^  •  tXffi  •  opooifie 

5  -  Ui.ui*A  daSCftli'TJUMI  -  ouhroutine  AnrALw  is  used  in  order  to  treat 

area  loads.   u*e  method  is  used  to  divide  true  loading  among 
the  meaioers  covered  by  trie  area  load  and  then  call  SBBRQOfZtt 
Liit&Lu.  Only  ioads  with  perimeter  rectangular,  with  the  edges 
parallel  to  trie  edges  of  the  plate  can  be  treated. 

6  -  I*Hm*u.Huh  -  Icetran 

7  -  fIXM  UiAivi  -  As  anown  in  Fig,   (A5.16)  Appendix  A. 

8  -  **Auiu  -  CAU,  AR£*2*0  with  the  necessary  input  variables  in 

COiftSGH.  Input  variables  are? 

IGOft*  fB0Gft,  MHO,  WPtJO,  IPTOfB,  ;.>,  WLOTi,  with 

their  usual  meaning. 
aa  ■  x-eoordinate  of  lower  left  corner  of  area  load 

(WmMi  to  m) 

1A  *  y-ccordirate  of  lower  left  corner  of  area  load 
IB  a  x-coordinate  of  upper  right  corner  of  area  load 

(mmumiB  to  n) 

XB  *  y-coordinate  of  upper  right  corner  of  area  load 
(ttUIVAUUD  to  MOM 
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LB&STA  «  Starting  loading  magnitude  (at  lower  or  left  edge) 
LTXPE    *  A  variable  that  indicates  the  kind  of  the  load; in 
particulars 

Uniform  load 

Load  linear  along  the  x-axis 
*  •  ■  y-axis 


Icetran  instructions 
Binary  cards 


LTXPE  »  0 
LTYPB  =  1 

ltxpe  *  z 

10-  LENGTH  - 
»>ource  Program  -  i^q 
Object  Program  -            xqo, 

11-  CALLS  -  Subroutine  LXNELO 

12-  ki PLICATIONS  -  SUBSXST&H  PLATES 

1>  LQm  DoisCHXPTiOH  -  The  members  covered  by  the  load  are  first 
located.  If  LTXPE  is  0  or  1  the  load  is  distributed  a:  ong 
horizontal  members*  If  LTXPE  ■  2  the  load  is  distributed  MM| 
vertical  nembers.  Then  subroutine  LliiELC  is  called  in  order 
to  punch  the  necessary  output  statements. 


66. 


5.17.     ■    .     -JriNE  OUT 

1  -  NAME  .  Subroutine  OUT 

2  -  AUTHOR  -  John  Kolliniatia 

3  -  DATS  .fliy3i  1966 

4  -  TXi-E  -  Specific 

5  -  SHORT  DESCRIPTION  -  This  subroutine  examines  the  value  of 

KIOUT ,  which  is  set  by  the  uservs  data  or  COMMAND  (if  CDL 
is  used)  and  according  to  its  value  the  necessary  steps  are 
performed  by  CUT  or  other  subroutines  in  order  to  produce  the 
specified  OUTPUT. 

6  -  LANGUAGE  -  Icetran 

7  -  FLOW  CHART  -  As  shown  in  Fig*   (A5.1?)  Appendix  A. 

8  -  USAGE  -  Call  rCBGC  with  the  necessary  input  variables  in 

COMHO&  Input  variables  ares 

KIOUT  indicates  desired  output.  In  particular  if s 

KIOUT  ■  2  -  LIST  REACTIONS 

KIOUT  *  3  -  LIST  POINT  LOADS  XA..     XA.. 

KIOJT  «  k  -  LIST  LINE  LOADS  XA...  XA...  XB...  IB 

KIOUT  «  5  -  LIST  LOADS 

KIOUT  «  6  -  LIST  POINT  DISPLACEMENTS  XA...  XA 

tttfl  »  ?  -  LIST  LINE  DISPLACEMENTS  XA...  XA...  XB...  XB 

IIUUT  •  8  -  LIST  DISPLACEMENTS 

OUTXA  »  x  and  y  coordinates  of  point 

OUTXA  «  where  output  is  requested 
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GOTXA    ) 

OUTXA 

UJTXB 


x  and  y,  start  and  end  coordinates 
>  of  line  along  which  output  is  desired 


anro 

9  -  OUTPUT  -  the  output  is  in  a  form  of  punched  STRDQL  statements 
that  take  care  of  any  desired  collective  or  selective  output 
tnat  trie  USSR  desires. 

10-  MM 

Source  rrogram  -     ki  Icetran  instructions 

Object  program  -    27      Binary  cards 

11-  CALLS  -  LiKSLO,  PGIi&O 

12-  APi  LIGATION  -  SUBSIST©*  PLAT2S 

13-  LONG  DESCRIPTION  -  Flow  diagram  and  short  description  indicate 
completely  how  this  subroutine  functions. 


ft. 


Several  problems  have  been  runned,  for  the  purpose  of  deb- 
buging  individual  subroutines  and  examining  the  interaction  of  the»* 
Ihe  results  have  been  excellent  regarding  the  generated  outfit. 
Only  the  input  and  output  for  a  very  simple  problem  was  Included 
in  Appendix  F»  It  covers  the  case  of  a  square  plate  with  two  clamped 
and  two  hinged  o&ges  under  linear,  lateral,  area  load*  It  has  been 
selected  since  is  one  of  the  solved  problems  in  (11)* 

there  is  nothing  that  can  be  said  on  the  overall  accuracy 
of  the  method  since  the  non-availability  of  3TKUJDL  prevented  the 
autftor  from  solving  the  problem  completely,  however  judging  from 
the  results  of  (1),  where  a  grid  has  been  generated  manually,  the 
larger  deviation  from  values  measured  by  experiment  are  of  the 
order  of  15$. 

After  the  displacements  have  been  found  the  moments  and 
shears  can  also  be  evaluated  (-.sore  correctly)  from  the  following 
relations  given  in  (2). 


^--^%^%)  *" 


6y. 


^-<%^t) 


^'^T?'    3»9J^  (6.3) 


where 

w  *  lateral  displacement 


D*  ^^ -7T-  (6.5) 


ilC'-K1) 


Before  we  state  some  recommendations  for  future  expansion 
of  this  work,  there  are  several  other  points  that  must  be  said 
regarding  its  completion.  First  the  output  from  this  program  has 
to  be  tested  in  STKUQL  when  it  becomes  available*  Second  the  CBL 
program  (a  first  draft  of  which  is  shown  In  Appendix  D)   has  to  be 
compiled,  tested  with  dummy  subroutines  and  the  subroutines  of  this 
program.  This  will  make  the  program  for  rectangular  plates  in  fle- 
xure an  integrated  one. 

In  the  very  restricted  category  of  rectangular  plates  the 


70. 


bending  problem  is  not  the  only  problem.  There  also  axial  loads, 
instability  or  combinations  of  thesu  A  further  expansion  of  this 

work,  therefore,  would  be  to  write  similar  programs  for  these  cases 
and  then  put  everything  together.  It  would  really  be  a  great  contri- 
bution to  engineering  to  be  aWLe  to  analyze  any  plate  with  any 
combination  of  lateral  and  axial  loads  in  a  matter  of  -inutes. 
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SUBROUTINE  OUT 
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A?^m  ? 


WTPAfa  m^  of  PRttffiW 


9^. 


V-JOB KOI  I  INTATTS 


/RESTORE  OLD  69 
/LIBRARY  ICF.SI  TR 
/DUMP 
*»Q  2  6 


DYNAMIC  ARRAY  HOCOR » VECOR »NU JO , JORENU , KI NDJO,NUPL JO, NUMI JO >NUMEMB 

_DYNAMIC  ARRAY  K INDMF ,NUJQ1 .HOfOR 1  . VFTOR1 

REAL  LDMSTA»LDMEND 
.^COMMON  HOLE  Hj  V  FLFISUTHI  C  K  »  M  A  NU  JO  ,  YOUNG  ,  PRAT  I  0  ,  HOLEN  1  ,  HOL  EN2  ,  VFLEN1 
1>VELEN2»EX,WHY,KIND1 » EX 2 » WHY  2 »PCLOMl »KI L I LO ,LTYPE »NUHO Jl ,NUVE Jl » 
2NJ,HOCOR( P) »VECQR(P) »NUJO(P) > JORENU(P) ,HOLEN3 , VELEN3 ,NUHOJ2 , NUHOJ3 


3»NUVEJ2»NUVEJ3>NUPLJO(P) »NUMIJO( P) »NUMEMB(P) »<INDME(P> »NUM.E»NUVEMb 
4,NUH0MB,NURIDI,NULFnT  , NUPOLO, NU JO)  (P) ,HOCORl(P)  ,VFCORl (P)  ,NSIZE, 
5NULDS,ENDLD,LDMSTA»LDMEND»KINDJ0(P) »KIOUT»NLOTY 
__EQU  I  VALENCE  .  (  EX  ,  EX  1 ,  EX LOAD, X  A  >OUTXA  )  ,  (  WHY  ,'*HY1  ,  WYLOAD,  Y  A>QUTY  A)  , 


KKIND1»KIND2)  »  (  EX  2  »XB  »OUTXB  )  »  (  WHY  2  »  YB  »OUTYB  )  »  (  STALD  »  POLOM1  ) 
JfORMAT  (5F10.3) 

3  FORMAT  ( 16) 

4  FORMAT  (  I4,4F10.3) 

5  FORMAT  (214) 

_6JF_ORMAT__LIAj!  4F  1 0.3)  


7  FORMAT  (314) 

8  FORMAT  ( I5,3F10.3»I5,3F10.3 ) 


9  FORMAT  ( I5.4F10.3  ) 
CALL  INIT(QQDUB) 


READ  (5»2)  HOLEN, VELEN, THICK, YOUNG, PRATIO 
WRITE(6 ,2)  HOLEN, VELEN, THICK, YOUNG, PR AT  10 


READ  (5.3)  MANUJO 
WRITE  (6  »3_L.  MANUJO 


READ  (5»4>  NHOL»HOLENl»VELENl»HOLEN2»VELEN2 
WRITE (6,4)  NHOL, HOLEN 1,VELEN1,H0LEN2,VELEN2 


GO  TO  ( 50,60) ,NHOL 
50  CONTINUE 


THERE  IS  NO  HOLE 
HOLENl=HOLEN 


VELEN1=VELEN 
CALL  DIPONH 


10  CONTINUE 
CALL  GRID 


11  CONTINUE 

CALL  MEINNH 


12  CONTINUE 
GO  TO  70 


60  CONTINUE 
IHERJE_L3_^_H0l  F 


CALL  DIPOHO 
13  CONT INUE 


CALL  MUGRID 
14  CONTINUE 


CALL  MEINHO 
70  CONTINUE 


READ  (5,5)  NUPSU,NULSU 
WRITE(6,5)  NUPSU,NULSU 


IF  (NUPSU)  85,85,80 
80  CONTINUE 


DO  85  I=1»NUPSU 
HEAD   (5,6 )  KIND1 , FX ,WHY 


WRITE  (6,6)  KIND1,EX,WHY 
NUPSUP_=1 


95. 


.CAU-.-PX5BQC 


85  CONTINUE 

IF  (NULSU)  15*15,90 


90  CONTINUE 

„DO  15  1  =  1, NULSU 


READ       (5,6)     KIND2»EX1 ,WHY1,EX2»WHY2 
talR.II  E  _  L6j  6 1  .XI  ND?,FX1,WHY1  , FX?,WHY7 


NULSUP=1 
CALL  LIBOC 
15  CONTINUE 

CALL  JOINRE 


16  CONTINUE 
CAI  I  J4EMPRO 


100 

CONTINUE 
READ  (5,7) 

NUPOL, 

NULIL' 

»NUARL 

WRITE(6,7)  NUPOL, NULIL. 
JLF  (NUPOL)  110,110,105 

•NUARL 

105 

CONTINUE 
DO  110  1=1. 

NUPOL 

KIOUT=l 
NUPOLO=l 

READ   (5,9) 
WRITE  (6,9) 

NLOTY 
NLOTY 

,EXLOAD»WYLOADi 
,EXLOAD»WYLOAD< 

.POLOMl 
.POLOMl 

110 

CALL  POINLO 
CONTINUE 

IF  (NULIL)  120,120,115 
115  CONTINUE. 


DO  120  1  =  1, NULIL 
KIOUT=l 


READ  (5,8  )  NL0TY,XA»YA»3TALD,LTYPE»XB»YB,ENDLD 
_W  S ITE(6,8)  NLOTY,XA,YA,STALD,LTYPE,XB,YB,ENDLD 


CALL  LINELO 
120  CONTINUE 


IF  (NUARL)  130,130,125 
125  CONTINUE  


DO  130  1=1, NUARL 
KIOUT  =  l 


READ  (5,8)  NLOTY, X A » Y A, LDMST A »LTYPE » XB » YB » LDMEND 
WRITE(6,8)  NLOTY,XA,YA,LDMSTA,LTYPE,XB,YB,LDMEND 


CALL  AREALO 
130  CONTINUE 


READ   (5,5)  NUOUT 
WRITE  (6,5)  NUOUT 


DO  140  1=1, NUOUT 

READ  (5,9)  K I  OUT , OU TXA,OUTYA , OUTXB ,OUT YB 


WRITE<6,9)  KlOUT,OUTXA»OUTYA»OUTXB,OUTYB 
CALL  OUT  


140  CONTINUE 

141  FORMAT  (T8,8H  FINISH  ) 


WRITE  (7,141) 
DESTROY  HOCOR 
DESTROY  VECOR 
DESTROY  NUJO 


DESTROY  NUMIJO 
DESTROY  NUPLJO 


DESTROY  NUMEMB 
DESTROY  KINDME 
DESTROY  KINDJO 
DESTROY  JORENU 


J 


• 


96. 


DESTROY  NUJQ1 


DESTROY  HOCOR1 
DESTROY  VECOR1 
END 


57. 

26 

SUBROUTINE  DIPONH 

DYNAMIC  ARRAY  HOCOR»VECQR,NUJO»JC      .  I      -  ,      -    IIJOtNUl  EMB 

DYNAMIC  ARRAY  KINI    •       •   :ORl,VECC 

COMMON  HOLEN.VELEN,THICK,MANUJO,YOUNG»PRATIO,l     ll»HOLEN2»Vf 
1»VELEN2»EX>WHY,KIND1 »EX2 » WHY2 »POLOM 1  ,KlLILO,LTYl   ,1   HOJl  ,NUVEJ1, 
2NJ  ,HOCOR  (  P  )  ,VECOR  (  P  )  ,i\UJO  (  P  )  , ..      ..:  (  Pi  ,HOLEN3 , VELEN3 ,NUHOJ2 , NUH0J3 
3»NUVEJ2»NUVEJ3»NUPLJO<P) ,NUMIJO(P) ,        ( P ) ,KINDME ( P ) ,      .... 

>NUHOMB»NURIDI ,NULEDJ , NUPOLO » NUJOl [Pi  ,HOCORl ( P ) , VECOR1 ( P ) tNSlZE» 
5NULDS»ENDLD»LDMSTA,L[      ,  k  1  iMDJO  (  P  )  >rs  IOUT  ,  NLu  T  Y 

IF  ((->,   N  -1./3.  J  2>2,1 

1  WRITE  (6,4) 
CALL  EXIT 

2  CONTINUE 

4  FORMAT  (41H  POISSON'S  RATIO  LARGE » ANALOGY  NOT  VALID  ) 
RATIO=VELENl/HC 

ANUUO=MANUUC 
AUVEJ1=SQRT)/ 
AUHOJl=AUVEJl/RATIO 
NUVEJ1=AUVEJ1 
NUHOJ1=AUHOJ1 
INCVER=(NUVEJ1  +  D      Jl 
INCH0R=NUVEJ1*(NUH0J1+1 ) 
IF  (MANUJO-INCVER)  10,5,10 

5  CONTIi 

:      MAX  NUMBER  OF  UOINTS  ACHE  I'VE 
NUVEJ1=NUVEJ1+1 
GO  TO  5C 
10  CONTINUE 
:       EXAMINE  POSSIBILITY  OF  INCREASING  HORIZ. JOINTS 
IF  (MANUJO-INCHOR)  25*20,30 
20  CONTINUE 
:      MAX  NUMBER  OF  JOINTS  ACHE  IV 
NUH0J1=NUH0J1+1 
GO  TO  l 
25  CONTINUE 
:      IMPOSSIBLE  TO  INCREASE  EITHER  THE  HORIZ  OR  VERT  JOINTS 
GO  TG  : 
30  CONTINUE 
I      EITHER  THE  HORIZ  OR  VERT  JOINTS  CAN  BE  I NCREASED .MAKE  LARGER  CHANGE 
IF  ( INCHOR-INCVER )40»35,35 
35  NUH0J1=NUH0J1+1 

GO  TO  50 
40  NUVEJ1=NUVEJ1+1 
50  CONTINUE 

NUHOMB=  (NUHGJl-1  )  *NUVEJ1 

NUVEMB= (NUVEJ1-1 )*NUH0J1 

NURIDI= (NUHOJl-1 ) *(NUVEJ1-1 > 

NULEDI=NURIDI 

NUME=NUHOMB+NUVEMB+NUR"IDI+NULE 

NJ=NUVEJ1*NUH0J1 

DEFINE  KINDME,NUMI  ,   LL»HIGH 

DEFINE  NU  I    •      ,FULL,HIGH 

DEFINE  NUMIJ0,l     ,FULL,HIGH 

DEFINE  NUPLJ-*      ,FULL»HIGH 
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DEFINE 

DEFINE 

DEFINE 

DEFINE 

DEFINE 

DEFINE 

DEFINE 

DEFINE 

MANUJO= 

RETURN 

END 


VECOR1 ,NJ,FULL»HIGH 
HOCOR1 ,NJf FULL, HIGH 
NUJOl  ,NJ, FULL, HIGH 
KINDJO,NJ,FULL»H 
JORENU,NJ, FULL, HIGH 
NUJO,NJ,FULL,Hl 
VECCR  ,NJ, FULL, HIGH 
HOCOR  ,NJ, FULL, HIGH 
C 
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026 


1HL  iv       UF 

uF  JuIi^TS  lb 


JOINTS 


VALID  ) 


SUBROUTINE  JIPOHO 

DYNAMIC  ARRAY  HOCOR»VE(  ■  •      •        •        »iMUPL  Ju  ,NUMI  JU  j 

DYNAMIC  ARRAY  K I NDME »NU JO  1 >HOCOR 1 » VECOR1 

10LEN»VELEN,THICK,MANUJ0,Y0UNG»      [0»HOLENl »HOLEN2 » VELEN1 
1 >VELEN2»EX»WHY»KIND1»EX2»WHY2»P0L0M1,  ILILO,LTYl   »        •   IVEJ1, 
2NJ,:      (P)  »VECOR(P  )  »iMUJO(P)  »  J(J         •        ,  VELEN3  »NUHOJ  .. 
3 »NUVEJ2»NUVEJ3»NUPLJO( P ) sNUMI JO(P ) ,N      H P ) »K INDME ( P ) »N 
4»NUH    -        •      DI»NUPOLO»NUJ01 (P> jHOCORl (P) > VlCORI (P ) ,N  I 
5NULDS»ENDLD»LDMSTA»Lt      > iU  iMuJO  (  P  )  »K  I  OUT  ,NLOtY 
DIMEN       J( 6 )  »N  (6  )  »AN (6 ) ,A( 6  ) 
REAL  D 
C      THIS.     DUTINR  IS  CALLED  bY  F06R]  .IT  biVES 
C      FOR  EACH  PART  uF  THE  PLATE  IF  THI     '  N   BER 
C      SPECIFIED. THERE  IS  A  HOLE 
IF  ( PRAT  I  0-1./ 3. )  2»2»1 

1  W RITE  (6,4) 
CALL  EXIT 

2  CONTINUE 

4  FORMAT  (41H  POISSON'S  RATIO  LARGE > ANALOGY  NOT 
HQLEN3=HOLEN-HOLENl-HOLEN2 
VELEN3=VELEN-VELEN1-VELEN2 
D( 1 ) =H0LEN1 
D(2 )=H0LEN2 
0(3 ) =H0LEM3 
D( 4) =VELEN1 
D(5 ) =VELEN2 
D(6 ) =VELEN3 
DO  5  K=l,6 
A(K)=D(K)/D(D 

5  CONTINUE 
ALPHA  =  A(4)+A(5)+A(6j+A(2)-MA(4>+MZ)-^(6)+a(3)*A(4)+A(3)*A(5)  + 

1A(3)*A(6) 
BETA=1.+A(3)+A(4)+A (6) 
ANUJO=MANUJO 
DO  8  K=li6 

AIMUJ  =  4*ALPHA*  (  (  A(  1  )  /A(M  )**2  ) 
IF  (ANUJO-ANUJ)  7,8»o 

7  ANUJO=ANUJ 
MANUJO=ANUJO 

8  CONTINUE 
DO  10  K=l»6 

AN ) 1*=BETA/ALPHA.SQRT )  >   ) BETA/ ALPHA***2* • ANUJO/ALPHA* 
AN(K)=AN( 1)*A(K) 
N( K) =AN ( K ) 

10  CONTINUE 

N1=N  (  1)+N(  2  )+i\l(  3  )-2 

N2  =  N (4)+N ( 5)+N(6 )-2 

NJ=N1*N2-N(  2)  *N  (  5  )+2*(  N  (  2  )  +N  (  5  )  )  -4 

NUH0J1  =  N (  1  ) 

NUH0J2=N(2 ) 

NUHOJ3=N( 3 ) 

NUVEJ1=N(4) 

NUVEJ2=N( 5  ) 

NUVLJ3  =  M  6) 


100. 


NUHO 
NUVE 
NUVE 

NURI 

NULE 

DEFI 
DEFI 
DEFI 
DEFI 
DEFI 
DEFI 
DEFI 
DEF  I 
DEFI 
DEFI 
DEFI 
DEFI 
MANU 
RETU 
END 


J0=NUH0J1+NUH0J2+NUH0J3~2 
JO=NUVEJl+NUVEj2+NUVEJ3-2 

=  (NUVEJ0-1  )  *NUHOJG-(  IMUVEJ2-1  >  -  U\iUriuJ2-<£  ) 
4B=  (  NUHO  JO- 1  )*  (NUVEJO)-(!MUVEJ2-2)       J  J  2-1) 
DI=  (NUVEJO-1  )*(NUHOJO-l)-(  NUHUJ2-  1  )       I  J  2-1  > 
DI=NUR I 
=  NU'     t  .        IURI  DI+NULEDI 


NE 
NE 
NE 
NE 
NE 
NE 
NE 
NE 
NE 
NE 
NE 
NE 


KIN 

NUM: 

NUMIJO 

NUPLJO 

VEC 

H0C0R1 

NUJ01 

KINDJO 

J^RENU 


VECOR 
HUCOR 


NUME  ,FULL,HlGri 
NUME ,FULL,H 

,  HULL, HIGH 
NUME ,FULL »HIGH 
NJ»FULL,HIGH 
NJ» FULL, HIGH 
NJ, FULL, HIGH 
NJ»FULL»H IGH 
NJ »FULL»H I GH 


NUJO»NJ» FULL, HIGH 


NJ, FULL, HIGH 
NJ, FULL, HIGH 


J0  =  o' 

RN 


id. 

**0?A , ^ 

SUBROUTINE  GRID 

DYNAMIC  A RRAY—  id 0  C  Q  R  »  V  E CO R  » NU JO, JORFNU » K J NpJO.NUPL JO, NUMT JO, NUMFMB 

DYNAMIC  ARRAY  K INDME , NUJOl , HOCOR1 , VECOR 1 

COMMON  H0LEN3VELEN,THICK»MANUJ0jY0UNG,PRATl0,H0LENl»HOL£N2,VELENr_ 


1»VELEN2»EX,WHY,KIND1 »EX2»WHY2  »POLOMl » KI LI LO , LTY PE » NUHO Jl ,NUVEJl ♦ 
_2HL».HQCOR  (  P  )  %VFrOR(P)  ,Nt)JO(P)  ,  JORFNIK  P)   .HOI  FNH.VFI  FM^  .Nl  IHOJ  ?  .  Nl  IHO  J  "3 


3>NUVEJ2»NUVEJ3*NUPLJO(P) >NUMI JO( P) »NUMEMB( P) »KINDME( P) »NUME»NUVEMB 
4,NUHOMB,NURIDI  ,NUI  FOT  ,  NUPOLO, NUJOl  <P>  ,HOCORl  (P)  .VFCOR1  (P)  ,NSI7F, 
5NULDS,ENDLD»LDMSTA»LDMEND»KINDJO(P) , KIOUT,NLOTY 
EQUIVALENCE  (  NUHCJ1  »NUHOJO)  >  (  NUVFJ]  ,NUVFJO) _^_ 


*THIS  SUBROUTINE  GENERATES  A  GRID, GIVEN  A  RECT.  PLATE  WITHOUT  A  HOLE 

NVESPA  =  N1IVEJQ-1 

NHOSPA=NUHOJO-l 

_JdOSPAN=NHQSP-A 


VESPAN=NVESPA 
_X5PAC=HOLEN/HOSPAN 


YSPAC=VELEN/VESPAN 


NVESPA=NVESPA+1 
NHQSPA=NHOSPA+l 
DO  30  IO=l,NHOSPA 
1=10-1      


DO  20  JO=l,NVESPA 
J=JO-l 


K=I*NUVEJO+J+l 
AI  =  I 


IF  (I-NUHOJO+1)  12»10,12 
10  HOCOR (  KJ  =HOLEN 

GO  TO  13 
1?  CONTINUE 


HOCOR(K)=AI*XSPAC 
13  CONTINUE 
AJ  =  J 
IF  (J-NUVEJO+1)  15»14»15 


14  VECOR(K)=VELEN 
GO  TO  16 


15  CONTINUE 
_V£CQR(K)=AJ*YSPAC 

16  CONTINUE 
NUJO(K)=K  


20  CONTINUE 
J  =  0 


30  CONTINUE 

DO  31  J=1»NJ 


JORENU( J)=0 
31  CONTINUE  


RETURN 
END 


'     •  * 


./   4i  r 


102. 


11 


11 


112 


I 

■ 

2>WHY2>I       -     -  ■  ;    •       •      il, 

... 

•   .-LEi^l  •        •  -  -        [P  )  ,NSIZE  9 

.  •  .  .    ,  N  L  U  T  Y 

E  C  T  .  P  L  [  1       J  L  E 
Jl 

J  3 

I  (J  1-1.) 

JAC2  =  .'2-1.  ) 

XSPAC3=r  H0J3-1. 

Jl 

. 
A'.  J  3 

>AC1=VELEN1/(AUVEJ1-1. ) 
YSPAC2=VELEN2/ ( 

yspac:  .  .  ) 

J  = 

f  1 

• 
1  =  1     -1 

1-1 
A  I  =  I 
J  = 
I  F  {  :       P  <\ )  5  »  5  ,  7 

cc 

TO  THE  LEFT 
X=  I       CI 

-2  ) 
GO  TO  15 
CC      JE 

AT  Tl     )LE  OR  TO  THE  RIGHT  OF  IT 
NHOSPA=NUHOJ]        2-1 
IF  ( i - 1 ;  i      l ]  ? ,  9  9 1 1 

A  T     T i  i  I TS 

J 1-1.  I  .  J 1  +  1. )*XSPAC2 

L*(MJVEJl+NUVtJ;  3-2 '+( 1-NUnUJ] 

GO    TO     15 
C0NT1 
TO     1 

I0J1+NUH*  -3 

IF     (  J- 

113 

X=  (  J l-i.  )*:  3 J 2- 1. "J  -  •  -  . 


103, 


1XSPA 

113  C. 

= (NUVEJ+1 )*( I- 

Y  CC 

/EJ+1 
I  7  JNC     .  •    '.X 
J  =  .      -1 

1-1 
l  =  J 
IF  (J-NVE5PA)  17,17*19 
17 

Y  =  A  J  *  Y  S  P  A  C  1 

GO  TO  25 
lr-         JE 

AT  i  i      LE      :  .    _LY  OR 
IF  (  I-NUH0J1 )  194,  | 
192  CONTINl 

THE  RIGHT  OF  LEFT  El 


193 


194 


21 


23 


25 


AT  THE  LOl 


OVER  IT 


OF 


251 
252 


27 


IT  ( I-NUHOJ1-NUHOJ2+2)  193,194,1 
CONTINUE 

VERTICALLY 
IF  ( J-NUVEJ1-NUVEJ3 J  1  35»30»30 
CONTIN 

I  .  i         ,        PER  hDGE  OF 

Y=( AUVEJ1-1. )*YSPACl+VELEN2+( AJ-AUVEJ1 ) *Yir„ij 
GO  TO  2 
CONT  : 

out;      jle  x-  :  ITS  l 

NVESPA=NUVEJ]     EJ2-1 
IF  (J-NVESPA+1)  21,21,23 
CC 

A  T  T  i , .  I S E , BU T  0 L  I 

Y= ( AU      -I. )*YSPAC1+ ( AJ-AUVEJ1+1. ) *YSPAC2 
GO  TO 

T 

Y=l      J 1-1 
1Y5PAC3 
CONTIN 
FILL       f S 
- 
PERER=( l./l    . :  /ELEN 

)  251  •  -  .  • 
Y=VLLtN 

K=N      ^+J+1 
KI  =  K 

NUJ0(K)=K1 
HOCOR (K)=X 
VECOR(K  )=Y 
CON!  ' 
CONTIf 

J=  1  »f 


:    ■ 


■1. )*YSPAC2+( AJ-AUVEJ2-AUVEJ1+2. )* 


104. 


U(  J 

31  CONTI 
41  CONTI 

EN 


105. 

26 

SUBROUTINE  MEINNH 

DYNAMIC  AH  » VECUR  »NUJu»  JuRENU  »i<  I      -  , 

IIC  ARRAY  KIND!   »NUJ01 ,HOCOR 1 > VEO 

COI      HQLEN»VELEN,THICK»ivlANUJ  •       »l  <AT 10 >HQLEN 1 >MOLEN2 » VELEN1 
1 »VELEN2 »EX»WHY •      1  » EX2 , WHY2 , P0L0M1 , K I L I LO , LTY   ■  UHUJ 1 ,NUVEJ 1 , 
2NJ,i-     (p)  tVECOR(P)  ,NUJO(l   ,  I,     |   ,\       ,NUH0J2 ,NUhOJ3 

3  »NUVEJ2,i.       »NUPLJO(P)  »NUM"IJU.(P)  i       .<P  )  ,KIN   It  (P  )  , 
A,NUHOMB»NURIDI  ,NULEDI  ,NUPOLO  »  N:U  J01  (  P  )  ,H      L  (  P  )  »  VEC0R1  (  P  )  ,NS  I  ZE  , 
5NUL0S»ENDLD»LDiviSTA.        ,1  [  NDJO  (  P  )  ,K  IOUT  »  NLOTY 

EQUIVALENCE  ( NUHOJi ,NUHOJO ) , (      Jl.NUVEJO) 

DO  4  I=1»NUME 

KINDME (  I  )=o 

4  CONTINUE 
NUJ=NUHOJO*NUVEJO~l 
K  = 

DO    2.     I=1,NUJ 

IF     (HOCORt I )-HOCOR( 1+1) )     10,5 » 10 

5  CONTINUE 
J=I-K 

|S|UMIJO(  J)  =NUJO(  I  ) 

NUPLJO(U) =NUJO( 1+1 ) 

NUMEMb( J ) =J 

GO  TO  20 
10  CONTINUE 

K  =  K+1 
20  CONTINUE 
I      HOR I Z. MEMBERS 

K=0 

NUJ=NUVEJO*(NUHOJO-l ) 

DO  40  1=1,NUJ 

J=I+ ( NJVEJO-1 ) *NUHOJO 

NUMI J0( U) =NUJO( I ) 

K=I+NUVEJO 

NUPLJu( U ) =NUJO< K ) 

NUMEMBI J)=J 
40  CONTINUE 
:       DIAGONAL  I      -:S 

NUJ=NUVEJO*(NUHOJO-l )-l 

K  = 

DO  6  0  I=1»NUJ 

RATIO=VECOR( I ) / VECOR ( NUVEJO J 

IF  (RATIO-1.)  45,50»45 
45  CONTINUE 
:      THERE  IS  A  DIAGONAL 

J=I+ (NUVEJO-1 )*NUHOJO+ (NUHOJO-1 I       JO-K 

NUMI JO ( J)=NUJO( I ) 

NUB= I+NUVEJO+1 

NUPLJO( J) =NUJO(NUB) 

NUMEMB( J ) =J 

KINDME( J) =1 

GO  TO  6 
5j    CONTINUE 

NO  DIAGONAH  JOINT  I:     PERIMETER) 

K  =  K+1 


60  CONTINUE 

NUJ=NUVEJO*NUHOJO-l 

NUJ1=NUVEJ0+1 

K  = 

DO  80  I=NUJ1»NUJ 

RAT IO=VECOK( I ) / VECOR (NUVE JO ) 

IF  (RATI  0-1.)  63,70»b'j 
65  CONTINUE 

THERE  IS  A  DIAGONAL 

J=I-K+3*NUVEJO*NUHOJO-3*NUVEJO- 

NUMI JO ( J ) =NUJO( I ) 

NUB= I-NUVEJO+1 

NUPLJOC J ) =NUJO(NUI 

NUMEMBt J ) =J 

KINDME( J) =1 

GO  TO  ;; 
70  CONTINUE 

NO  DIAGONAL 

K  =  K+1 
80    CONTINUE 

RETURN 

END 


10*, 


2*IMUHG  JQ+  i 


107. 

26 

SUBROUTINE  MEINHO 

DYNAMIC  ARRAY  HC     -.     * > NU JO » JORENU » K  I      »NUPL JO ,NUM I  JO , 

DYNAMIC  ARRAY  K INDME ,NUJ01 »HOCOR 1 » VI 

HOLEN,VELEN,THIC  •    I  JO,YOl    •   AT  10 »HOLEN 1 >H0LEN2 » VELt 
1  »VELEN2»EX»WHV  •      1  ,  EX2  ,WHY2  ,P0L0M1»KILI  LO  »  LTYPE  .IMUHOJ 1  »NUVE Jl » 
2N.  •       (P)  »VECOK(P  J  ,iMUJO(  P)  »JOKENU  (P  )  »Hul     » VELEl\3  »NUH0J2  ,     JJ3 
3»NUVEJ2»N   =  $,NUPLJO(P)  ,NUMUO(P)  ,1       (  P  }  »K  INDMt  (  P  )  ,N 

•  JLEDI»NUPOL0»IMUJ01  (P)  »HC     L  (  P  )  »  VEC0R1  {  P  )  ,  NS  I  ZE  , 
5NULDSjENDLD»LDMSTA,L|      •      JO(P) tKIOUTtNLOTY 
NUVEJO=NUVEJl+NUVEj2+NUVEJ3-2 
NUHOjO=NUHOJl+NUHOJ2+NUHOJ3-2 

3  FORMAT  (315) 
DO  4  I=1,NUME 
KIND ME ( I ) 

4  CONTINUE 

C      VERTICAL  MEMBERS 
NUJ=NJ-1 
K  =  U 

DO  20  I=1»NUJ 

NALPHA=NUVEJO*l  IVEJ1-D 

NBETA=NALPHA+(NUHOj2-3 ) * ( NUVE J 1+NUVEJ3 )+2 
IF  (  I-NALPHA)  12  ,  12  ♦  11 

11  CONTINUE 

IF  (I -N BETA)   13 » 12*12 

12  CONTINUE 

C      REGULAR  PROCEDL 

IF  (HOCOR( I )-HOCOR( 1+1 ) )  122,121,122 

121  CONTINUE 
J=I-K 

NUMI J0( J ) =NUJO( I ) 
NUPLJUl J ) =NUJO( 1+1 ) 
NUMEMB( J)=J 
GO  TO  20 

122  CONTINUE 
K  =  K+1 

GO  TO  20 

13  CONTINUE 

IF  ( VECOR (  I  )-VECOR(NUVEJl )  )  131»122,131 
131  CONTINUE 
GO  TO  12 

20  CONTINUE 

C      HORIZ. MEMBERS 

NALPHA=NALPHA-( NUVEJO-1 ) 

NUJ=NJ-NUVEJO 

K-b 

DO  30  I=1,NUJ 

IF  (  I-NALPHA)  22,22,21 

21  CONTINUE 

IF  (I-NBETA)  23,22,22 

22  CONTINUE 

C      REGULAR  PROCEDURE 
J=NUVEMB+I-K 
NUMI JO( J ) =NUJO( I ) 
NSUB=I+NUVEJO 


108. 


;NUJO( NSUB) 
J 


233,233,231 


NUPLJu( J) 
NUMEMB( J ) 
GO  TO  3. 
23  CONTINUE 

IF  ( VECOR (  I  )-VECOR(NUVEJl  )  ) 

231  CONTINUE 
NUB=NUVEJ1+NUVEJ2-1 

IF  (  VECOR ( I  J-VECOR(NUB)  )  232,233,233 

232  CONTINUE 
K  =  K+1 

GO  TO  3o 

233  CONTINUE 
J=NUVEMB+I-K 

NUP  = I+NUVEJO-( NUVEJ2-2 ) 
NUMI JO < J ) =NUJO( I ) 
NUPLJu( J ) =NUJO(NUP ) 
NUMEMB( J ) sj 
GO  TO  3'. 

30  CONTINUE 

RIGHT  DIAGONAL  MEMBERS 

NALPHA  =  NA!_PHA-1 

NUJ=NUJ-1 

K=0 

DO  40  1  =  1,  NUJ 

IF  ( I-NALPHA)  32,32»31 

31  CONTINUE 

IF  (  I-NbETA)  3JOZ  ,32 

32  CONTINUE 
REGULAR  PROCEDURE 
NUP= I+NUVEJO+1 

323  CONTINUE 

IF  ( VECOR (  I  )-VECOR (NUVEJO) 

321  CONTINUE 
J=NUVEMB+NUHOMB+I-K 
NUMI JO ( J) =NUJO( I ) 
NUPLJo( J ) =NUJO( NUP ) 
NUMEMB( J) =J 
KINDME ( J ) =1 
GO  TO  ^ 

322  CONTINUE 
K  =  K+1 
GO  TO  4. 

3  3  CONTINUE 

VICINITY  OF 
IF  ( VECOR ( I 

331  CONTINUE 
NUB=NUVEJl+NUVEJ2-2 
IF  ( VtCOR ( I )-VECOR (NUB ) ) 

332  CONTIimUE 
NUP= I +NU V E JO-NU VE J  2  +  3 
GO    TO    323 

40    CONTINUE 

LEFT  DIAGONAL  MEMBERS 

NALPHA  =  NUVEwL-     IOjl  +  (  NUVEJ1-1  ) 

NBETA=NBETA+NUVEJO 


321  ,322  ,321 


HOLE 
-VLCuR(NuVcJl  )  ) 


jjZ)  j>3  1  »  j:J  1 


5l<Li5lli'l2<L 


> 


*. 

109, 


NUJ=NJ-1 

NUJ1=NUVEJ0+1 

K  =  C 

DO    50     I=NUJ1*NUJ 

If      ( I-N  I     ^2,42,41 

41  C. 

IF(I-NBETA)     4 3, 42, 42 

42  CONTINUE 
REGULAR  PROCEDL 
NUP= I-NUVEUO+1 

423  CONTIN 

IF  ( VECOR( I )-VECOR(NUVEJO) )  421-»422>421 
421  CONTINUE 

J=NUVEMB+NUHi  EJO+I-K 

WRITE  (6,3)  I  ,  J  ,  K 

NUMI JO( U ) =NUUO( I ) 

NUPLJO< J ) =NUJO( NUP ) 

NUMEMB( J) =J 

KINDME( J)=l 

GO  TO  bJ 
42  2  CON T  I.N 

K  =  K+1 

GO  TO  '- 

43  CONTINUE 
VICINITY  OF  HOL 

IF  ( VECOR ( I )-VECOR(NUVEJlT)  432,431,4^1 

431  CONTINUE 
NUB=NUVEJl+NUVEJ2-2 

IF     (  VtCOR  (  I  )-VECUR(  NUB)  )     4  22  ,  ^l  2  » <+  32 

432  CONTINUE 
NUP=I-NUVEUO+NUVEJ2~l 
GO    TO    423 

50    CO-NTINUE 
RETURN 
END 


110, 


_**£L2-k- 


SUBROUTINE  POBOC 
DYNAMIC  ARRAY  HOC-OR»V£CLQR»NUJ 
DYNAMIC  ARRAY  K I  NOME » NUJ01 »H0 
COMMON  HOLENj^£LENjlIHICK,MANU 
1»VELEN2»EX»WHY»<IND1 *EX2»WHY2 
-2iJJ-aHQC0R(P)  WFCQR(P)  »NLUO(P)  * 


0 » J0R£NiUKINDJ0tNUPLJQ*NUMlJQ^NUMEMB 
C0R1.VEC0R1 

J0>Y0UNG»PRATlQ»-H0LENl»HQLEN2>VELFNI 
»P0L0M1»KILIL0,LTYPE»NUH0J1>NUVEJ1» 
JORENI.K  P  )  »HOLFN3»VELEN3»NI.)Hn.J2»NUHn.n 


3>NUVEJ2>NUVEJ3>NUPL 

4,NUH0MB,NURIDI  *NIII  F 

5NULDS»ENDLD>LDMSTA> 

_IHIS  SUBROUTINE  JLAK 


COORDINATES»EX( I )AN 
TT  WITH  THF  NFARF^T 


JO(P) »NUMI 
DJUNUPOLO^ 
LDMEND,KIN 
FS  EACH  PA 
D  WHY( I )  A 
GRID  JOIN 


JO(P) »NUMEMB(P) »KINDME(P) »NUME»NUVEMB 

NUJ01  (  PJ .3.RQCQR1  (P)  tVFCORIC P)  iNSIZ.F.  » 

DJO(P) »KIOUT»NLCTY 

R  TCUI  AR  POTNT  SUPPORT  W I T  H 

ND  OF  KINDtKINDK I )  AND  MATCHES 

T  THF  NIIMRFR  OF  WHICH  T  £ 


NUJO(ALPHA)  WITH  COORD,  HOCOR(ALPHA)  AND  VECOR ( ALPHA ) . 


D  I  ST  1  =  (  E* -HOCOR(  1  )  )**?  +  {  WHY -VFCOR(  1  )  )  **2 

NALPHA=1 

DO  100  J=2,NJ 


DIST=(EX    -HOCOR( J) )**2+(WHY 
tf  inT.^Ti-ni.sTnn.io.^o 


-VECOR( J) >**2 


10  CONTINUE 

DISTA!iC£_JERQM_ 
GO  TO  100 

30    CONTINUE 


NEW  POINT  GREATER  THAN  OR  EQUAI TO  PREVIOUS 


DISTANCE  FROM 
PI  ST  1=01  ST 


NEW  POINT  LESS  THAN  PREVIOUS 


NALPHA=J 

00  CONTINUF 

NEAREST  JOINT 
If 


FOUNT 
IJ<  I N  D  J  0(NALPHA)-KIND1 


)  150,700)200 


150  KINDJ0(NALPHA)=KIND1 
200  CONTINUE 


RETURN 
_ENQ 


111. 

**0?6 

SUBROUTINE  L I BOC 

DYNAMIC  ARRAY  HOCOR > VECOR ,NU JO, JORENU, KI NDJO, NUPL JO , NUM UO,NUMEMB 

DYNAMIC  ARRAY  K I NDME »NUJ01 » HOCOR 1 » VECOR 1 

COMMON  HOLEN,VELEN,THICK,MANUJO,YOUNG,PRATIO,HOLENl>HOLEN2,VELENl 
1»VELEN2»EX,WHY,KIND1 » EX  2 , WHY  2 » POLOM1 » K I L I LO , LTYPE »NUHO Jl * NUVE Jl » 
2NJ>HOCOR(P) »VECOR(P) >NUJO(P) > JORENU( P) ,HOLEN3 , VELEN3 >NUHOJ2 , NUHOJ3 
3»NUVEJ2»NUVEJ3»NUPLJO(P) »NUMI JO ( P ) »NUMEMB(P) ,KINDME(P> »NUME»NUVEMB 
4,NUH0MB,NURIDI  »NULEDI  , NUPOLO*  NU  JOl  (  P  )  »HOCORHP)  »VECORKP)  »NSIZE» 
5NULDS»ENDLD»LDMSTA»LDMEND»KINDJO( P) »KIOUT»NLOTY 

EQUIVALENCE  (EX.EX1  )  »  (  WHY,  WHY  1)  ,  (  K  IND1_»KI  ND2J 


C      THIS  SUBROUTINE  TAKES  A  SET  OF  LINE  SUPPORTS  AND  APPROXIMATES 

Jl THFM  BY  A  SFT  OF  FQUIVALENT  GRID  JOINT  SUPPORS 

DO  100  J=1,NJ 
_C FIND  EXACTLY  MATCHING  COOR. 


IF  (EX1    -EX2    )  20,10,20 
10  CONTINUE 


c 

X-COORD.  MATCHES 

IF  (EX1    -HOCOR(J))  11,12,11 

c 

11  CONTINUE 

NOT  A  JOINT  ON  BOUNDARY 

GO  TO  100 
12  CONTINUE 

C 

A  JOINT  ON  BOUNDARY 

IF  (WHY2    -VECOR( J)  )  13*14,14 

13  CONTINUE 

JOINT  OUTSIDE  Y-RANGE 


GO  TO  100 
14  CONTINUE 


JOINT  PROBABLY  INSIDE  Y  RANGE, CHECK  Yl 
IF  (WHY1    -VECOR(J))  16,16,15 


15  CONTINUE 

JOINT  OUTSIDE  Y-RANGE 


GO  TO  100 
16  CONTINUE 


JOINT  DEFINITELY  WITHIN  Y-RANGE 
JORENU( J)=JORENU( J)+l 


GO  TO  80 
20  CONTINUE 


Y-COORD.  MATCHES 

IF  (WHY1    -VECOR(J))  21,22,21 


21  CONTINUE 

NOT  A  JOINT  ON  BOUNDARY 


GO  TO  100 
22  CONTINUE 


A  JOINT  ON  BOUNDARY 

IF  (EX2    -HOCOR(J))  23»24>24 


23  CONTINUE 

JOINT  OUTSIDE  X-RANGE 


GO  TO  100 
24  CONTINUE 


JOINT  PROBABLY  INSIDE  X-RANGE , CHECK  XI 
IF  (EX1    -HOCOR(J))  26,26,25 


c 

25  CONTINUE 

JOINT  OUTSIDE  X-RANGE 

GO  TO  100 

26  CONTINUE 

C 

JORENUt J)=JORENU( J ) +2 

JOINT  DEFINITELY  WITHIN  X-RANGF 

- 


112. 


80    IF     (KINPJCK  J) -KIND?  )     90*100,100 


90  KINDJO( J)=KIND2 
100  CONTINUE 


PICK  ANOTHER  JOINT 
LPO^JIONTINUE 


PICK  ANOTHER  BOUNDARY  CONDITION 
RETURN 


END 


JUtSlZB- 


U3. 

SUBROUTINE  POINLO 

DYNAMIC  ARRAY  HOCOR»VFC  OEUALU  J  Q  *UQRE.NU+K  I-ND-J-Q-,  NUPL  JO,NU  M  I  J  O-^J^IUME-MB- 
DYNAMIC  ARRAY  K I NDME , NU JOl » HOCOR 1 , VECOR 1 
_CQMMQN__HOLEN.  ,  VELEN  ,  TH  I  CK  »MANU  JO,  YOUNG*  PRAT  I  Q,HULFN  1  ,  HOLFN2  ,  VFLFN1 


1»VELEN2,EX,WHY,KIND1 » EX 2 , WHY 2 tPOLOMl , K I L I LO ,LTYPE ,NUHO J  1 , NUVE Jl » 
.2ALU  HDCQR1R  J_*  V  ECOR ( P  )  ,N.UJO(P)  , JQRFNIK  P)  ,HOl.FN3,VFI  EN3  ,NUHOJ?  ,  NUHOJ3 

3,NUVEJ2»NUVEJ3»NUPLJO(P) ,NUMI JO( P) ,NUMEMB(P) ,KINDME( P) , NUME,NUVEMB 
_4_*NUH0MB  ,  MURUU  tNUI.  EDI  *  NUPOL  ,Qi  NU  JO]  (P)  ,HOCORl  (P)  ,  VECTOR  1  (  P  )  *J\1£1ZEj. 


5NULDS»ENDLD»LDMSTA»LDMEND,KINDJO(P) »KI0UT» NLOTY 
IU I  VALENCE  IEXj  EX LOAD*  t ) WHY,WYLOAD* 
THIS  SUBROUTINE  TAKES  EACH  POINT  LOAD  AND  FINDS  THE  NEAREST  NODE 

-C OF-JT-HE  GRID  WHTCH  TS  NOT  A  SUPPORT. TO  RF  rONrTDFRFH  AS  POINT  OF 

C      APPLICATION  OF  LOAD. 
. _7_n__F_QRMAJ  (T8,1ZH  JOINT  1  PADS) 


75  FORMAT  (T8,I4,9H  FORCE  Z  ,F15.3) 
JJLJEQRMAT  (T  8,I4,10H  MOMFNT  X  ,F15.3) 


83  FORMAT  (T8»11H  LIST  LOADS  ,14) 

85  FORMAT  (Tfl.U.lflH  MOMFNT  Y  .  F  1  5  .  1  ) 


86  FORMAT  (T8,20H  LIST  DISPLACEMENTS  ,14) 
F-ORJVL.AAL  ARRAY  OF  NON-SUPPORT  JOTNTS 


CALL  NOSUJO 

IF     (NUPQLO-1)  _22jl21a22. 


21  CONTINUE 
WRTTF  (7.70) 


2  2  DIST1=(EXL0AD-H0C0R1( 1 ) ) **2+( WYLOAD-VECOR 1 ( 1 ) )**2 
NALP H A  = 1 


DO  TOO  J=2,NSIZE 

DIST    =(EX LOAD-HOCORK J ) ) **2+( WYLOAD-VECOR 1 ( J) )**2 

IF    (DIST1-DIST)     100,100,30 

30    D I  ST  1=0  1ST 

NALPHA=J 
.JLOO_CDfcLT_IilLIE 


NEAREST  JOINT  FOUND 

GO  TO  (  1 0  1,101  ,103,103,103,106)  ,KIOUT 


103  CONTINUE 

NUJ0q  =  NUJ0l  (NAI  PHA) 


WRITF  (7,83)  NUJ09 
GO  TO  200 


106  CONTINUE 

NUJ09= NUJ 01 (NA LPHA) 


WRITE  (7,86)  NUJ09 
GO  TO  ?0Q 


101  CONTINUE 
LE_  (NLOTY^LL92,91,92 


91  CONTINUE 
N LO I Y  CO  N  T  A  J  N  S  A  N  l.LOAD  IS  A  Z  FORCE 


NUJ09=NUJ01 (NALPHA) 
WRTTF  (7.75)  NUJOQ.POIOMI 


GO  TO  200 

92  CONTINUE 


NLOTY  DOESN'T  CONTAIN  AN  1. CHECK  IF  2  DOES 
IF  (NLOTY -2)  94,9  3,94 


93  CONTINUE 

Nl OTY  CONTAINS  A 


NUJ09=NUJ01 (NALPHA) 

WRITE  (7,80)  NUJ09»RQLQM1. 

GO  TO  200 


.9.4_CQMTINUE 


■ 


il*. 


.NLQIX—CQUXAmS-^-3- 


NUJ09=NUJ01 (NALPHA) 

WRITE  (7,85)  NUJO-2-i.eaijQMl. 

200  CONTINUE 

. NUPOLO=NUPOLO+l 
RETURN 

END 


115. 


j**QZ6- 


IPLJ0»NUMIJ0»NUMEMB_ 


SUBROUTINE    NOSUJO 

DYNAMIC    ARRAY    HQCORLjVECOR  ,NU  JO,  JORENLL 

DYNAMIC    ARRAY    K I NDME, NUJOl »HOCOR 1 , VECOR1 

COMMON    H0LEN»VELEN»T_HLLCK_?MANUJ0»Y0UN5»PRATI0»H0LEN1»H0LEN2»VELEN1 
1»VELEN2»EX,WHY,KIND1»EX2»WHY2»P0L0M1,KILIL0,LTYPE»NUH0J1,NUVEJ1» 
■2NJj.HQ.CQRl£J_i.VFrQR(P)  iNUJQ(P)  »  JORFNtH  P)  .HOI   FN^VFI   FI\R  .Nl  1HOJ7  ,  Nl  iHOJ^t 
3»NUVEJ2»NUVEJ3»NUPLJO(P)  »NUMI JO( P) ,NUMEMB(P>  >KINDME(P) »NUME»NUVEMB 
4,NUH0MB,NURIDI  .NULEDI  >  NUPOLO*  NUJOl  (EJ  >HOCOR  l.(  P  )  »  VECQRIXP-LaNSLZE  » 
5NULDS,ENDLC»LDMSTA»LDMEND,KINDJO(P) »KIOUT»NLOTY 

K  =  0 


DO  50  1=1, NJ 

IF  UCLNDJOU_)-0) 


-lD.»5ilO 


CONTINUE 

NOT  A  SUPPORT 

J=I-K 

NUJOl  (  J)=NUJ.Qi_Ll_ 


HOCORK  J)=HOCOR(  I  ) 
V£ COR  1  (  J  )  =VECOR  (  T  ) 


JLQ. 


50 


GO  TO  50 
CONTINUE 
A  SUPPORT 
K  =  K+1 

CONTINUE 

NSIZF=J 


JOINT  SKIP. 


RETURN 
END 


116* 


**f>26 — 

**OBJECT 

SUBROUTINE  JOIMRE 


DYNAMIC  ARRAY  HOCOR , VECOR ,NUJO , JORENU » KI NDJO, NUPL JO, NUMI JO , NUMEMB 
DYNAMIC  ARRAY  K I  NOME »NU JOl , HOCOR 1 , VECORl 


DYNAMIC  ARRAY  NU 

C  QMMQ&LJdQi_EAUV£.l  EN  ,  TH  I  CK  ,MANlJJO  , YOUNG  ,  PRAT  I  0,  HOI  ENl»HOI  EN2,VEI  EN1 


1>VELEN2,EX,WHY,KIND1»EX2»WHY2»P0L0M1,KILIL0,LTYPE,NUH0J1,NUVEJ1, 
2NJ,HOCOR( P  )  ,VECOR  IP)  ,NU  JO  (  P  )  ,  JORFNU  (  P  l.HOl  FN3.VFI   FN3 »NUHOJ? . NUHOJ3 


3,NUVEJ2»NUVEJ3,NUPLJO(P) »NUMI JO( P) »NUMEMB(P> ,  KINDME(P) ,NUME,NUVEMB 
4 , NUHOMB ,  NUR I D I ,NULEDI ,NUPOLO, NU JQ1 (P) ,HOCORl (P) ,VFCORl ( P ) , NS I ZEi 
5NULDS,ENDLD»LDMSTA»LDMEND»KINDJO(P) »KIOUT»NLOTY 
l_EQRMAI     (T8,70H    TYPF     PIAIMF    GR  T  D     XY      ) 


4  FORMAT     (T8,15H    CONSTANTS    E     »     ,  F 1 1 . 3 » 5H, ALL     ) 

5  FORMAT     (  T  8  , 15H    CONSTANTS    G    ,     ,F11.3,5H,ALL     ) 


6    FORMAT    (T8,18HJOINT    COORDINATES     ) 
8    FORMAT     ( T8> I4»2F10.4»9H    SUPPORT     ) 


10    FORMAT     (T8,  I4»2F1C4) 

12  FORMAT  (T8,15HJOINT  RFLEASES  ) 


14  FORMAT  (T8,l4»12H  MOMENT  X  Y  ) 

15  FORMAT  (T8,I4,10H  MOMENT  X  ) 


16  FORMAT  (T8»I4»10H  MOMENT  Y  ) 

17  FOR MAT  ( T 8  » 1 9H  UNI TS> INCHES, K IPS 


WRITE  (7,3) 

WRITE  (7,4)  YOUNG 


G=YOUNG/(2.*( l.+PRATIO) ) 
WRITE  (7,5)  G 


WRITE  (7,17) 
WRITE  (7,6) 


K=l 

_Q£FINE  NU,_NJ. 


DO  100  I=1,NJ 
KI=KINDJO( I )+l 


GO  TO  (80,50,70,50) ,KI 
5  0  CONTINUE 


JOINT  IS  HINGED, RELEASE  MOMENTS 
NU(K)=NUJO(  I) 


K  =  K+1 

NUJ9=NUJO( I ) 
H0C09=H0C0R( I ) 
VEC09=VEC0R( I ) 


WRITE  (7,8)  NUJ9,H0C09,VEC09 
GO  TO  100      


70  CONTIPUE 

JOINT  IS  FIXED 
NUJ9=NUJO( I ) 

H0C09=H0C0R( I ) 


VEC09=VECOR( T ) 

WRITE  (7,8)  NUJ9,HOC09,VEC09 
GO  TO  100 
80  CONTINUE     


JOINT  IS  FREE 
NUJ9=NUJ0( I ) 
HOC09=HOCOR( I ) 
VEC09=VEC0R( I) 


WRITF  (7,10)  NUJ9,H0C09,VEC09 
100  CONTINUE  


IF  (K-l)  110,110,120 
110  CONTINUE      


. 


• 


■ 
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TJriF&E-^-RE—NO— JOINT  -RELEASES 


GO  TO  200 

. 120  CONTINUE 


THERE  ARE  K-l  JOINT  RELEASES 

WRITE  (7,12) 

L  =  K-1 

DO  150  K  =  l,[ 


NOCO=QQFDIl(NU,K) 
KI=QQFDI1 ( JQRENU,NOCO)+l 
GO  TO  (140,130,130)  ,KI 
__130  CONTINUE 


C      A  LINE  SUPPORT. RELEASE  ONLY  ONF  MOMENT 
NAJ_PHA=NUIja 


JOR=JORENU(NALPHA) 

GQ  TO  ( 131,132,140) ,JOR 


131  N9=NU(K) 

WRITE  (7,16)  N9 

GO  TO  150 
1  37  N9=NU(K) 


WRITE  (7,15)  N9 
GO  TO  150 


140  N9=NU(K) 

_ WRITE  (7,14)  N9 


150  CONTINUE 
200  CONTINUE 


201  FORMAT  (T8,19H  MEMBER  INCIDENCES  ) 
WRITE  (7,2  01_) 


202  FORMAT  (T8,4I4) 
DO  300  I=1,NUME 
NUMEM9=NUMEMB( I ) 
NUPLJ9=NUPLJ0( I ) 


NUMI J9=NUMIJ0( I ) 
__WRJ  TE  (7,202)  NUM  EM9  »  NUMI J9, NUPLJ9 
300  CONTINUE 
DESTROY  NU 


RETURN 
£JAD 


j**0_2f>_ 
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SUBROUTINE  MEMPRO 

DYNAMIC  ARRAY  HOCOR  ,VFrOK  .NU. JO.. JOKFNU.K  TNpJD.NMPl  .  ID  .  Nil  IM  T  JO  .  Nil  IMEMR 

DYNAMIC  ARRAY  K I NDMF » NUJOl » HOCOR 1 , VECOR 1 

DYNAMIC  ARRAY  MEMNU.RE I NER » TOINER 


COMMON  HOLEN»VELEN,THICK,MANUJO»YOUNG,PRATlO»HOLENl»HOLEN2»VELENl 
X>^ELFJ^2^_EJl^l.fciY_tXXNDl  ,EX2,WHY2  ,P0L0M1  »K  1 1  II  Of  I  TYPE  »  N1.IHO.J1  »NIJVE.I1  « 


2 N J, HOCOR ( P) , VECOR (P ) ,NUJO(P ) » JORENU( P) , H0LEN3 » VELEN3 »NUHO J 2 » NUH0J3 
3  »  NUVE  J2  »NUVE  J  3  tMJPI ,  J.Q (  PJ  _a.NUMI  JO  (  P  I  .NUMFMB(P)  ,KINHMF(  P)  ,  NUMF « NUVFMR 
4,NUH0MB»NURIDI »NULEDI , NUPOLO, NUJOl (P ) »HOCORl(P) »VECORK P)  »NSIZE> 
5NULDS,ENDLD,LDMSTA,LDMEND,KINDJ0(P)  ,  K  I  OUT  ,  Nl.  OTY 


DEFINE  MEMNU  »NUME 
DEFINF  RFINFR,NI.)ME 


DEFINE  T0INFR»NUME 

FORMAT  (T8»28HMEMBEfi.  PROPFRT  T  FS  %PR  TSMAT  1C 


6  FORMAT  (T8,I4»4H  IY  ,F15.5,4H  IX  ,F15.5) 
_B_..FORMAT  (T8, 1 3H  I OADING  M'  ) 


FACT0R=(THICK**3) /( 24.*(  1 .-PRAT  10**2 )  ) 
JrfRUF  (7,4) 


DO  10  I=1,NUME 
REINER( I1=*Q 


C 

10  TOINER( I )=0 

ARRAYS  SET  EQUAL  TO  ZERO 

I  l=NUHOMB+NUVEMB 
I?=I1+NURIDI 

DO  200  1  =  11  ,12 
N  =  0 

DO  100  J=l  ♦  I  T 

NALPHA  =  NUMIJO(  I  ) 


NBFTA  =  NUPLJO(  I  ) 
NALP_hL=N_U-M IJQ(  J) 


NBET=NUPLJO( J) 

IF  (NALPHA-NALPH)  100,50,100 


50  N=N+1 
MEMNU (N )=NUMFMB( J) 


100  CONTINUE 
NAO=N+1 


GO  TO  (200,200,110)  ,NAO 
CONTINUE     


DO  120  K=1,I1 
NALPH  =  NUMU„0_(  K) 


NBET  =  NUPLJO(K  ) 


IF  (NRFTA-NRFT)  I ?C 

»111 

,120 

111 

N  =  N  +  1 

MEMNU(N)=NUMEMB(K) 

120 

CONTINUE 

DOJL30  K=I7,NUME 

NALPH=NUMIJO( K ) 
NRET=NUPLJO(< ) 

NCAT=MEMNU( 1 ) 

IF  ( N  B  E  T -NUP L  J 0 ( N  CAT))  130,121,130 


121  N=N+1 

MEMNU(N)=NUMEMB(K) 


130  CONTINUE 
N  =  N  +  1 


MEMNU(N)=NUMEMB( I ) 

SIX  BEAMS  OF  PANNEL  FOUND 


VESIDE=VECOR(NBETA)-VECOR(NALPHA) 
HOS I DE=HOCOR ( NBET A ) -HOCOR (NALPHA) 


- 


: 


•        • 


(    i 


' 


tifc 


JUXtA&sSQRTt VFSIDE**2+H0  5IDE*»2) 
HOVERA=HOSIDF/VESIDE 
.niVERA  =  AD.IAG/VESIDE    


DO    180    K=1»N 
NGAMMA=MENlNUlJ<a 


GO    TO     ( 161,162,161  ,162,163,163)  ,K 
161     COE£I=LDdQ^.E&A^P.RATin/HnVFRA)*VFSTPF 


COEF2=(2*HOVERA*(  1  .+PRATIO)*<  1 . -3 .*PRAT 10 )  )*VESIDE 

GO    TO    170 , 


162  C0EF1=( l.-( HOVER A**2)*PR AT IO>*VESIDE 

COEF2=(  2.*(  l.  +  PRATIQl*l_l.-3.*PRATIO)  >*V£SIDEL 


GO  TO  170 
163  COFF1  =  (PRATIO/Hn\/FRA)*(P)TVFRA»*^n*\/F^TnF 


COEF2=0 
-U_0__R£INER  (  NGAMMA  >  =COEF1 *FACTOR+RFI NFR  (  NGAMMA  ) 


TO  I  NER ( NGAMMA )=COEF2*FACTOR+TO I NER( NGAMMA) 
JL81L  CON  T  I N  U  E 


200  CONTINUE 
DO  300  K=l .NUMF 


NUMEM9=NUMEMR(K ) 
_R£XNE9=REINER(K) 


TOINE9=TOINER(K) 
MRU  E_„  (7t6)  NUMEM9  ,RE  INE9  ,  TOI  NF9 


300  CONTINUE 

DFSTROY  MFMNU 


DESTROY  REINER 
DFSTROY  TOINFR 


WRITE  (7,8) 
RETURN 


END 


•    ' 


**076 
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SUBROUTINE  LINELO 

DYNAMIC    ARRAY     HQCQR , VFCQR ,NUJO, JORFNU , K T NDJO«NUPl  JO.IMUMI  JQ.NUJMFMB 

DYNAMIC    ARRAY    K I NDME »NU J01 » HOCOR 1 , VECOR 1 

DYNAMLC- ARRAY__ALQ£AL  s  Y  LOCAL  , XGLLD  »  YGLLD  j  NULQMF  ,  K  I  NDLO,  NULOJQ 


COMMON    H0LEN,VELEN»THICK,MAi\UJ0,Y0UNG,PRATI0»H0LENl»H0LEN2»VELENl 

1  *VFI  FN?,EX,WHY,K1ND1  ,  EX  ?  ,  WHY  2  ■>  P0L0M1  ,  K  I  L  I  LOM  T  YPE  t  N ,  J bLD ■  J 1  < NI.IVE.il  * 

2NJ»H0C0R(P  )  ,VECOR(  P)  ,NUJO(P  )  » JORENU(  P  )  , H0LEN3  » VELEN3  ».\UH0J2  ♦  NUH0J3 
3>i_NUMEJ2,NUVFJ3,NUPLJQ(P)  ,NUMI JO(P)  »NUMFMB(P)  ,  KINDMF(P)  ,  NUMF  ,  NIJVFM8 


4»NUH0MB»NURIDI  ,NULEDI  » NUPOLO,  NU  J01  (P  )  ,HOCORHP)  »VEC0R1  (P)  »NSIZE» 
5_NU  LD£  aENDLD  >_L_DMS  T  A^LDMEND_j_K  LND_JCLLPJ_t  KIQUT,NLOTY 


EQUIVALENCE     ( EX ,X A ) * ( EX 2 » XB ) » (WHY,YA) » (WHY2»YB) » ( ST ALD » P0L0M1 ) 

DFFINF    XI  OCAI.. ,NI1VFMR+NUH0MR 

DEFINE    YLOCAL     »NUVEMB+NUHOMB 

DEFINE    XGLL-D »-NU VEMB.+NUHOMB 


DEFINE    YGLLD       ,NUVEMB+NUHOMB 
.DEFINE    Ml  J I  OMF     ,NUVFMB  +  f 


DEFINE    NULOJO»NUVEMB+NUHOMB 
DEFTNE    KTNDI  0*NI  IVEMR  +  N1 IHOMR 


K  =  0 


IF    (XA-XB)    40,41,40 
40    CQH11 RUE 


41 

IF  (YA-YB)  49 
CONTINUF 

,42,49 

c 

LOADING  LINE 
00  411  1=1 ,NJ 

PARALLEL 

TO 

Y-  AXIS 

41  1 

IF  (XA-HOCOR( I ) )  411, 
CONTINUF 

412  i 

.411 

c 

LOADING  LINE 
K=0 

DOESN'T 

COINCIDE 

WITH  A 

SET 

OF 

MEMBERS 

N  =  0 
_GO_JL0    100 


412 
C 

CONTINUE 

LOADING  LINE  COINCIDES  WITH  A 

SET  OF 

MEMBERS 

K  =  0 

DO  419  J=ltN.Uy.EMB_ 

NALPHA=NUMIJO( J) 
NBETA=NUPLJO( J) 

IF  (HOCOR(NALPHA) 
413  CONTINUF 

-XA) 

419 

•  413 

,419 

414 

IF  ( VECOR( NBFTA)- 
CONTINUF 

rA) 

419, 

419,- 

414 

415 

IF  (VECOR(NALPHA) 
CONTINUE 

-YB) 

415 

,419 

,419 

IF     ( VECOR(NALPHA)-YA)     416,416,417 

)NTINUE  _ 

K  =  K+1 

XI  OCAL(K)=YA-VECOR(NALPHA) 


XGLLD(K)=YA 

GO  TO  418 
417  CONTINUE 

K=K+1   


XLOCAL(K)=0 

XGLLD(K)=VECOR(NALPHA) 
418  CONTINUE 

NULOME  (  K  )  =NUMEMB.(_J_)_ 


KINDLO(K)=2 

IF  (VECOR(NBETA)-YB)  4161 >416 1 ,4162 


■ 


■ 
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6  161  CONTINUE — 

YGLLD(K)=VECOR(NBETA) 
YLOCAL  (  K  )  =V-£CQRi-NB£TA  )-VECOR  (JiALPHAJ 


GO  TO  419 
IT  TNI  IE 


YGLLD(K)=YB 
YLOCAHK)=YB-VECOR(N ALPHA) 


419  CONTINUE 
GO  TO  200 
42  CONTINUE 
LOADING  LINE   PARALLEL  TO  X-AXIS 


DO  421  1  =  1, NJ 

TF     ( YA-VECORf T ) )     421,422,421 


421    CONTINUE 
GO    TO    U9^ 


422    CONTINUE 

IOADIN    I    INF    COINCIDFS    WITH    A    SFT    OF    MFMRFRS 

NUV=NUVEMB+1 
NUVE  =  NI.IVFMB  +  NIJHOMB 


DO    429    J=NUV,NUVE 
iiAJ_P  HA  =NU  M 1 JOLJJ 


NBETA=NUPLJO( J) 

IF  ( VFCOR(NALPHA)-YA)  479,473,479 


423  CONTINUE 
IF  (HOCOR(NRFTA)-XA)  479,479,474 


424  CONTINUE 
LE_LHOCOR(.MALPHA)-XB)  47^,479,429. 


425    CONTINUE 
LLIitOXOR  CNALP  HA)-XA)     476,476,477 


426  CONTINUE 
K  =  K  +  1 


XLOCAL ( K ) =XA-HOCOR ( N ALPHA ) 


GO  TO  428 
427  CONTINUE 


K  =  K+1 

XI OCAL(K)=0 


XGLLD(K)=HOCOR(N ALPHA) 
A2B„CQN1LNJJ£ 


NULOME(K)=NUMEMB( J) 
J£XNDLQJJLi^2 


IF     (HOCOR(NBETA)-XB)     4261,4261,4262 
4761     CONTINUF 


YGLLD(K)=HOCOR(NBETA) 

YLQC  AL  tK)  =  HOrOR(NRFTA)-HOCOR(NAI_PHA) 


GO  TO  429 
AZh2    CONTINUE! 


YGLLD(K)=XB 

Y| nCAl  <K )=XR-HOCOR(NALPHA) 


429  CONTINUE 
GQ_  TO  2D0 


49  CONTINUE 

LOADING    LINE^NOT    OF    THE_FOEM    X=CONST  .OR  ,Y=CONSXJL 

M  =  0 
K=Q 


COEF=(YB-YA)/(XB-XA) 

DO  100  I=1,NUVEMB 
NALPHA=NUMI JO ( I ) 
NBETA=NUPLJO( I ) 


!?£# 


J_E_ LN-Ol      LQk.UQ1.^Qt+ 


494    CONTINUE 

DO    493    J= 1 9 N 

IF  (NALPHA-NULOJO(N) )   493*100,493 
49  3  CONTINUE 


Y=YA+COFF*(HOCOR( NALPHA)-XA) 
TEDPL  =  ABS(VECO-R(NBETA)-Y) 


TED=(VECOR(NBETA)-VECOR(M ALPHA) )/1000. 
ABSY  =  ABS ( Y^ 


IF  (VECOR(NALPHA)-O)  491,496,491 
4  96  CONTINUE 


IF  (ABSY-TED)  57,57,491 
491  CONTINUE 


IF  (TEDPL-TFD)  56*56*492 
4-92  CONTINUE 


IF  (XA-HOCOR(NALPHA) )  50,50,100 
5  0  CONTINUE 


IF  (XB-  HOCOR(NALPHA) )  100,51,51 
51  CONTLNUE 


INTERSECTION  POSSIBLE  BUT  NOT  CERTAIN 
IF  ( VECOR(NAl  PHA)-Y)  57,100,100 


52  CONTINUE 

INTERS ECTION  STILL  POSSIBLE 


IF  ( VECOR(NBETA)-Y)  100,56,55 
55  CONTINUE 


INTERSECTION  CERTAIN 

K=K+1 


NULOME(K)=NUMEMB( I ) 
JLLQCAL  (  K  )  =Y-VECOR  (  NALPHA  ) 


XGLLD(K)=HOCOR (NALPHA) 
YG1.LD(K)=Y 


KINDLO(K)=0 
_G£L  TO    100 


56    CONTINUE 

PLUS    JOINT    OF    MEMBER    LIES    ON    LOAD    LINE 


IF     ( KINDJO(NBETA)-l )     561*562*562 
561     CONTINUF 


JOINT    NOT    A    SUPPORT 
KL=  K  + 1  


NULOME(K)=NBETA 
„XGLLD(K)=HOCOR  (MB  ETA) 


YGLLD(K)=VECOR(NBFTA) 
KINDLO(K)=l 


N  =  N  +  1 
NULOJO(N)=NBETA 


IF  (YA-YR)  100,200,100 
562  N*N+1 


NULOJO(N)=NBETA 
GO  TO  100 


57  CONTINUE 

MINUS  JOINT  LIES  ON  LOAD  LINE 

IF  (KINDJO(NALPHA)-l  )"  571,572*572 
571  CONTINUE         


JOINT  NOT  A  SUPPORT 
K=K+1 


NULOME(K)=NALPHA 
XGLLD(K)=HOCOR( NALPHA) 
YGLLD(K)=VECOR (NALPHA) 
KINDLO(K)=l 


- 


12;. 

N  =  N  +  ] . 

NULOJO(N)=NALPHA 
IF  (YA-YB)  100.70n.lQfi 
572  CONTINUE 

N  =  N  +  1  . 


NULOJO(N)=NALPHA 
.1110  CONTINUE 


C      ALL  INTERSECTION  POINTS  WITH  VERT  I CALMEMBERS  FOUNT. THERE  ARE  ALSO 

-C N  JOINT  INTERSECTIONS 

C      HORIZONTAL  MEMBERS 

00-2-00  INDEX  =  1,NUH0MB ■ ■  . 


I=INDEX+NUVEMB 
■i^L£HA=MiMI-J0  (  I  ) 


NBETA=NUPLJO( I ) 
XEL  (N-0)  130.13XXT10 


130  CONTINUE 
XHE-LIW OF  THF  I  PAD  INK  PASSFS  THRU  SOMF  JQTNTS 


_ 


DO    137    J=1,N 

TF     (NAI.PHA-NUI.Qjnt.l)  )     135,200,135 


135    CONTINUE 

_L£_    (NBFTA-NUI  0J0(  J)  )     137,700, H7 


137  CONTINUE 
M EMBER..  I  DOESNOT  PASS  THRU  A  LOADED  JOINT 


IF  (YA-YB)  138»200,139 
138  CONTINUE 


IF  ( YB-VECOR(NALPHA) )  200,1381,1381 
±1&1^.  0  NT  I N  U  E 


IF  ( YA-VECOR(NALPHA)  )  151,151,200 
139._C0J>LT_I_M11E 


IF  ( YA-VECOR(NALPHA) )  200,1391,1391 
JL391  CONTINUF 


IF(YB-VECOR(NALPHA) )  151,151,200 
1  CONTINUE 


INTERSECTION  POSSIBLE  BUT  NOT  CERTAIN 
X=XA+t ( VECOR(NALPHA)-YA)*(Xb-XA) )/( YB-YA) 


IF  (HOCOR(NALPHA)-X)  152*200,200 
1  5?  mNTINUF 


INTERSECTION  STILL  POSSIBLE 

IF  (HOCOR ( NBETA)-X)  20Q>200»L5S 


155  CONTINUE 
INTERSECTION  CERTAIN 


K  =  K+1 

NUI OMF(K)=NUMFMR( I ) 


YGLLD(K)=VECOR(NALPHA) 
_XGU_D(K)=X 


XLOCAL(K)=X-HOCOR (N ALPHA) 
■  KINPLQCKJ^O 


200  CONTINUE 
N1)LDS  =  K 


CALL  LINE  ($KlNDLO,$XGLLD,5»YGLLD,$NULOME»iXLOCAL,$YLOCAL) 
DESTROY  XLOCAJ 


DESTROY  YLOCAL 
DESTROY  XGLLD 


DESTROY    YGLLD 
HF^TRnv    Mill  OMF 


DFSTROY    KINDLO 
DESTROY    NULQJXL 
RETURN 

END 


: 


J  2*. 


**  0-2-6- 


SUBROUTINE  L I NE ( K I NDLO,XGLLD» YGLLD >NULOME »X LOCAL ♦ Y LOCAL ) 

aY-NAMlC  ARRAY  HOrOR  ,  VfCOR  %NI  l.JO  .  .JDRENI 

DYNAMIC  ARRAY  K I  NOME , NU JOl » HOCOR 1 , VECOR 1 

DYNAMIC  ARRAY  XLOCAL  ,  Y  LOCAL  ,XGLLD»YGLJ_,NUl_M£_JLLNDLO 


DYNAMIC  ARRAY  RAD  I  US » STALME , P , JO  I  NT 

COMMON  HOI  FN.VFi  FM  .  TH  T  CK  ,  M  AN1I  I.  ID  .  YDI  IN  C  .  PR  AT  I  O  .  HOI  FNl.Hni  FN7.\/FIFM1 

1»VELEN2»EX,WHY,KIND1 » EX2 » WHY 2 tPOLOMl > K I L I LO ,LTYPE »NUHO Jl ,NUVE Jl » 
___  , HOCOR  LP)  ,  VECOR  (  P)  ,NUJQ  (  P  >  >  JORFN'J{  P  )  ,HOLEN3  .  VFL  FN3  .NUHC  J7  ,  NUHQJ3 
3,NUVEJ2»NUVEJ3»NUPLJO(P) »NUMI JO( P) ,NUMEMB( P) , KINDME( P) ,NUME»NUVEM6 
4,NUH0MB,NURIDI ,NUIEDI ,  NUPOLO, NU JQ1 (P) , HOCOR 1 (P) ,VFCOR1 (P) ,NSI7F, 


5NULDS,ENDLD»LDMSTA,LDMEND,KINDJO( P) , KIOUT»NLOTY 

FQ.UIVAI  FNfF  (FX,XA),(FX7,XR),(WHY,YA)  .  (WHY7«YR)  «  (  £T  Al  n  .  PHI  OM1  ) 

C      THIS  SUBROUTINE  IS  A  LOGICAL  CONTINUATION  OF  LINELO.IT  WAS  FORMED 

X AS  A  SEPARATE  PROGRAM  SINfF  PRFCOMPIlFR  AT  THIS  TTMF  CANNOT  ATfFPT 

C      SUCH  A  LARGE  NUMBER  OF  SYMBOLS 

C SHORT  LIST  OF  PHLNJ.S ACCORDING  TO  THFIR  DISTANCF  FROM  THE 

C      LEADING  EDGE  OF  THE  LINE  LOAD 
DFFTNF  RAPTUS  .  Nl  IWFMR  +  NI  IHDMR 

DEFINE  STALME  »NUVEMB+NUHOMB 
HE  FINE  __P_     .  «JSUiVEMfi±NUH OMB 

DEFINE  JOINT   ,NUVEMB+NUHOMB 
10  FORMAT  (T8,12H  JOINT  LOADS) 


11  FORMAT  (T8,I4,9H  FORCE  Z  ,F15.3) 
17  FORMAT  (T8,I4,10H  MOMFNT  X  .FH.^) 


13  FORMAT  (T8,I4»10H  MOMENT  Y  ,F15.3 

14  FORMAT  (T8».13-hL  MFMBEf_JLClAJD__ 


15  FORMAT  (T8,I4»23H  FORCE  Z  CONCENTRATED  P,F15.3»2H  L,F15.3) 

16  FORMAT  (T8,I4 , 24H  MOMENT  X  CONCENTRATED  P,F15.3,2H  L,F15«3) 

17  FORMAT  (T8,I4»24H  MOMENT  Y  CONCENTRATED  P»F15.3,2H  L.F15.3) 

18  FORMAT  (T8,I4,19H  FORCE  Z  UNIFORM  W  ,F15.3,4H  LA  ,F15.3,3H  -  / 
1T8,4H  LB  »F15.3) 

19  FORMAT  J T 8,I4,20H  MOMENT  X  UNIFORM  W  ,F15.3,4H  LA  ,F15.3,3H  -  / 


1T8,4H  LB  ,F15.3) 
20  FO R  M  A  T  (  T  8  ,  I  4  »  2 OH  MOMENT  Y  UNIFORM . W  ,F15.3,4H  LA  ,F15.3,3H  -  / 


1T8,4H  LB  ,F15.3) 

21  FORMAT  (T8,I4,19H  FORCE  Z  LINFAR  WA  , F 1 5 . 3 , 4H  WB  ,F15.3,3H  -  / 
1T8»4H  LA  ,F15.3,4H  LB  ,F15.3) 

22  F  ORMAT  (  T  8  ,.  I  4  .  2  0  H__MOMEN  T  X  J_ I  N  FAR  WA  ,F15.3,4H  WB  ,F15.3,3H  -  / 
1T8,4H  LA  »F15.3»4H  LB  »F15.3) 

2  3  FORMAT  ( T  8 , I  4  »  20H  MOMENT  Y  LINEAR  WA  ,F15.3,4H  WB  ,F15.3,3H  -  / 


1T8*4H  LA  ,F15.3»4H  LB  ,F15.3) 
74  FORMAT  (T8,11H  I  1ST  LOADS  ,14) 


25  FORMAT  (T8,20H  LIST  DISPLACEMENTS  ,14) 
DO  201  I=1,NULDS 


IF  (KINDLO( I )-l )  2009,2009,2006 
7006  CO N T I N U E 


IF  (XA-XB)  2008,2007,2008 
2007  CONTINUE 


RADIUS( I )=XGLLD( I )-YA 
GO  TO  2004 


2008  CONTINUE 

RADIUS( I )=XGLLD( I )-XA 


GO  TO  2004 

2009  CONTINUE 

RAD  I  US (  I  )=SQRT(  (XGLLDf I ) -XA ) **2+ ( YGLLD ( I  )-YA)**2 ) 

2004  CONTINUE 


IF  (LTYPE-1)  2002»2003,2002 
2002  CONTINUE 


IZj. 


JUUUJ\L    LOADING    UNIFORM 


2003 


STALME( I )=STALD 

GO- TO  201 

CONTINUE 

AFN  =  SORT(  (  XB-X  A  )  **2-*-(  YB-Y A)  **2  ) 


INITIAL  LOADING  IS  LINEAR 
-SXALMFf  I)=STALn+(RADIMSCI)/AEN)*<  ENDLD-STALD) 


201 


CONTINUE 
IF  (  KINDLI 
2005  CONTINUE 

DO  210  J=1,NULDS 


_2-0J)5v2O05»35  0 


DO  210 
T=L--1- 


L=2*NULDS 


.2  02 


IF  (RADIUS(D-RADIUS(  I  )  ) 

CONTINUE  


202*210,210 


NULOM=NULOME( I ) 
XGLL  =  XGLLD11) 


YGLL=YGLLD( I ) 
XLOCA  =  XLOCAL(  I  ) 


KINDL=KINDLO(  I  ) 
-R  AD  I  U= RADIUS  (  I  ) 


STALM  =  STALME(  I  ) 
miLDM EJJJ^NilLD MEM   ) 


XGLLD( I )=XGLLD(L) 
YGI  I   D(  T  )=YG1  I   D(l    ) 


XLOCALf I )=XL0CAL(L) 
JUND LQJ  I  1=K LND  10(1) 


RADIUS! I )=RADIUS(L) 
_SJLAL_M£iJLl^S  T  A  L  MF(L) 


NUL0ME(L ) =NULOM 
■XnLLDd   )=XGI  1 


YGLLD(L)=YGLL 
XLOCAL(L)=XLOCA 
KINDLO(D=KINDL 
.RADIUS  (  L  )  =RADIU 


STALME(L)=STALM 
210    CONTINUE 


LIST    SHORTED 
2101     CONTINUE 
GO 


_21.QL3_ 


TO    (2102»2  102»2102,2103»2102»2102>2103) »KI0UT 
CONTINUE 


DO    2104     I=1,NULDS 
JUNK=KTNDI 0( I )-t-1 


2105 


GO  TO  (2105*2106*2107)  »JUNK 
rQJMT  INUE 


LINE  INTERSECTS  A 
MEMBER  =  NULOME(  I  ) 


MEMBER 


NGAMMA=NUMIJO( MEMBER) 

NDFI  TA  =  NUPI  JO(MFMRFR) 

RLEME=SQRT( ( HOCOR ( NDELTA ) -HOCOR ( NGAMMA ) ) **2  + ( VECOR ( NDELTA ) -VECOR ( 
JLNJaAMMA )  1  **T  ' 


IF 


.21Q, 


(XLOCAL( I )-RLEME/2.) 
!T(1  )=NUJO( NGAMMA) 


2108,2108*2109 


NADD=0 

GO  TO  2104 


2109  JOINT( I )=NUJO(NDELTA) 

NADD=Q 


GO  TO  2104 
.2106  CONTINUE 


126. 


■  IQTNT  (  I  )=N1II  DME(  I  ) 


NADD=0 
-  -GO-  TO  21-CH— 
2107  CONTINUE 

M  EMBER  =  NU  .LOMF-t-U- 


NGAMMA=NUMI JO ( MEMBER ) 
MDELTA-MUPL JO ( MEMBER ) 


JOINT ( I )=NUJO(NGAMMA) 
-JOINT  C  I  +  ll  =  NLUa 
NADD=1 
2104  CONTINUE   


NUL=NULDS+NADD 
-LF-1KI0UT-4)  2110*2111*2110 


2110  JOIN9=JOINT( 1 ) 

WRITE  (7,-25)  ^J0IN9 
DO  2112  I=2»NUL 

LE^LJJIIMJ 


2±±1±2U  2*2113 


2113  J0IN9=J0INT( I ) 
WRTTF  (7,75)  JOTN9 


2112  CONTINUE 
G_Q  Ta  3  60 


2111 


J0IN9=J0INT(1 ) 
WRITE  (7*24)  JOIN9 


DO  2114  I=2»NUL 

IF  ( JOINT( I )-JOINT( 1-1 ) ) 


2115*2114*2115 


2115  JOIN9=JOINT( I ) 
WRITE  (7*24)  JOIN9 


2114  CONTINUE 
GO  T0_3  60 


2102    CONTINUE 
TF    (NUI  D.Szil 


)  ?1?*711 ,71? 


211  CONTINUE 

^MLY,QNE  CONCENTRATED. 


J_QAJX 


217 


FACTOR=SQRT( (XA-XB)**2+( YA-YB)**2) 
J2.i_ii„=F_A„C  T OR* ( (  ST  A I  D+ENDLD)/2.  ) 
GO    TO    300 

XOMULNJJJZ 


MORE  THAN  ONE  CONCENTRATED  LOADS 

FACTOR=SQRT ( (XA-XGLLD( 1) ) **2  +  ( YA-YGLLD( 1) ) **2 ) 

P(l)=FACTOR*((STALD  +STALME ( 1 ) )/2. ) 

DO  213  L=2»NULDS 


I  =  L-1 

FACTOR=RADIUS(L)-RADIUS(  I  ) 

C=FACTOR*(2.*ABS(STALME(L) ) +ABS ( STALME ( I ) ) )/(3.*( ABS(STALME(L ) )+ 
„1AB_S(  STALME (  I  )  )  )  ) 


P(L)=C*(  (STALME(D+STALME(  I  )  )  /2.  ) 
_P  (  I  )  =P  (  I  )  +F  ACTOR*  ( ...(  STALME  (L>+STALME(  I  )  )/2»  )-P(L) 


213  CONTINUE 

FACTOR=SQRT( ( YB-YGLLD ( NULDS ) ) **2+( XB-XGLLD ( NULDS ) )**2) 


P(NULDS)=P(NULDS )+F ACTOR* ( ( ENDLD 
300  CONTINUE 
TH^RE 
IOP  =  0 


+STALME(NULDS>  )/2.) 


ARE  NULDS  CONCENTRATED  LOADS 


DO  320  K=l, NULDS 
IF  (KINDLO(K)-l ) 


370*301*301 


301  CONTINUE 

JOINT  LOADS 


3011 


IF  (IOP-1)  3011,3012*3012 
WRITE  (7,10) 


127. 


30  12  IOP=IOP+l 


GO    TO     (302,303,304)  ,NLOTY 
-302    N  U  L  OM9 =  N1LLQM  F  (  K  ) 
P9  =  P( K) 


GO    TO    320 
303    NU-LOM9  =  Nl.ILQM£XtU- 


P9=P(K) 

WRITE    (  7  ♦  1.2  )    WJLQm  ,P9 


GO    TO    320 

.04    NULQR9=NU LOMEIK) 

P9=P( K) 
WRITE     (7,13)     NULQM9,P9 


20    CONTINUE 
OCONJ—LOADS  .XlaiU— M£MB£5-L£)AaS-N0W 


IOP  =  0 

DO    340    K=l 

•  Mil  ns 

321 

IF     (KINDLO(K)-I ) 
rnNTTN|iiF 

321 

,340,340 

c 

MEMBER    CONCENTRATED 
TF     (TOP-1)     3211»3212 

LOADS 
,321? 

3211    WRITE    (7,14) 
1QPJE.J.QP_+] 


322 

GO    TO    (322,323,324) 

NUI  0M9  =  MJI  OMF(K  ) 

,NLOTY 

P9=P(K) 

XI  OCA9=XI  OfAl    IK) 

WRITE     (7,15)     NULOM9 
GO    TO    340 

,P9 

>XLOCA9 

323 

NUL0M9=NUL0ME(K) 
P9=P(K) 

XL0CA9=XL0CAL( K) 
MRU F    t7«1 6.)    NUL0M9,P9,XL0CA9 


GO    TO    340 

MJJ_  0  M  9_=Mi  LQMEJJO. 

P9=P(K) 

XI  OfA9  =  XI  OfAl    (K) 


WRITE     (7,17)     NULOM9,P9,XLOCA9 
lONT  I  NUE 


THIS    LOADING    ARRANGED. TRY    NEXT 
GO    TO    370 


50    CONTINUE 

CD    TO  (  3501  ,?501  ,3501  •2103,  3501  ,?501  ,7103,3501  )   •KTOUT 


35 


01  CONTINUE 

MEMBER  LINE    LjQADS 

NULD=NULDS+1 

ST  ALME  (  NULD  )  =£NDLD_ 


IOP  =  0 

DO  360  l  =1  .NUI  PS 


I=L  +  1 

IF  (  IOP-1)  3.5 1.»352»  35.2 


3 

3 


51  WRITE  (7,14) 
52_ I0P=I0P+1 


GO  TO  ( 353,354,355)  ,NLOTY 
5  3  CONT-IMIE 


35 


IF     (LTYPE-1)     3531  ,3532, 3531 
31     NUL0M9  =  NUL0ME1U 


STALM9=STALME(L) 

XL0CA9=XL( 


128, 


3532 


yi  nrAQ  =  Yi  nrAi  t\  \ 
WRITE    (7,18) 
GO    TO    360 
NULOM9= 
STALM9= 
STALM8= 
JCUKA9  = 


NUL0M9,STALM9,XL0CA9,YL0CA9 


NULOMF(L) 
STALME(L) 


STALMF(  I  ) 
XLOCAL(L) 


YLOCA9= 
WRITE  { 
GO  TO  3 


Jl5_4 


CONTINU 
IF     (LTY 
-3541     NI.IL0M9= 


S  T  ALM  9 , S TALM8  ,  X LOCA9  >  YLQCA9^ 


PE-1)     3541,3542,3541 
NULOMEi-U 


STALM9= 
JCL£>CA9  = 


YL0CA9= 
WRITE  ( 


STALMF(L) 
XJ_QCAL(L) 
YLOCAL(L) 
Zjl19)     Nl 


, S  T  ALM 9  ,  X  LQCA9_^LOCA2_ 


GO    TO    3 
-15  4  2    NULQMSLs 


60 

Mill  OME(  L) 


STALM9= 
STALM8= 
XLOCA9= 
YL0CA9= 
WRITE  ( 
GO  TO  3 


STALME(L) 

STALMEt I ) 


XLOCAL(L) 

YLOCAKU 

7,22 )     NULOM9,STALM9,STALM8,XLOCA9,YLOCA9 
bQ 


355 


3551 


CONTINU 

rr  (lty 

NULOM9= 
STALM9= 


E 
PE-1J_ 


3  5  5 1 ,3552,3551 


XL0CA9= 
YI  0CA9^ 


NULOMEf  L) 
SIALME(L) 
XLOCAL(L) 
YJ-OCAI   (I  ) 


WRITE 
_GQ.  TO 


7,20) 
60 


NUL0M9,STALM9,XL0CA9,YL0CA9 


3552 


NUL0M9= 

STALM9= 

STALM8= 

_XJ_Q£A9je 


NULOME(L) 
STALMFt L) 

STALMF( I ) 
XI  OCAI  (I  ) 


YLOCA9=YLOCAL( L> 

W R I IE_J  7,23)  N U L 0 M 9  ,STALM9 ,ST Al  M8 , XLOCA9 , YLOCA9 


360 

3  70 


CONTINUI 
CONTINU 
DESTROY 
DFSTROY 


RADIUS 
STAI  ME 


DESTROY 
DESTROY 
RETURN 
EMD 


JOINT 
JP_ 


**D.Z£l. 


SUBROUTINE    AREALO 

DYNAMIC    ARRAY     HOCQR  ,  V  FCQR  ,  NU  JO  ,  JOR  FNIU  K  T  NpJO  «  Nl  IP!   .  in*  Nl  IM  T  JO  .  Nl  IMFMR 

DYNAMIC    ARRAY    K I NDME , NU JOl » HOCOR 1 > VECOR 1 

DYNAMIC    ARRAY     YMFMB  »  XST  A M  F  »  ST  A|  DM.FNpl  DM 


COMMON    HOLEN,VELEN,THICK,MANUJO»YOUNG,PRATlO,HULENl,HOLEN2,VELENl 
>i/EJ-EN2^FJL^HY,KIND1  ,  FX?  »  WHY  ?  ,  POI  OM1 «KII    II  0*1   T  YPF  .  NUHO  J  1  .N1IVFJ1  , 
2NJ,hOCOR(P)  ,VECOR(P)  .NUJO(P)  , JORENU(P)  »H0LEN3 »VELEN3 »NUH0J2 »NUH0J3 
3  ,  NUVE J2 ,  NUVE  J3  j.NUPLJO  (  P)  »NUM  I  JOl  P  !  «NUMFMR(P)  «KT  NDMEiPJ  »_NUME  »MJV  EMB_ 
4,NUHCMB,NURIDI jNULEDI »NUPOLO, NUJOl IP ) ,H0C0R1(P) , VECOR 1 (P) »NSIZE» 
.-5NULDS»ENDLD»LDMSTA»LDMEND»KlNDJ0(PJ_i.janUI*ML0TY 


REAL    LDMSTA,LDMEND 

FQUTVAI   FNCF     (FX,XA) *  (WHY»YA)  «(FX?.XR1.(WHY?.YR1.(<;T4I  p.  POI  DM]   ) 

DEFINE    YMEM3»NUH0MB+NUVEMB 

DEFINE    XSJ\AME^NUHQM£+NUVEMJ3 


DEFINE    STALDM,NUHOMB+NUVEMB 
D £  FINE    EN  DL  DM_aNUHOM£+NJJ  V  E  M  B 


IF     (LTYPE-1)     10,10,20 
1Q    CONTINUF 


LOADING    UNIFORM    OR    CHANGING    WITH    X. DISTRIBUTE    AMONG    HOR IZ . MEMBERS. 
K  =  0 


DO    100     INDEX=1 ,NUH0MB 
I=INDFX+NUVFMR 
NALPHA=NUMIJO( I ) 
NRFTA=NUPI  JOt  T  ) 


IF     (XA-HOCOR(NBETA) )     11,100,100 
IUEL 


IF     (XB-HOCOR(NALPHA) )     100,100,12 
,1?    CONjLLMUE 


_ 


OVFRLAP    POSSIBLE    BUT    NOT    CERTAIN 
TF     (YA-\/FfnRfNRFTA))ntn.inO 


13    CONTINUE 
OVERLAP    SXILL    POSSIRLF 


IF     ( YB-VECOR(NALPHA) )     100,14,14 
14    CONTINUF       


OVERLAP    CERTAIN    , RECORD    Y. 
K  =  K+1 


YMFMB(K)=VFCOR(NBETA) 
XSTAME  (  K  )  =HOCQR  (  NALPMM. 


100    CONTINUE 
NaT 


_ 


DO    110     I=1,K 
I  =T  +  1 


IF     (XSTAME(D-XSTAME(  I  )  )     129,108,129 

.  LQ  8    H  =  N  + 1 


110  CONTINUE 

1 79  ADFNO=YB-YA 

Y1  =  YA 

Y2=YR 


130  NUMBER=N 

IF  (  N  UMBE  R  - 1 J^L3:L»JL12  a  1 3 1 


132  CONTINUE 

^STALDM ( 1 ) =LDMSTA*ADENQ 
ENDLDMf 1 ) =LDMEND*ADFNO 
GO  TO  ISO 


131  CONTINUE 

DO  15lQ  I  =  1,N 


L=I  +  1 

m=_l^-:l 


130, 

TF     (T-l)     134»133tl3/i 

133  CONTINUE 

STALDMC  !)=(((  YMEMB(  I  ) -Y 1  )  +A  YMEMB  (  L  )  --YJHEMB  M  )  1./2-*.)  »LDMSTA) 
ENDLDMt I )=( ( ( YMEMB( I ) -Y 1 )+( YMEMB(L )-YMEMB( I  )  ) /2 . ) *LDMEND ) 

GO  TO  150 ; . 


134  CONTINUE 
IF  (I-N)  135*137*135 


135  CONTINUE 

STALDM(  I  )  =  (  ( YMFMBt  L  )  -YJ4EMBLM  L)  72.  )  *LDMSXA_ 
ENDLDM( I ) =( (YMEMB( L)-YMEMB(M) )/2. )*LDMEND 

GO  TO  1,50 


137  CONTINUE 

SXALDNK  I )  =  (  ( YMEMB(  I  )-YMEMB( M)  )  /2 .  + ( Y2-YMEMB (  I )  )  )*LDMSTA 

ENDLDMt  I  )  =  (  < YMEMB(  I  )-YMEMB(  M)  )/2.  +  (  Y2-YMEMB(  I  )  )  )*I_DMEND 
-15D--CON-T  I  N  U  E  


LOADS  READY  TO  CALL  LINELO 
IF  (  L  T  Y  PF-1  )  3fU?0.  5_(L 


30  CONTINUE 
UNIFORM  I  PAD  OR  L T NEAB_ Al  ONG  X. 


DO  31  I=1»NUWBER 
YA=YMEMB( I ) 


YB=YMEMB( I ) 

S  T  A  L  D  =  STALDM(  I ) 


ENDLD=ENDLDM( I ) 
CALL  LINFLO 


31  CONTINUE 
_GO  TO  300 


50  CONTINUE 
LINEAR  ALONG  Y-AXIS. 


DO  51  I=1»NUMBER 
XA  =  YMFMB(  I  ) 


XB=YMEMB( I ) 
^TALD=STALDM( I ) 
ENDLD=ENDLDM( I ) 
L  T  YPE=1 


CALL  LINELO 
51  CONTINUE 


GO  TO  300 
20  CONTINUE 


LOADING  LINEAR  CHANGING  WITH  Y. DISTRIBUTE  AMONG  VERT. MEMBERS 
K  =  0  


DO  200  I=1,NUVEMB 
_NAJLP.HA  =  NUMIJO(  I  ) 


NBETA=NUPLJO( I ) 

IFJYA -VFCOR(NBETA) )  21*200*200 


21  CONTINUE 
IF  (YB    -VFrOR(NALPHA) )   200>200>22 


22  CONTINUE 

OVERLAP  POSSIBLE  BUT  NOT  CERTAIN 


IF  (XA    -HOCOR(NBETA) )  23*23*200 

23  CONTINUE  


OVFRLAP  STILL  POSSIBLF 

IF  (XB    -HOCOR(NALPHA)  )_200_»24»24 


24  CONTINUE 

OVERLAP  CERTAIN  RECORD  X. 


K  =  K  +  1 
YMEMB(K)=HOCOR(NBETA) 


XSTAME(K)=VECOR(N ALPHA) 
200  CONTINUE      


N  =  Q . 

DO  210  1=1, K 

IF  (XSTAMEt  I  )-XSTAMF(  1)  )  2  1 0  ,  208^2011 
208  N  =  N  +  1 

XSTAMF(N)=XSTAMF< I ) 
YMEMB(N)=YMEMB( I ) 
210  XONT  I  NUE 


229  ADENO=XB-XA 

Y1=XA 

Y2  =  XB 
GO  TO  130 


300  CONTINUE 

BI-SXROY— ¥  HEMS- 
DESTROY  XSTAME 
DESTROY  STALDM 
DESTROY  FNDLDM 
RETURN 
END 


026 


1&. 


SUBROUTINE  OUT 

DYNAMIC    ARRAY    HOCOR,VE<         ,      JJO,JC  .      [IN    JO  ,NUPLJU ,NUMI JO ,NI 

RAY    K .  _  -  .  .  OR1 

•  •       NUJO, YOUNG, PKATIO,HOLENl,HOLEN2»VELENl 

1»VELI    .     •        ,WHY»KIN     1,EX2,      IY2-,  ■  •      ryi       •         IOJ1  »(MUVEJ  1 » 

2NJ,        ZOR(P.)  tVECOR(p)  tNUJO(P)  tJOl  '  >i-iULEN3  , VELEN3  »NUHOJ2  ,N       OJ3 

3»NUVEJ2»NUVEJ3»NUPLJO(P)  ,NUMIJU(I    ■•  I    (  P  )  »K  IlMuMt  (  P  )  ,NUMt  ,NUVt 

4»N  -  |  .    ..      JL     DI,NUPOl      »IM    J01(P)>H  L  ( P  )  »  VECORl  ( P  )  ,NSIZE  , 

NULDS»ENDLD>LDMSTA,LC  •     -         JO(P),  -     LoTY 

1  FORMAT  (T8,15H  LIST  REACTIONS) 

2  FORMAT  (T8»11H  LIST  LOADS) 

3  FORMAT  (T8»19H  LIST  DISPLACEMENTS) 
A  FORMAT  (T8,20H  STIFFNESS  ANALYSIS 

6  FORMAT  (TS,17H  LOADING  LIST  »1») 

IF  (MANUJO)  1U,5»10 
5  CONTIi 

WRITE  (7,6) 
;  I  T  E  (7,4) 

MANUJO=MANUJO+l 
10  CONTINUE 

GO  TO  (2C,   •    •  ,    •    >   ,KI0UT 

20  CONTINUE 

REACTIONS  WANTE 
•  I T  E  (7,1) 

GO  TC  1  . 
30  CONTINUE 

POINT  LOADS  OR  DlSPL,  ITED 

NUPOLO=2 

LOADMA=0. 

CALL  POINLO 

GO  TO  1- 
40  CONTINUE 

LINE  LOADS  WANTED 

NLOTY=l 

CALL  LINELO 

GO  TO  1C 
50  CONTINUE 

ALL  LOADS  WANTED 

WRITE  (7,2) 

GO  TO  100C 
80  CONTINUE 

ALL  DISPLACEMENTS  WANTED 

WRITE  (7,3) 
1000  CONTINUE 

RETURN 

END 


• 
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SAMPI F    PRORI  FM-DATA 


/DATA 


90. 
65 


90. 


3000 


1 

4 

1 

1 

0.0 
90. 

0.0 
0.0 

0.0 
90. 

90.0 
90. 

2 
? 

0.0 
0.0 

0.0 
90. 

90.0 
90. 

0.0 
90.0 

1 

1 

0.0 

0^,0 

0.0 

1 

90. 

90. 

-10.0 

7 

7 

3 
4 

10. 
20. 

15. 
25. 

20. 

45. 

5 
6 

60. 

10. 

7 
8 

60. 

10. 

80. 

30. 

**EOF 
/END    OF 

FILE 

SAMPLE  PROBLEM. 

SQUARE  PLATE  90.X90.  INCHES. 
THICKNESS=0.5  INCHES. 
TWO  CLAMPED, TWO  HINGED  EDGES. 
YOUNG'S  MODULUS=3000O.  KIPS 
POISSON'S  RATIO=0.3 
HYDROSTATIC  LOADING. 
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TYPE 

PLANE  GRID 

XY 

CONSTANTS  E  , 

3000U.000,ALL 

CONSTANTS  G  , 

11538. 465, ALL 

UNITS 

-♦INCHES, KIPS 

JOINT 

COORDINATES 

1 

0.0 

0.0 

SUPPORT 

2 

0.0 

12.8571 

SUPPORT 

3 

0.0 

25.7143 

SUPPORT 

4 

0.0 

38.5714 

SUPPORT 

5 

0.0 

51.4286 

SUPPORT 

6 

0.0 

64.2857 

SUPPORT 

7 

0.0 

77.1429 

SUPPORT 

8 

0.0 

90.0UU0 

SUPPORT 

9 

12.8571 

0.0 

SUPPORT 

10 

12.8571 

12.8571 

11 

12.8571 

25.7143 

12 

12.8571 

38.5714 

13 

12.8571 

51.4286 

14 

12.8571 

64.2857 

15 

12.8571 

^7. 1429 

16 

12.8571 

90.0000 

SUPPORT 

17 

25.7143 

0.0 

SUPPORT 

18 

25.7143 

12.8571 

19 

25.7143 

25.7143 

20 

25.7143 

38.5714 

21 

25.7143 

51  .4286 

22 

25.7143 

64.2857 

23 

25.7143 

77.1429 

24 

25.7143 

90.0000 

SUPPORT 

25 

38.5714 

0.0 

SUPPORT 

26 

38.5714 

12.8571 

27 

38.5714 

25.7143 

28 

38.5714 

38.5714 

29 

38.5714 

51.4286 

3  0 

^0.5714 

6h.^857 

31 

38.5714 

77.1429 

32 

38.5714 

90.0ou0 

SUPPORT 

33 

51.4286 

0.0 

SUPPORT 

34 

51.4286 

12.8571 

35 

51.4286 

25.7143 

36 

51.4286 

38.5714 

37 

51.4286 

51.4286 

38 

51.4286 

64.2857 

39 

51.4286 

77.1429 

1 

• 

* 

r                                                        -                         —     , 

• 

• 

- 

' 

• 

, 

■ 

■» 

. 

. 

, 

, 

» 

• 

i??. 


40 

5  1  .4286 

90.0000 

SUPPORT 

41 

64.2857 

.0 

SUPPORT 

42 

64.2857 

12.8571 

43 

64.2857 

25.7143 

44 

64.28  5  7 

38.5714 

4b 

64.2857 

51.4286 

46 

64.2857 

64.2057 

47 

64.2857 

77.1429 

48 

64.2857 

90.0000 

SUPPORT 

49 

77.1429 

0.0 

SUPPORT 

50 

77.1429 

12.8571 

51 

77.1429 

25.7143 

52 

77.1429 

38.5714 

53 

77.1429 

51.4286 

54 

77.1429 

64.2857 

55 

77.1429 

77.1429 

56 

77.1429 

90.0000 

SUPPORT 

57 

90.0000 

0.0 

ojPPuRT 

58 

90.0000 

12.8571 

SUPPORT 

59 

90. Ou 

25.7143 

SUPPORT 

6.; 

90.  Ol 

38.5714 

SUPPORT 

61 

90.00 

51.4286 

SUPPORT 

62 

90.0000 

64.2857 

SUPPORT 

63 

90.00u0 

77.1429 

SUPPORT 

64 

90.0000 

90.0o00 

SUPPORT 

JOIN! 

'  RELEASES 

2 

MOMENT  Y 

3 

MOMENT  Y 

4 

MOMENT  Y 

5 

MOMENT  Y 

6 

MOMENT  Y 

7 

MOMENT  Y 

58 

MOMENT  Y 

59 

MOMENT  Y 

60 

MOMENT  Y 

61 

MOMENT  Y 

62 

MOMENT  Y 

63 

MOMENT  Y 

MEMBER  INCIDENCES 

1 

1    2 

2 

2    3 

3 

3    4 

4 

4    5 

5 

5    6 

6 

6    7 

7 

7    8 

8 

9   1-G 

9 

10   11 

10 

11   12 

11 

12   13 

12 

13   14 

13 

14   15 

14 

15   16 

15 

17   18 

16 

18   19 

■  -    • 


17 

19 

18 

20 

21 

19 

21 

22 

20 

22 

23 

21 

23 

24 

22 

25 

26 

23 

26 

27 

24 

Zl 

28 

25 

28 

29 

26 

29 

3u 

27 

3u 

31 

28 

31 

32 

29 

33 

34 

30 

34 

35 

31 

35 

36 

32 

36 

37 

33 

37 

38 

34 

38 

39 

35 

39 

40 

36 

41 

42 

37 

42 

43 

38 

43 

44 

39 

44 

45 

4  0 

45 

46 

41 

46 

47 

42 

47 

48 

43 

49 

5  0 

44 

50 

51 

45 

51 

52 

46 

52 

53 

47 

53 

54 

48 

54 

55 

49 

55 

56 

50 

57 

58 

51 

58 

59 

52 

59 

6o 

53 

6  0 

61 

54 

61 

62 

55 

62 

63 

56 

63 

64 

57 

1 

9 

58 

2 

lo 

59 

3 

11 

6( 

4 

12 

61 

5 

13 

62 

6 

14 

63 

7 

15 

64 

8 

16 

65 

9 

17 

66 

10 

18 

67 

11 

19 

68 

12 

20 

69 

13 

21 

70 

14 

22 

71 

15 

23 

vy« 


72 

16 

24 

73 

17 

25 

74 

18 

26 

75 

i  ) 

27 

76 

20 

28 

77 

21 

29 

78 

22 

3  J 

79 

23 

31 

24 

32 

81 

25 

33 

82 

26 

34 

83 

27 

35 

84 

28 

36 

85 

29 

37 

86 

3o 

38 

6  7 

31 

39 

88 

32 

40 

89 

33 

41 

90 

34 

42 

91 

35 

43 

92 

36 

44 

93 

37 

45 

94 

38 

46 

95 

39 

47 

96 

4  0 

48 

97 

41 

49 

98 

42 

5  u 

99 

43 

51 

10  0 

44 

52 

101 

45 

53 

102 

46 

54 

103 

47 

55 

104 

48 

56 

105 

49 

57 

106 

50 

58 

107 

51 

59 

108 

52 

6  0 

109 

53 

61 

110 

54 

62 

111 

55 

63 

112 

56 

64 

113 

1 

1  J 

114 

2 

11 

115 

3 

12 

116 

4 

13 

117 

5 

14 

118 

6 

15 

119 

7 

16 

1  2  0 

9 

18 

121 

1C 

19 

122 

11 

2C 

123 

12 

21 

124 

13 

22 

125 

14 

23 

126 

15 

24 

140. 


127 

17 

26 

128 

18 

27 

129 

19 

28 

130 

20 

29 

131 

21 

3U 

132 

22 

31 

133 

23 

32 

134 

25 

34 

135 

26 

35 

136 

27 

36 

137 

28 

37 

138 

29 

38 

139 

30 

39 

31 

40 

141 

33 

^z 

142 

34 

43 

143 

35 

44 

144 

36 

45 

145 

3  7 

46 

146 

38 

47 

147 

39 

48 

148 

41 

50 

149 

42 

51 

15. 

43 

52 

151 

44 

53 

152 

45 

54 

153 

46 

55 

154 

47 

56 

155 

49 

58 

156 

50 

59 

157 

51 

60 

158 

52 

61 

159 

5  3 

62 

15 

54 

63 

161 

55 

64 

162 

9 

2 

163 

10 

3 

164 

11 

4 

165 

12 

5 

166 

13 

6 

167 

14 

7 

168 

15 

8 

169 

17 

10 

170 

18 

11 

171 

19 

12 

172 

2  0 

13 

173 

21 

14 

174 

22 

15 

175 

23 

16 

176 

25 

18 

177 

26 

19 

178 

27 

2  0 

179 

2d 

21 

180 

29 

22 

181 

30 

23 

141. 


182 

31 

24 

183 

33 

26 

184 

3  4 

27 

185 

35 

28 

186 

36 

29 

187 

37 

3j 

188 

38 

31 

189 

39 

32 

19C 

41 

34 

191 

42 

35 

19? 

43 

36 

193 

44 

3  7  ' 

194 

45 

38 

195 

46 

39 

196 

47 

40 

197 

49 

42 

198 

50 

43 

199 

51 

44 

200 

52 

45 

201 

5  3 

46 

202 

54 

47 

203 

55 

48 

204 

57 

50 

205 

58 

51 

206 

59 

52 

207 

60 

5  3 

208 

61 

54 

209 

62 

55 

2  It 

63 

56 

MEMBER  PROPERTIES»PRISMATK 

1 

IY 

0.05151 

IX 

0.01913 

2 

IY 

0.05151 

IX 

.01913 

3 

IY 

0.05151 

IX 

.01913 

4 

IY 

0.05151 

IX 

0.01913 

5 

IY 

0.05151 

IX 

0.01913 

6 

IY 

0.05151 

IX 

0.01913 

7 

IY 

0.05151 

IX 

0.01913 

8 

IY 

0.103U2 

IX 

.03627 

9 

IY 

0.10302 

IX 

0.03627 

10 

IY 

0  .103uZ 

IX 

0*03827 

11 

IY 

0.103U2 

IX 

0.03O27 

12 

IY 

0.103U2 

IX 

0.03827 

13 

IY 

0.10302 

IX 

0.03027 

14 

IY 

U.  103u2 

IX 

0.036^7 

15 

IY 

0.103u2 

IX 

0.03827 

16 

IY 

0.10302 

IX 

0.U3827 

17 

IY 

0.1C3U2 

IX 

0.03027 

18 

IY 

0.10302 

IX 

.03827 

IV 

IY 

0.10302 

IX 

0.03o27 

20 

IY 

0.10302 

IX 

0.03827 

21 

IY 

0.103o2 

IX 

0.03627 

22 

IY 

0.10302 

IX 

0.03627 

23 

IY 

0.l03o2 

IX 

0.03627 

24 

IY 

0.103o2 

IX 

o. 036^7 

25 

IY 

u.  10302 

IX 

0.03027 

^  ■-..!. 
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26 

[  Y 

0.10302 

IX 

0.03827 

27 

[  Y 

.  103u2 

IX 

0.03627 

28 

I  Y 

0.10302 

IX 

0.03827 

29 

[  Y 

0.10302 

IX 

.03827 

30 

[  Y 

.103  02 

IX 

0.03827 

31 

[  Y 

C. 10302 

I  x 

0.03027 

32 

tY 

0.10302 

IX 

0.u3OZ7 

33 

[  Y 

C. 10302 

IX 

0.03827 

34 

IY 

.103 j2 

IX 

0.03827 

35 

[  Y 

0.103^2 

IX 

O.03b^7 

36 

I  Y 

0.10302 

IX 

0.0jo^7 

37 

[Y 

0. 10302 

IX 

0. 03o^7 

38 

[  Y 

.10302 

I  X 

0.03027 

39 

EY 

0.10302 

IX 

0.u3o27 

EY 

0.103  02 

IX 

0.03oz7 

41 

EY 

0.10302 

IX 

0.03827 

42  J 

[Y 

. L0302 

IX 

0.03827 

43  ] 

EY 

0.10302 

IX 

0.03827 

44 

[Y 

0.1U302 

IX 

O.O^b^  ( 

45  1 

EY 

0.10302 

IX 

0»03o^7 

46  J 

Y 

0.103U2 

IX 

0.03827 

47  J 

LY 

0.10302 

IX 

o.Oio^  / 

48  J 

Y 

0.10302 

IX 

0.03o^  I 

49  ] 

Y 

0.103o2 

IX 

0.03827 

50  ] 

Y 

0.05151 

IX 

0.01913 

51  ] 

[Y 

.05151 

IX 

0.01913 

52  ] 

Y 

1  .05151 

IX 

0.01913 

53  ] 

Y 

0.05151 

IX 

0.01913 

54  ] 

Y 

0.05151 

IX 

0.01913 

55  1 

.Y 

C. 05151 

IX 

0.01913 

56  J 

EY 

0.05151 

IX 

0.01913 

57  ] 

Y 

0.05151 

IX 

0.01913 

58  ] 

Y 

0.10302 

IX 

0.03027 

59  ] 

lY 

0.1o3u2 

IX 

0.03o^7 

60  j 

[  Y 

0.10302 

IX 

0.OJOZ7 

61  J 

EY 

0.10302 

IX 

0.030^7 

62  J 

Y 

0.10302 

IX 

0.03027 

63  J 

[Y 

.  10302 

IX 

0.03827 

64  ] 

Y 

.05151 

IX 

0. 01913 

65  1 

IY 

0.05151 

IX 

0.01913 

66  1 

EY 

0.10302 

IX 

0.03827 

67 

[  Y 

0.103u2 

IX 

0.03027 

68  . 

[  Y 

C. 10302 

IX 

0.03027 

69  1 

Y 

0. 10302 

IX 

0.03827 

70  . 

EY 

0.103U2 

IX 

0.03827 

71  i 

EY 

U.i03u2 

IX 

0.03827 

72  1 

EY 

0.05151 

IX 

0.01913 

73 

[  Y 

0.05151 

IX 

0.01913 

74 

[  Y 

0.103O2 

IX 

0.03027 

75 

EY 

0.103^2 

IX 

0.03827 

76 

[  Y 

0.103^2 

IX 

0.03827 

77 

EY 

0.10302 

IX 

0.03827 

78 

EY 

0.10302 

IX 

0.03O27 

79 

EY 

.10302 

IX 

0.03027 

60 

EY 

.05151 

IX 

0.01913 

« 
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8i   : 

[Y 

0.05151 

IX 

0.01913 

82  j 

[Y 

;  .  )  3  0  2 

IX 

0  .0382  7 

83  ] 

Y 

0.1 03 0 2 

IX 

0.03827 

84  ] 

Y 

.103  02 

IX 

0.0.5^7 

85 

IY 

,i^5v2 

IX 

0.03527 

86  ] 

[Y 

.I03o2 

IX 

0.03827 

87  ] 

Y 

0.10302 

IX 

.03827 

88  1 

[Y 

0.0  5.151 

IX 

.01913 

89  ] 

;y 

0.05151 

IX 

0.01913 

Y 

0.10302 

IX 

0.03827 

91  ] 

Y 

0.10302 

IX 

0.03527 

92  1 

Y 

.10302 

IX 

0.035,:  1 

93  J 

Y 

0. 10302 

IX 

0.036^7 

94  ] 

Y 

0.10302 

IX 

0.035<i7 

95  ] 

Y 

0.10302 

IX 

0.03827 

96  ] 

Y 

0.05151 

IX 

0.01913 

97  ] 

Y 

0.05151 

IX 

0.01913 

98  ] 

Y 

.10302 

IX 

0.03027 

99  ] 

Y 

0. 10302 

IX 

0.03527 

100  ] 

Y 

0.10302 

IX 

0.03527 

101  1 

Y 

,103o2 

IX 

0.03827 

102  1 

Y 

0.103w'2 

IX 

0.03827 

103  ] 

Y 

0,10302 

IX 

0.03527 

104  1 

Y 

0.05151 

IX 

0.01913 

105  1 

Y 

0,05151 

IX 

0.0191 5 

106  ] 

Y 

.  Iu3u2 

IX 

0.03827 

107  ] 

Y 

0.10302 

IX 

0.03527 

108  J 

Y 

o.io302 

IX 

0»u35^7 

109  j 

lY 

0,10302 

IX 

.03527 

110 

lY 

C. 10302 

IX 

0.035^7 

111  ' 

EY 

0.10302 

IX 

0.03527 

112 

t  Y 

0,05151 

IX 

0.01913 

113 

IY 

0.06244 

IX 

0.0 

114 

[  Y 

0.06244 

IX 

0.0 

115 

[  Y 

0.06244 

IX 

0.0 

116 

[  Y 

.06244 

IX 

.0 

117 

EY 

0.062^+4 

IX 

0.0 

118 

Y 

0.06244 

IX 

0.0 

119  . 

[  Y 

0.06244 

IX 

0.0 

120 

[Y 

0.06244 

IX 

0.0 

121 

tY 

0.06244 

IX 

.0 

122 

[  Y 

.062  44 

IX 

0.0 

123  ] 

Y 

0.06244 

IX 

0.0 

124 

t  Y 

0.06244 

IX 

o.o 

125 

[Y 

,06244 

IX 

.0 

126 

[  Y 

0.06244 

IX 

0.0 

127 

[  Y 

0.06244 

IX 

O.o 

128 

[Y 

0.06244 

IX 

0.0 

129 

[  Y 

0.06244 

IX 

0.0 

130 

[  Y 

0.06244 

IX 

0.0 

131 

[  Y 

0.062  44 

IX 

o.  o 

132 

[  Y 

C. 06244 

IX 

0.0 

133 

[  Y 

0.06244 

IX 

0.0 

134 

EY 

0.06244 

IX 

o.o 

135 

EY 

..6244 

IX 

0.0 
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136 

IY 

.06244 

IX 

0.0 

137 

I  Y 

0.06244 

IX 

0.0 

138 

I  Y 

0.06^ 

IX 

. 

139 

I  Y 

C  .  0  6 

IX 

o.o 

I  Y 

0.06244 

IX 

0.0 

141 

I  Y 

0.06244 

IX 

. 

142 

IY 

0.06244 

IX 

0.0 

143 

I  Y 

.06244 

IX 

0.0 

144 

I  Y 

. .062  44 

IX 

0.0 

145 

I  Y 

0*06244 

IX 

146 

IY 

.06244 

IX 

0.0 

147 

I  Y 

0.06244 

IX 

0.0 

148 

I  Y 

0.06244 

IX 

0.0 

149 

I  Y 

0.06244 

IX 

0.0 

15C 

I  Y 

.062  44 

IX 

0.0 

151 

IY 

0.U6244 

IX 

0.0 

152 

IY 

0.06244 

IX 

0.0 

153 

I  Y 

Q.U6^<+h 

IX 

. 

154 

I  Y 

u.u62^4 

IX 

0.0 

155 

I  Y 

0.06244 

IX 

0.0 

156 

I  Y 

0.06244 

IX 

0.0 

157 

1  Y 

0.06244 

IX 

0.0 

158 

I  Y 

.0  62  44 

IX 

0.0 

159 

I  Y 

0.06244 

IX 

0.0 

160 

I  Y 

0.06244 

IX 

0.0 

161 

IY 

0.06244 

IX 

o.o 

162 

I  Y 

,06244 

IX 

0.0 

163 

I  Y 

0.06244 

IX 

0.0 

164 

I  Y 

0.06244 

IX 

0.0 

165 

IY 

0.06244 

IX 

0.0 

166 

IY 

0.06244 

IX 

0.0 

167 

I  Y 

0.06244 

IX 

0.0 

168 

IY 

0.06244 

IX 

0.0 

169 

IY 

0.0624<+ 

IX 

0.0 

170 

IY 

C. 06244 

IX 

0.0 

171 

IY 

0.06244 

IX 

0.0 

172 

I  Y 

0.06244 

IX 

0.0 

173 

IY 

.06244 

IX 

CO 

174 

I  Y 

0.06244 

IX 

0.0 

175 

I  Y 

0.06244 

IX 

0.0 

176 

I  Y 

0.06244 

IX 

0.0 

177 

IY 

0.U6244 

IX 

. 

178 

IY 

0.06244 

IX 

0.0 

179 

I  Y 

0.06244 

IX 

0.0 

18C 

IY 

0.U6244 

IX 

0.0 

181 

IY 

0.06244 

IX 

0.0 

182 

I  Y 

0.06244 

IX 

0.0 

183 

I  Y 

.06244 

IX 

.  '■ 

184 

IY 

0.06244 

IX 

0.0 

185 

I  Y 

0.06244 

IX 

0.0 

186 

IY 

0.06244 

IX 

CO 

187 

IY 

0. 06244 

IX 

0.0 

188 

I  Y 

0.06244 

IX 

o.o 

189 

IY 

0.06244 

IX 

o.o 

190 

I  Y 

0.06244 

IX 

0.0 
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191 

IY 

.06244 

IX 

0.0 

192 

IY 

.0  62  44 

IX 

0.0 

193 

IY 

0.06244 

IX 

0.0 

194 

IY 

.06244 

IX 

. 

195 

IY 

0.062  44 

IX 

0.0 

196 

IY 

.06244 

IX 

0.0 

197 

IY 

0.0< 

S244 

IX 

0.0 

198 

IY 

J<4A 

IX 

. 

19V 

IY 

0.06244 

IX 

0.0 

IY 

.06244 

IX 

C  .0 

201 

IY 

0.06244 

IX 

0.0 

202 

IY 

.  J6244 

IX 

0.0 

203 

IY 

0.06244 

IX 

0.0 

204 

IY 

0. 06244 

IX 

0.0 

205 

IY 

C. 06244 

IX 

0.0 

206 

IY 

.06244 

IX 

0.0 

207 

IY 

0.062  44 

IX 

0.0 

208 

IY 

0.06244 

IX 

0.0 

209 

IY 

0.06244 

IX 

0.0 

IY 

.06244 

IX 

0.0 

LOADING  •] 

L« 

MEMBER  LOADS 

57 

FORCE 

Z 

LINEAR 

WA 

0.0 

WB 

-9, 

,  104  - 

LA 

o. 

Ld 

12  .057 

65 

FORCE 

z 

LINEAR 

WA 

-9.184 

Wb 

-18, 

,3  67  - 

LA 

0.  . 

LB 

12  .o57 

73 

FORCE 

z 

LINEAR 

WA 

-18.367 

wb 

-27, 

,551  - 

LA 

• 

Lb 

12.857 

81 

FORCE 

z 

LINEAR 

W  A 

-27.551 

-36, 

,7D3  - 

LA 

0.0 

LB 

12.857 

89 

FORCE 

z 

LINEAR 

WA 

-36.735 

wb 

-45, 

.918  - 

LA 

0. 

Lb 

12  .o37 

97 

FORCE 

z 

LINEAR 

WA 

-45.918 

a/B 

-35, 

,102  - 

LA 

0.0 

Ld 

12.857 

105 

FORCE 

z 

LINEAR 

W  A 

-55.1U2 

Wb 

-64, 

,2o6  - 

LA 

0.0 

Ld 

12 .857 

MEMBER  LOADS 

58 

FORCE 

z 

LINEAR 

WA 

0.0 

WB 

-18, 

,367  - 

LA 

C. 

lb 

12.857 

66 

FORCE 

z 

LINEAR 

WA 

-18.367 

WD 

-36, 

,7^5  - 

LA 

o. 

Lb 

12.857 

74 

FORCE 

z 

LINEAR 

WA 

-36.735 

WB 

-33, 

,102  - 

LA 

0.0 

LB 

12.837 

82 

FORCE 

z 

l INEAR 

WA 

-55.102 

W  B 

-7  3, 

,469  - 

LA 

0.0 

Lb 

12.857 

90 

FORCE 

z 

LINEAR 

W  A 

-73 .Hb9 

WlJ 

-91, 

.8  37  - 

LA 

o. 

Lb 

12.857 

98 

FORCE 

z 

LINEAR 

WA 

-91.8  37 

wb 

-110, 

,204  - 

LA 

0.0 

Lb 

12.857 

106 

FORCE 

z 

LINEAR 

WA 

-110.204 

-128, 

,371 

LA 

0.0 

Ld 

12.857 

MEMBER  LOADS 

59 

FORCE 

z 

LINEAR 

WA 

0.0 

WB 

-18, 

,367  - 

LA 

0.0 

LB 

12  .857 

67 

FORCE 

z 

LINEAR 

WA 

-18.367 

Wb 

-36, 

,735  - 

1 
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LA  0.0  Lb 

75  FORCt  Z  LINEAR 

LA  0.      Ld 

8  3  FORCE  Z  LINE 
LA  0.0    LB 

91  FORCE  Z  LINEAR  WA 
LA  0.     Ld 

99  FORCE  Z  LINEAR  WA 
LA  0.0    LB 

107  FORCE  Z  LINEAR  WA 
LA  0.0  LB 
MEMBER  LOADS 

60  FORCE  Z  LINEAR  WA 
LA  0.0    Ld 

68  FORCE  Z  LINEAR  WA 
LA  0.0    LB 

76  FORCE  Z  LINEAR  WA 
LA  0.0    Ld 

84  FORCE  Z  LINEAR  WA 
LA  0.0    Lb 

92  FORCE  Z  LINEAR  WA 
LA  0.0  LB 
10U  FORCE  Z  LINEAR  WA 
LA             0.0    LB 

108  FORCE  Z  LINEAR  WA 
LA  0.0  Lb 
MEMBER  LOADS 

61  FORCE  Z  LINEAR  WA 
LA  0.0    LB 

69  FORCE  Z  LINEAR  WA 
LA  0.0    LB 

77  FORCE  Z  LINEAR  WA 
LA  0.0    LB 

85  FORCE  Z  LINEAR  WA 
LA  0.0    Lb 

93  FORCE  Z  LINEAR  WA 
LA  0.0    Lb 

101  FORCE  Z  LINEAR  WA 
LA  0.0    Ld 

109  FORCE  Z  LINEAR  WA 
LA  0.0  LB 
MEMBER  LOADS 

62  FORCE  Z  LINEAR  WA 
LA  0.0    LB 

7  FORCE  Z  LINEAR  WA 
LA  O.v,    Ld 

78  FORCE  Z  LINEAR  WA 
LA  0.0    LB 

86  FORCE  Z  LINEAR  WA 
LA  0.0    LD 

94  FORCE  Z  LINEAR  WA 
LA  d .  C    L  3 

102  FORCE  Z  LINEAR  WA 
LA  0.      Lb 

110  FORCE  Z  LINEAR  WA 


1^.857 

-36. 735 

-55.102  - 

12. 

-55.102" 

WB 

-73.469  - 

12.857 

-73.469 

WB 

-91.837  - 

12.657 

-91.837 

-110.204  - 

12.657 

110.204 

WB 

-128.571  - 

12.857 

0.0 

WB 

-18.367  - 

12  .857 

-18.367 

WD 

-36.735  - 

12.857 

-36.735 

WB 

-55.102  - 

12 .637 

-55*102 

WB 

-73.469  - 

12.657 

-73.469 

WB 

-91.837  - 

12.857 

-91.837 

WB 

-110.204  - 

12.857 

110.204 

WB 

-128.571  - 

12  .057 

0.0 

WB 

-18.367  - 

12  .837 

-18.367 

WB 

-36.733  - 

12.857 

-36.735 

-33.102  - 

12  .857 

-55.102 

Ift  B 

-73.469  - 

12  .037 

-7  3.')  69 

Wb 

-91.637  - 

12.857 

-91.037 

Wb 

-110.204  - 

12  .637 

110.^04 

Wb 

-128.371  - 

12.857 

0.0 

Wb 

-18.367  - 

12  .857 

-18.367 

Wb 

-36.733  - 

12.657 

-36.735 

WB 

-55.102  - 

12.857 

-55*102 

WB 

-73.469  - 

12.637 

-73.469 

WB 

-VI. 8 37  - 

12.857 

-91.837 

Wb 

-110.204  - 

12  .637 

110.^04 

-126.371  - 

• 
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- 


«?• 


LA 

0.0 

LB 

12.857 

MEMBER  LOADS 

63 

FORCE 

Z 

LINEAR 

WA 

0.0 

LA 

o.o 

Ld 

12.857 

71 

FORCE 

Z 

LINEAR 

WA 

-18.367 

LA 

0.0 

LB 

12.857 

79 

FORCE 

Z 

LINEAR 

WA 

-36.735 

LA 

0.0 

LB 

12.857 

87 

FORCE 

z 

LINEAR 

WA 

-55.102 

LA 

0.0 

LB 

12.857 

95 

FORCE 

7 

LINEAR 

WA 

-73.469 

LA 

0.0 

LB 

12.857 

103 

FORCE 

z 

LINEAR 

WA 

-91.837 

LA 

0.0 

LB 

12.857 

111 

FORCE 

z 

LINEAR 

WA 

-110.204 

LA 

0.0 

LB 

12.857 

MEMBER  LOADS 

64 

FORCE 

z 

LINEAR 

WA 

0.0 

LA 

0.0 

LB 

12.857 

72 

FORCE 

z 

LINEAR 

WA 

-9.184 

LA 

0.0 

LB 

12.857 

80 

FORCE 

z 

LINEAR 

WA 

-18.367 

LA 

0.0 

LB 

12.857 

88 

FORCE 

z 

LINEAR 

WA 

-27.551 

LA 

0.0 

LB 

12.857 

96 

FORCE 

z 

LINEAR 

WA 

-36.735 

LA 

0.0 

LB 

12.857 

104 

FORCE 

z 

LINEAR 

WA 

-45.918 

LA 

0.0 

LB 

12.857 

112 

FORCE 

z 

LINEAR 

WA 

-55.102 

LA 

0.0 

LB 

12.857 

LOADING  LIST  ' 1 ' 

STIFFNESS 

ANALYSIS 

LIST 

REACTIONS 

LIST 

LOADS 

10 

LIST 

LOADS 

19 

LIST 

LOADS 

20 

LIST 

LOADS 

LIST 

'  DISPLACEMENTS 

42 

LIST 

DISPLACEMENTS 

42 

LIST 

DISPLACEMENTS 

51 

LIST 

DISPLACEMENTS 

FINISH 

WB 


WB 


WB 


WB 


WB 


WB 


WB 


WB 


WB 


WB 


WB 


WB 


WB 


WB 


-18.367 
-J6.735 
-55.102 
-73.469 
-91.837 
•110.204 
■128.571 

-9.184 
-18.367 
-27.551 
-36.735 
-45.918 
-55.102 
-64.286 
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Off,  i-R0F9S^  IJSTJW 


/JOB  1673   KOLLINIATIS  1^9. 

iSIGN   3=TCES 


/RESTORE  CDP1 
CDP 


INITIATE  • PLATES*  » SOLUT TON  • 


ADD  •SPAN'  REAL  »HOLEN»  REQUIRED 
FILE 


ADD  'THICK •  REAL  'THICK*  REQUIRED 
PILE 


-•"YOUNG*  REAL  *YOUNG^  ST-ANDARD~3iK>0O. 


RE  Atr  +PRA  T I/O-*- STANDARD  >3 
-pftE 


^ADD  -*MAX»  INTEGER  *MANUJO*  STANDARD  300 
FILE 


'AT  X*  REAL  *H0LEN1*  REQUIRED 


TO  *AND  Y*  REAL  *VELEN1 •  REQUIRED 


-1-0— •-SPAN*  REAL  *HOLEN?«  ^EOUIRED 


-i-D-'WID-TH*  REAL  *VELEN2  •  REQUIRED 


-EXECUTE  ^DIPOHO»  -»-MUGRT-D-«  » METNHO+ 
FILE 


ADD  'POINT  S' 


REPEAT  TABULAR 


-FD 'X'  REAL  'EX'  REQUIRED 


•Y'  REAL  'WHY*  REQUIRED 


T 

■ 

I 

f  ,     ,,      ■  ,••.-■ 

'  •  V\  I     •  X  AM  * 

i  ■-<  •  • 

i    •  .  , 

,  .    i  .    i         >    . 

JH 

T 

■■.;      i   .       . 


EXISTENCE  *FREE»  'HINGED*  'FIXED*  SET  »KIND1*  STANDARD  0 


EXECUTE  *POBOC* 


END  REPEAT 


FILE 


ADD  'LINE  S* 


REPEAT  TABULAR 


ID  *X1*  REAL  *EX*  REQUIRED 


ID  *Y1*  REAL  'WHY*  REQUIRED 


ID  *X2*  REAL  *EX2*  REQUIRED 


ID  *Y2*  REAL  *WHY2*  REQUIRED 


EX  I  STENCE  •  FREE*~*HTnGED*  •  FIXED"*~3ET"«KINDr»_ STANDARD  0 
EXECUTE  *  L  I  BOO- 
"END  REPEAT" 

FT  LI- 


ADD  'FORM  G* 


EXECUTE  *JOINRE*  'MEMPRO* 


FILE 


iDD  •POINT  L* 


PRESET  INTEGER  'KIOUT*  EQUAL  1 


PRESET  INTEGER  «NUPOLO*  EQUAL  1 
"REPEAT  TABULAR 
EXISTENCE  'NOTHING*  *Z  FORCE1  *X  MOMENT*  »Y  MOMENT*  »SET  *NLOTY*  - 


REQUIRED 


ID 

•X* 

REAL 

•EX* 

REQUIRED 

ID 

•Y* 

REAL 

•WHY 

»  REQUIRED 

ID 

.pi 

REAL 

•STALD     QUI 

EXECUTE  'POINLO' 


tnd~"r"epeat 
ttlf; 


try, 

|l  !  » 

1 

I 
,  |  f  '  " 

_         ,  .         ■  •  II  '  •  I  •  I 

«  >  .  t     ,     » 

'  q  •    1 1 

T  A  ' 


151. 

PR E-SET- INTEGER  'KI0UT* EQUAL  1 
REPEAT  TABULAR 


EXISTENCE  'FORCE  Z  LINEAR*  'FORCE  Z  UNIFORM*  'MOMENT  X  LINEAR'- 
-'FOMENT— ^UNIFORM'  'MOMENT  Y  LINEAR'  'MOMENT  Y  LINEAR' 


'MOMENT— X— UNIFORM'  •  MOMENT  JY—4H NEAR  •  'MOMENT  Y  UNIFORM'SET  HCFbTrtO^ 

IREO 

ID     'AT    X'    REAL     'EX'    REQUIRE-D— 
ID     'AND^M   4R  EAL     'WHY'    REQU  I  RE  D 


4T>    'UP    TO    Xt-  REAL     *EX2<    REQUIRED 


ID     'AND    UP    TO    Y'    REAL     'WHY2'    REQUIRED 


^OND  TXTON    •  K  Tt-TtO  *  ^NT^GE  R-EQUA  L 
_P_R_E-SEI-.JNT  EGER     ♦  NLQTXi-  EQUAU-4— 


PRESET    INTEGER     'LTYPE'EQUAL    1 


■EXECUTE    UXNEL-CL' 


OR    C0NO1TI  ON    '  K I L IUX*    I  NT  EGER-  -E  QUAL    1 


-PRESET    INTEGER     *NLOTY»    EQUAL    1 


-PRE SEX  4  N TEGER^LXYPE'  EQUAL    0 
-EXECUTE  -M.JNELCM 


OR    CONDITION    'KILILO'    INTEGER    EQUAL    2 


iSEZ- XNT EGER     •  NLO TYJL  ^QUAL    2 


JPRESEX   INTEGER     •  LTYPE'  EQUAL_X 
EXECUTE  -LLLMEi-OJ 


J^Et  COND I  T  I  ON  _UU  L  I LO  !_  I  N T EGER-  EQU  A4__3- 
PRE SET    INTEGER     'NLOXY'    EQUAL— 2— 


PRESET    INXEGER    'LT-YP-E-'-EQUAL    0 


•LINELO' 


OR  rONDJTION    'KILILO'     INTEGER    EQUAL    4 
i>RE£EX-IlNTi£GER-    'NLOTY'     EQUALS 


preset  integer  •  ltype 'equai l 

execute:  'lineloj 
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-OT-WER-W-ISE 


PRESET    INTEGER     'NLOTY'     EQUAL-3^ 


PRESET    INTEGER     'LTYPE'    EQUALS 
-EXECUTE    LINE  1.0 


END    REPEAT 


LLE 


ADD,  'AREA— L-i- 


RRESET    INTEGER     'KTOUT'    EQUAL.  4- 
REPEAT    TABULAR 


EXISTENCE    'FORCE    7    UNIFORM'     'FORCE    7    LINEAR    ON    X*     'FORCE    Z    LINEAR    ON   Y» 
_'_MOMENT    X  UNI FORM  •     •  MOM FNT    X    LTNFAR    ON    X •    AMQMENT^X   LINEAR jQN_  _Y. J^_ 
JJ4QMENT    Y    UNIFORM'     'MOMENT    Y    1  INEAR    ON   X«     'MOMENT    Y    IINEAR-ON    Y  « -  - 

SET     'KT1.ILD1    RFQUTRFD 


XD— 1XAX  REAL    'EX'    REQUIRED 


I D     '  YA-L_  RE  AJ LW  H  Y  •    R  E  QU I R  ED 


ID     «WA«    REAL    J-L-DMSTAJ-REQUXRED 


JJX__' XRL  REAL    •  EX2  •    REQU  LREa 


XCL^YB'    REAL-  «WMY_2L-REQUJ RED 


ID    'WB'-REAL— LLDMENO'    REQUIRED^ 


^ONT-I-T-T0N^-MCT-L4L0'     INTEGER 
PRESET    INT-EGER   J  WLOTYJ-EQUAL-i- 


-PRESET    INTEGER    »LTYPE-'  EQUAL    0 


-EXE-CUTE  -*  A^EAta  l 


SM-C  OND 1 T-I  ON    •  K  I  L  IX.Q -•-  I NJ-E-6ER-  E^UJAJ 1- 

PRESET  -INTEGER -*MJ3T-Y<    EQUAL    1 


PRESET    I NEGER     'LTYPE'    EQUAL— L- 
EXECUTE     'AREALO' 


lLT-I  ON  -  'X1L-I-L0 -t  J  NTEG  ER    E QU  Al  .£. 


'RESET    INTEGER     'NLOTY'     EQUALS 


PRESET    INTEGER     'LTYPE*     EQUAL    2^ 


■ 
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I     I 

i  X  3 ' 

i  •     . 

•     ■.     i  •        - 

■  -    .  ,  •        i 

-     ,     ■ 

I     • 
,  -     i 


-E-XE-CUTE  'AREALO1- 


OR  CONDITION  'KILILO'  INTEGER  EQUAL-^ 
PRESET  INTEGER  'NLOTY'  EQUAL  2 


PRESET  INTEGER— U.TYPE '--EQUAL  0 


EXECUTE  'AREALO' 


CONDITION  'KILILO'  INTEGER  EQUAJ V 

PRE-S E T -I  N T EGER- *-KLO-T-Y-i— EQUAL  2 


PRESET  INTEGER  'LTYPEJ— EQUAL  L 
^EXECUTE  »AREALQA_ 


_QR__CQNDI  T  ION  'KILILO'  INTEGER  EQUAL  5 


IESET  INTEGER  •  NL-QTXL-EQUAL-^2- 


'reset  integer   'ltype'  equals 
;auxi 


OR    COND I  T I  ON    •  K IL I  LCLl    I NTEGER-EQUAL-4- 
-PRESET  J NT EGER    'NLOtY^    EQUAL  3^ 


PRESET  INTEGER  'LTYPE'  EQUAL  0 
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XXE-CU TE-    'ARE A L 0  _• 


QR    CQMD I  T I  ON    •  K I  LlLCLi^XfcLTEGER  EQUAI Z_ 

_ ..  PRE  SET  -  IN-T  EGER  .'NLQJXL,  EQUAL    3 


PRESET    INTEGER    »LTYPE*    EQUAI t- 


EXECUIE_  •AREALO' 


OR    rONDTT  LON    >  K  1 1. 1 1  0  '     IN  T  EGER    EQUAI     8 


PRESET    INTEGER    'NLOTY'    EQUAI 3- 


_RRESEI-1N  T  EGER    'LTYP E  •    EQ U  AL—2- 
FXFrilTF     «ARFAIQ» 


END    REPEAT 


£LLE_ 


ADD    'LIST' 


REPEAT    TABULAR 


MODIFIER    'REACTIONS* 


■ 
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• 
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PRESET  INTEGER  «KIOUT-»  -EQUAL  2 


EXECUTE  •OUT* 


OR  MODIFIER  'POINT  LOADS* 


XD-i-XA-'-  REAL  «EX»  REQUIRED 


ID  »YA«  REAL  'WHY'  REQUIRED 


PRESET  INTEGER  'KIOUT'  EQUAL  3 
-EXECUTE— IQUTf 


OR  MODIFIER  'LINE  LOADS' 


JJ1    LXA'  REAL  'EX'  REQUIRED 


ID  «YA»  RFAJ 'WHY*  RFQUTR  ED 


iCL  ±XB'  REAL  !EX2'  REQUIRED 


ID  «YB«  REAL  'WHY2'  REQUIRED 


PRESET  INTEGER  'KIOUT'  FQUAL  4 


if*. 


ITE    'OUT' 


J3R„M_ODIFIER    'LOADS' 


_PRE_SET„  J N TEGER     'KTflUT  L  EQUAL    5 
EX  E  CUT  E_JilUU 


OR   MODIFIER    'POINT    DISPLACEMENTS! 
ID     'XA'  -REAL  _'FX'    REQUIRED 


n     «YA»     RFAI      » WHY  •    REQUIRED 


PRES£T,XNTEGER  'KIOUT'  FQUAL  6 
EXECUTE  'OUT* 


jQRJYLQIUFIER  'LINE  DISPLACEMENTS! 
ID  'XA«  REAI 'EX*  R EQU I R ED_ 


in  »YA*  RFAI   tWHYt  RFQUIRFD 


.ID— LXB-L  REAL  «EX2'  REQUIRED 


ID  «YB'  REAL  'WHY2'  REQUIRED 


PRESET  INTEGER  'KIOUT'  EQUAL  7 


:UTE  'OUT' 


OR  MODIFIER  'DISPLACEMENTS' 
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WE SET    I NTEGER    'KIOUT *-  EQUAL    8 


rXECUTE    •OUT* 


END    MODIFIER    REQUIRED 


END    REPEAT 


"ILE 


MNISH 


J-EUD   OF-  S1UL 


1 
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AFrSSDIX     E 
A  D1CTI03ARX  OF  TSR/3S  IN  |  »S30N 

JBH2LO  «     End  load  Magnitude  of  line  loads 

■        *    x-coordinate  of  point  supports 

EXZ      *     2nd  x-coordinate  of  line  supports 

hOCOR  *     x-coordinate  of  joint 

H0CQR1  •  x-coordinate  of  a  non  support  joint 

tiOtm  *    Span  of  plating  (parallel  to  the  x-axis) 

ttULEivl  »  Distance  along  the  x-axis  from  the  center  of  coordinates 

up  to  the  left  edge  of  hole 
xiOLEH2  ■  Span  of  hole  (parallel  to  the  x-axis) 
HOLES®  m  HOLEN  -  ttOMCNX  -  fiGL£N2 
JORENU  ■  Joint  release  number 
K1K&I0  *  Kind  of  joint 
UNQKS  ■  Kind  of  member 

KJhXST    »  A  number  that  specifies  necessary  output 
JUJUOi     =  Kind  of  point  support 
LliHSTA  =  Start  load  magnitude  for  area  loads 
UXi&ftD  *  End  load  magnitude  for  area  loads 
LTXPE    ■  A  number  that  specifies  kind  of  load 
HAlflJJG  »  Maximum  number  of  joints 
KJ  P    umber  of  joints 

aUJO    *     Identification  number  of  joint 
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NUJ01     ■  Identification  number  of  non-support  joint 

tiULDS    «  dumber  of  grid  loads  used  to  approximate  a  given  LIKE  - 

I    AD 

idJMKMB  ■  Identification  number  of  member 

KUMIJO  »  Identification  number  of  sinus  joint 

MfHI      *  'lumber  of  members 

MJ1-UC  =  Identification  number  of  plus  joint 

_is!DI  *  'lumber  of  left  diagonal  members 

SIORIDI  *  dumber  of  right  diagonal  members 

MHOMB  ■  Mimber  of  horizontal  members 

iJUVSMB  *  Number  of  vertical  members 
Wflfl 

■fH)        Momber  of  division  points  exjraained  in  DIPGBO 

jmoji 

MM 

HhOTL  -     dumber  indicating  kind  of  load 
PRaTIO  *  Poisson#s  ratio 

.  m  rcint  load  magnitude 
TalCX  »    Thickness 
VBCQR  *    y-coord.  of  joint 
VBCORl  •  y-coord.  of  non  suport  joint 
WM  ■    Width  of  plate 
VJGLKK1  *  y-coordinate  of  lower  edge  of  hole 
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VELE82  »  Width  of  hole 

EXQ  «  HUH  -  f&fift  -  ^BLEN2 

ittg         a  y-coordinatc  of  a  point  support 

WHX2     =    y-end  coordinate  of  a  line  supi«ort 

XQiii.G  ■     Iouag,«>  raodulua  of  elasticity 

SLILZU) 

*  Kxtra  variables  not  used  by  programs 
fcUPCLG 
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1  -  b*W  ^  pwi'r*®  ,frt  %%  to  w* 

The  following  instructions  are  particularly  useful  la  using 
the  ■,  rc:;r&si  as  it  is  now.  The  input  of  course  lias  1  iven  ac- 

cording to  FC.1''ATS  already  included  In  the  projran.  Integers  have 
to  be  light  justified  in  their  field.  Floating  point  to/bet*  can 
be  i  iti  In  their  field  as  long  as  the  point  is  punched,  The 

units  to  be  used  are  a  d  .'3, 

Necessary  cards  are* 
ford  1  fopst  (ff  ^ 

In  the  first  10  columns  pBt  the  span  of  the  plate. 
In  columns  11-20  put  the  width  of  the  plate. 
In  colusms  21-30  put  the  thickness  of  the  plate. 
In  columns  31-^0  put  Xoung's  todulus. 
In  columns  41-50  put  roisson^s  ratio. 

<frra  z  tmmt  MB 

In  columns  1«*6,  rirfit  juatiflad    put  the  desired  rmxirain  ;msnbor 
of  joints. 

In  eoluran  4  put: 

1  if  there  is  no  hole 

2  if  there  is  a  hole 
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In  case  there  is  a  hole  additional  information  is  necessary  in 

this  card,   i  Cy| 

Columns  5-1^  Distance  along  span  frost  the  left  edge  of  the  plat© 

to  the  left  ©djre  of  hole. 

Columns  15-24  Distance  along  width  froia  the  lower  edge  of  the 

plat«  to  the  lower  edge  cf  hole. 

Columns  25*3%  Isnftti  of  r-ole  (in  the  direction  of  the  span). 

Columns  35-^  V*idth  of  hole  (in  the  direction  of  the  width  cf  the 

plate). 

Cart  *    ^^  (2Ifl 

In  columns  1-5  rlnht  .fo  stifled*  pat  the  number  of  point  supports. 

In  aolu&ns  6-10  rir.ht  justified,  put  the  number  of  line  supports* 

ftttti  .?  form*  Wt  ^Qi,3> 

lliere  nast  be  as  many  such  cards  as  there  are  point  supports. 
In  column  ^  put  as 

1  if  support  is  hinged 

2  if  support  is  fixed 

In  columns  5~i**  pat  the  x-coordinate  of  partieular  point  support, 

(See  Fig.  4.1)  (fop  location  of  axes) 

In  columns  5-1^  put  the  y-coordinate  of  partieular  point  support, 

cm  |  |— |  ttsffcfll 

There  nast  be  as  esany  such  cards  as  there  are  line  supports* 

In  column  4  put  a: 

1  if  support  is  hinged 
Z  if  support  is  fixed 
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Xr.  -         at  5*&4  PW*  tlMI  x«coordinate  of  starting  point  of  a  lint 
support,   (In  all  the  subsequent  description  the  starting  point 
of  a  line  is  the  one  to  the  left  if  the  line  is  not  parallel  to 
the  y-axis,  or  the  one  downwards  if  the  line  is  parallel  to  the 

y-ivxis) 

In  c  .     .  it  I5«&fc  put  the  y-ccordinato  of  the  stari.aj  point  of  I 
line  support. 

In  col  5-3'*  put  the  x-coerdinaie  of  the  and  point  of  the  li- 

ne support. 

In  colurms  35-^  put  the  ir-c  ordinate  of  the  end  point  of  the  line 
support. 

In  columns  1-4  put  the  number  of  point  loads,  right  .justified. 
In  columns  5-3  *»*  the  number  of  line  loads,  tftfft  .^UAfifii' 
In  columns  9-12  put  the  number  of  area  loads,  right  justified, 

girt  ?  tm*  u%mM) 

In  coluren  5  put  a: 

X  if  load  is  a  a  force 

2  if  load  is  an  x  moment 

3  if  load  is  a  y  pioaent 

(The  system  of  coordinates  is  assumed  righWuanded  and  the  loads 
are  positive  if  their  vector  points  towards  the  positive  direction 
of  the  axes) 

In  columns  6-15  put  the  x-ccordinate  of  the  point  of  application 
of  the  load. 
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the  ^-coordinate  of  the  point  of  a,,. plication 
of  the  load. 

In  coIufjis  26-35  F**t  the  nagr&tude  of  the  load. 
'in  ere  are  as  ■amy  such  cards  as  there  are  point  loads* 
Car?  |   vornat  (Ig.ffrMi  l§a  ^10,3} 

There  are  as  ssany  such  cards  as  there  are  line  loads. 
In  column  5  put  a  nontber  (1#2  or  3}  with  the  same  neardng  as  in 
previous  card. 

In  coluans  6-\5  pi%  the  x-coordinate  of  the  starting  point  of  the 
line  load. 

In  colu.rjis  16-25  pot  the  y-coerdinate  of  the  starting  point  of 
the  line  load. 

In  toilmmi  26-35  $**  the  ro&gnitud©  of  the  load  at  the  starting 
point. 
In  column  kQ  put  a: 

0  if  line  load  is  uniform 

1  if  line  load  is  linear 

In  columns  41-50  put  the  x-coordinate  of  the  end-point  of  the 

line  load. 

In  colusms  51-&0  put  the  y •coordinate  of  the  end-point  of  the 

line  load. 

In  columns  6l-?0  put  the  magnitude  of  the  load  at  the  end  point. 

Card  7  Format  (15.  3F10.3.  15.  3F10.3) 

There  are  as  nany  such  cards  as  there  are  area  loads. 
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bar  (1,2  or  3)  with  the  same  cleaning  as  in 
jorevivras  card. 

In  felejmi  6-15  P»t  th©  x  coordinate  of  the  left  lower  corner  of 
the  aroa  load. 

-x*  coiunns  i6-l":j>  pat  th©  y  coordinate  of  the  loft  lower  corner 
of  the  area  load. 

In  aolaunt  26-35  put  the  Magnitude  of  the  load  at  the  starting 
poi 
In  telann  ;f0  pst  a: 

0  if  load  It  uniform 

1  if  load  Is  linear  along  th©  x-direction 
Z  if  load  is  linear  along  th©  y-direction 

In  columns  ^1-50  pit  the  x-coordinate  of  the  uppar  rij^it  corner. 
In  oolnoai  51-60  put  the  y-coordinate  of  the  upper  right  corner. 
In  <ttltwtt  61-70  put  the  raagnitudo  of  the  load  at  the  end  point. 
Card  6  ioraat  (21h) 

In  columns  1-4  r^ht    just^fie^    put  the  number  of  different  kinds 
of  output  desired, 

fort  ?   ?9m\  US  -  Wt3) 

There  are  as  many  such  cards  as  there  are  kinds  of  output 
desired. 
In  eolunn  5  pnt  as 

Z  if  all  reactions  are  desired 

3  if  all  loads  are  desired  at  a  particular  point 
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f  the  loads  are  assirod  along  a  line 

5  if  all  loads  ar«  desired 

6  if  ents  aro  desired  at  a  particular  ;-olnt 
?  if  disj.lacew©nts  are  desired  along  a  line 

8  if  all  di«  oents  are  desired. 

In  case  the  above  up  is  a  2,5  or  8  no  other  information 

is  necessary  on  this  card. 

In  case  the  above  number  is  a  3  or  a  6: 
In.  tolurans  6*1  J  pat  the  x-coorrHnate  of  the  point  where  output 
is  desired. 

In  ooltmna  1'-         it  the  y -coordinate  of  the  point  where  output 
is  desired. 

In  case  the  above  number  is  a  f4  or  a  ?: 
In  columns  6-15  put  the  x-coord lnate  of  the  starting  point  of  the 
line  where  output  is  desired. 

In  columns  1  Rt  the  y-coordinato  of  the  starting  point  of  the 

line  i?here  output  is  desired. 

In  columns  26-35  put  the  x-coordinate  of  the  end  point  of 
the  line  where  output  is  desired. 

In  0olaaM  36-45  pit  the  y-eoordinate  of  the  end  point  of  the  line 
where  ouUtut  is  desired. 
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s  a  proposed  user's  formulation  of  *  prcblerat 
if  the  fc-ras  i  .  .xlicated  in     Appendix  Dt  is  used, 

Tha  statements  have  to  bo  in  the  order  indicated.  She  under* 
.i«d  part  of  a  statement  iradicatos  that  this  is  the  way  this  state* 
raeut  is  stored  la  the  dictionary.  Dfttft  identifiers  can  be  omitted 
if  I  ta  follow  in  the  specified  ordere 

qoawnda  ca  .ar  one  after  tat  ether, 

-  ■  *-<•■■•  •  •  • 
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optional 

one  or  the 
other 

Free  can  be 
omitted 
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FREE  ^ 
hINGJiB  ,> 
FIXED 
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